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A  B  S  TRA  C  T 


The  feasibility  of  accomplishing  high  frequency  direction 
finding  against  short  duration  (l'dZ-1 ZZ'o  ms)  H?  s.ryvave 
signals  using  narrow  aperture  antennas  is  investigated.  Two 
statistical  nrocedures  for  estimating  the  signal  hearing  are 
proposed  and  compared.  Tnese  procedures  employ  time  averaging 
to  reduce  the  large  instantaneous  hearing  error  caused  by  t.ne 
phase  and  amplitude  distortion  of  t.ne  wavefront  due  to 
scattering  and  multipath  interference.  Results  are  presented 
using  data  collected  with  the  Southwest  Research  Institute 
Coaxial  Spaced  Loop  i?D7  system.  It  is  shown  that  for  a 
limited  sample  of  data  f non  this  system  the  standard 
deviation  of  the  bearing  estimate  for  a  222  ms  signal  varied 
from  15  to  59  degrees. 
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I.  IMTPODUCTIOM 


A.  HISTORY 

High  frequency  radio  direction  finding,  abbreviated  in 
this  report  as  HFDF,  has  been  a  topic  of  interest  since  the 
first  uses  of  radio.  The  science  of  HFDF  has  developed  in 
spurts  of  ingenuity  and  need.  Tae  advances  in  electronic 
devices  and  the  more  accurate  modeling  of  HF  propagation  have 
been  important  steps  to  developing  accurate  HFDF  systems. 
However,  it  has  been  military  necessity  that  nas  spurred  the 
most  important  developments  in  this  field.  The  greatest 
concentration  of  published  literature  on  tnis  subject  resides 
in  the  technical  reports  published  during  and  in  the  decaae 
after  .vorld  War  Two. 

The  major  stuii°s  of  HFDF  have  dealt  with  the  problems  of 
polarization  of  skywave  signals,  the  effects  of  multipath 
propagation,  the  statistics  of  HF  propagation  through  the 
ionosphere  and  the  development  of  HFDF  antennas  and  arrays  of 
antennas.  An  important  distinction  has  developed  from  these 


studies  . 

There  are 

thre 

e  types  of 

*t  iT\-n 

h  j  Dj; 

an t enr.a s  :  ( 1  ,  wire 

aperture 

f  2  'i  medium 

aper  ture 

and 

(3'  narrow  aperture 

antennas 

.  .vhether 

an 

antenna 

array 

cc^uares  pna s  e  or 

amplitude.  the  primary  measure  of  problems  attendant  tc  its 
accurate  operation  is  the  width  of  its  aperture.  T.ne  term 
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aperture  in  this  paper  refers  to  the  linear  spatial  extent  of 
an  antenna,  not  to  an  area.  The  unit  of  measurement  is  either 
meter  or  wavelength.  If  the  aperture  is  on  the  order  of  one 
quarter  or  less  of  a  wavelength,  it  oar.  he  considered  narrow 
aperture,  and  it  suffers  the  greatest  number  of  difficulties 
to  achieving  accurate  direction  finding  capability. 

During  World  War  Two  the  Allies  experienced  considerable 
success  with  the  landbased  medium  and  wide  aperture  systems 
and.  not  unexpectedly,  limited  success  with  narrow  aperture, 
shipboard  K7DF  systems.  One  of  the  major  targets  of  the 
shipboard  H.FDF  systems  was  German  submarines.  As  improvements 
were  made  to  all  types  of  E7D7  systems,  the  submarine's 
transmissions  became  more  and  more  vulnerable.  In  an  effort 
tc  maintain  communications  and  to  tnwart  HFDF  systems,  toe 
Germans  shortened  the  duration  of  transmissions  to  lower  tne 
probability  that  the  transmissions  would  be  interceptec  and 
subsequently  located  by  ”7 07 .  A  highly  effective  means  of 
shortening  transmission  time  was  to  record  the  information  on 
tape  and  then  to  transmit  via  the  radio  at  a  much  faster 
playback  speed.  When  this  method  was  coupled  with  the 
practice  of  economizing  on  the  amount  of  information  sent, 
signal  durations  were  shortened  by  more  than  an  order  of 
magnitude.  A  rJ-bcat  employing  such  measures  was  appreciably 
less  susceptible  to  H7D7. 

The  problem  of  locating  a  short  duration  signal  remains 
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today.  It  is  still  a  common  problem,  even  when  the  target 
transmitter  does  rot  attempt  to  compress  its  signal.  In  a 
tactical  situation  it  is  typical  that  the  communication  net 
control  station.  usually  co-located  with  tne 
of f icer-in-tactical-conmand,  will  act  as  a  broadcasting 
station.  and  the  outstations  will  not  transmit  or  will  only 
transmit  a  b^ief  signal.  In  the  case  of  manual  morse  this 
signal  may  be  an  "T'"  for  "rcger  voim  last  transmission ”  .  cr 
in  the  case  of  tactical  voice  c omnunica t i c”s  tne  outstati'n 
will  brief lv  key  the  microphone.  In  eitner  case  tne 
transmission  may  not  last  longe"  than  200  to  ~Z 2 
milliseconds . 

The  rapid  growth  of  digital  communications  has 
significantly  increased  tne  ease  with  which  a  burst  signal 
can  be  generated  and  ’’reliably’’  received.  Given  a  digital 
pulse  of  duration  "t"  and  a  total  signal  duration  of  ”T  . 
there  is  a  sinole  exoression  for  tne  amount  of  information  in 
bits  that  can  oe  transmitted. 


If 


Figure  1 


the 


Partial  short  duration 
number  of  tits  is 


signal  bit  stream 
"n”,  then  N=?/2t. 


2 1  is 


lb 


approx ima telv  the  reciprocal  of  the  bandwidth;  tnerefore, 
N=(T)(3«).  approxima tely .  In  the  high  frequency  range  a 
bandwidth  cf  1?  .-cHz  can  be  readily  achieved.  As  an  example, 
if  BW=lf*  khz  anc  T=503  ms ,  M  equals  5300  bits.  Five  tnousand 
bits  is  sufficient  to  provide  considerable  encoded  tactical 
information.  Ever.  if  tne  signal  duration  were  limited  to  212 
ms  and  the  fi^st  half  cf  the  signal  dedicated  to  alerting  tne 
destination  receiver,  there  would  remain  1220  bits  for 
information.  At  five  bits  per  symbol  and  an  average  of  five 
symbols  peu  word.  this  would  allow  forty  words  to  be 
communicated  in  the  soace  of  200  its.  For  a  shin  on  a  covert 
thirty  day  patrol  sending  a  single  daily  status  report  via 
burst  communications.  the  total  communications  transmission 
time  would  amount  to  six  seconds  (about  two  one  millionths  of 
the  pa  t  1  period). 

To  date  this  author  has  not  been  able  to  find  any  past  or 
current  research  on  the  capacilities  of  HFDF  systems  to 
exploit  short  duration  signals.  This  is  probably  due  to 
several  reasons.  The  primary  reason  is  that  HFDF  engineers 
are  still  absorbed  in  the  more  basic  problem  of  improving 
HFDF  against  medium  and  long  duration  skywave  signals. 
Especially  in  the  case  of  narrow  aoerture  HFDF  antenna 
systems,  there  remains  considerable  need  for  improvement.  In 
the  case  of  wide  aperture  HFDF  antenna  systems,  tne  problem 
of  short  duration  signals  seems  to  be  tractable.  However, 
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there  does  not  appear  to  exist  a  comprehensive  study 
addressing  this  problem,  and  to  conduct  such  a  controlled 
study  of  skywave  propagated  signals  would  be  expensive. 
Notwithstanding  these  difficulties,  the  snort  duration  signal 
could  easily  become  an  acute  tactical  problem  for  the  side 
that  carnot  exolcit  it,  and  the  problem  tnerefore  deserves 
immediate  attention. 

3.  FURPCSZ  OF  THESIS 

The  general  question  of  interest  is  how  good  are  existing 
HFDF  systems  at  determining  lines  of  bearing  on  sno^t 
duration  signals.  Any  complete  examination  of  an  HFDF  system 
requires  one  to  investigate  the  characteristics  of  either 
wide  or  narrow  aperture  antennas.  Additionally,  one  must 
examine  the  problems  of  site  location,  signal  acquisition 
system,  "eceiver  demodulation,  tearing  sense  circuits  and  t.ie 
noise  environment.  Performance  must  also  be  determined  f r" 
ground  versus  skywave  and  multipath  versus  single  patn.  This 
would  be  an  enormous  task  i?  all  current  HFDF  systems  ve^e 
considered.  The  scope  of  this  effort  is  much  mo^e  restricted. 

The  purpose  of  this  thesis  is  to  investigate  the  mean  and 
variance  of  tearing  estimates  for  short  duration  skywave 
signals  received  with  a  narrow  aperture  HFDF  system.  Several 
statistical  procedures  are  developed  and  compared  with  the 
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intention  of  discriminating  between  reliable  and  unreliable 
data  and  calculating  the  best  bearing  estimate  from  tne 
reliable  data. 

The  problem  of  how  to  acquire  a  snort  duration  signal  or 
how  to  interface  such  an  acquisition  system  to  an  HJTDF  system 
will  not  be  addressed  in  this  report. 
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II.  TH70PETICAL  COJ.S  I DZ?  AT  IONS 

A.  HI  I  :i  JPFQU2MCY  SKY'.’/AVZ  CHANNEL 

High  frequency  rHF)  is  the  region  from  three  to  thiry 
megahertz.  Its  complexity  is  due  to  many  natural  pnencnera 
which  are  interdependent,  complex  in  themselves  .  sore  poorly 
understood,  cosmic  and  microcosm ic  in  extent  and  difficult  to 
neasu^e.  The  primary  complexity  fo"  the  HFD?  engineer  resides 
ir.  the  ionosphere.  The  ionosphere  can  oe  considered  an 
inhomogeneous  plasma  that  surrounds  a  sphere  of  finite,  out 
variable,  conductivity  and  separates  the  sphere  from  free 
space  and  numerous  solar  ana  galactic  sources  of 
disturbances.  The  ionosphere  has  been,  and  continues  to  be.  a 
subject  of  considerable  research.  References  1  througn  1  are 
a  rich  resource  of  information  on  the  ionosphere  and  research 
reports  are  added  monthly,  but  th^°  scops  of  the  problem  is 
imrense.  To  predict  accurately  ionospheric  conditions,  one 
must  be  equipped  with  mors  than  the  physics  of  the 
ionosphere.  The  physics  provides  the  equations  of  tne  system, 
but  the  forcing  functions  and  the  boundary  conditions  must  be 
sufficiently  measured  to  forecast  accurately. 

The  forcing  functions  are  tne  solar  flux,  gravitational 
waves,  weather  related  dynamics  and  the  signal  of  interest. 
The  ionosphere  is  usually  modeled  in  terms  of  electron 
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concentrations*  t.heT'efC"e ,  phenomena  which  affect  tne 
concentration  or  the  excitation  of  the  electrons  drive  tie 
system.  The  solar  flux  is  primarily  a  diurnal  phenomenon;  its 
impact  is  strongest  in  tne  portion  of  the  plasma  illuminated 
by  tne  sun.  This  flux  is  made  up  of  electromagnetic  energy 
and  streams  of  particles.  In  the  case  of  sun  spots  and  solar 
flares,  there  are  often  increased  emissions  t.iat  tend  tc 
disturb  the  normal  structure  of  the  ionosphere.  (Perhaps  it 
would  be  mo~e  accurate  to  state  that  tne  disturbance  is  tc 
our  model  of  the  i oncspher,e .  1  The  solar  disturbance  evolves 
in  three  stages.  The  first  is  the  impact  of  electromagnetic 
energy  in  tne  ultraviolet  and  x-ray  ranges  teat  causes  an 
increased  electron  concentration  in  the  lowest  electron  layer 
(D-layer'.  The  second  effect  is  the  arrival  of  high  energy 
protons  and  alpha  particles  that  also  increase  tne  D-layer. 
The  duration  of  the  disruptions  due  to  these  two  erases  is 
limited  tc  several  hours.  Tne  third  chase  is  tie  arrival  of 
low  energy  protons  and  electrons  which  shower  the  earth  in 
patterns  molded  by  the  Earth's  magnetic  field.  In  this  pnase. 
which  may  last  as  long  as  several  days,  the  ionosphere 
experiences  magnetic  storms,  an  increase  in  the  D-layer . 
sporadic  conditions  in  the  next  higher  F-layer  and  tne 
spectacular  aurorae. 

Acousto-gravity  waves  constitute  a  forcing  iunction  of  a 
different  scope.  Periodic  variations  in  tne  dynamics  of  the 
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Earth-mocn-sun  gravitational  system  and  isolated,  anomalous 
gravitational  activity  on  the  Earth  combined  with  HE  acoustic 
waves  (The  Mt.  St.  Helen  eruption  was  a  recent  source  of  such 
waves)  exert  forces  that  distort  the  general  concentric 
spherical  form  of  the  electron  plasma  layers.  The  distortions 
are  not  only  static.  Traveling  ionospheric  disturbances  a~e 
not  uncommon,  and  their  effect  is  to  create  a  doppler  shift 
on  transmitted  signals.  If  the  ionospheric  disturbance  is 
tilted,  the  ray  trace  of  a  transmitted  signal  will  be  bent  in 
azimuth. 

The  third  mentioned  forcing  function  is  the  we  a  tner .  T.ne 
dynamics  of  the  weather  affect  the  pressure,  the  temperature 
and  the  nixing  of  the  atmosphere.  These  three  factors  in  turn 
have  a  significant  impact,  on  the  election  concentrations, 
particularly  the  concent "at  ions  at  the  lower  altitudes.  The 
weather  is  also  a  very  important  factor  in  high  frequency 
ranges  because  it  is  a  noise  source.  Much  of  the  high 
frequency  background  noise  is  attributed  tc  thunderstorm 
activity  which  is  continually  occurring  at  seme  point  on  the 
Earth.  (It  should  b^  noted  that  m^st  of  the  electromagnetic 
energy  cf  a  thunderstorm  is  in  the  V L?  region.) 

Manmade  signals  are  one  of  the  smallest  forcing  functions 
acting  on  t  ne  ionosphere,  but  they  a  re  naturally  of  great 
interest.  The  target  signal  injects  itself  into  the 
ionosphere;  it  operates  on  the  ionosphere  and  is  operated  on 
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by  the  ionosphere.  The  study  of  this  interaction  has  lead  to 
a  description  of  the  change  of  transmitter  antenna 
polarization  to  elliptical  polarization,  the  phenomena  of 
refracted  high  frequency  waves,  multipatn  interference  and 
the  concepts  of  maximum  useable  frequency  (MUF),  lowest 
useable  freauency  (LHP)  and  optimum  working  frequency  (POT'. 

The  ionosphere  is  a  system  vitn  largely  fluctuating 
boundary  conditions.  The  surface  of  the  Partin  is  tne  only 
boundary  that  can  be  considered  fixed  with  respect  to  daily, 
seasonal  and  eleven  year  solar  cycles.  Other  boundary 
conditions  are  much  more  dynamic.  Of  these,  the  layering  of 
electron  concentrations  is  primary.  The  inner  two  layers,  D 
and  P,  which  are  mostly  the  result  of  solar  electromagnetic 
radiation  have  been  mentioned.  Tne  outer  layer.  P,  vaicn 
often  is  subdivided  into  an  PI  and  72  layer  is  relatively 
more  stable.  It  remains  when  the  portion  of  the  ionosphere  of 
interest  rotates  into  the  sola”  umbra  and  tne  v  and  7  layers 
disperse.  The  D  and  7  lasers  during  iayligrit  a”e  responsible 
for  the  non-ieviat i ve  attenuation  of  much  of  tne  HP  spectrum 
of  interest  (3-12  MHz'.  The  dispersion  of  the  D  and  S  layers 
permits  the  P  layer  to  become  a  virtual  reflector  situated  at 
altitudes  typically  from  220  to  422  km.  (Tne  actual  mechanism 
of  propagation  through  tne  7  layer  is  refraction  whicn  can  be 
modeled  as  reflection  from  a  virtual  height  greater  t nan  the 
actual  zenith  cf  the  tending  ray.''  ?  propagation  opens  up  tne 
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evening  airwavs  to  long  distance  communications  and  attendant 
long  distance  HFDF  in  the  3-12  MHz  range.  For  the  engineer 
this  is  a  mixed  blessing. 

The  engineer's  interest  in  scywave  propagation  is  in 
the  difference  between  the  direction  of  arrival  of  the  target 
signal  and  the  g~eat  circle  hearing  to  the  target  and  in  tne 
variance  of  the  measurements  of  the  angle  of  arrival.  Aside 
f"Oir  the  equipment  limitations  and  site  location  distortions 
and  reflections,  -any  of  the  errors  ard  variances  that  r,°ed 
to  he  "esolved  to  improve  DF  are  due  solely  to  the 
i  onosphere . 

In  the  evening,  targets  of  interest  in  the  3-12  MHz  hand 
can  he  exploited,  hut  there  is  a  considerably  greater  chance 
of  interference  from  other  discrete  sources  or  from  general 


noise 

s ou roes  . 

Additionally,  tnere 

is 

increased 

complex! ty 

when 

s ignal s 

routinely 

arrive  aft 

er  t 

wo  or  three 

hops  which 

correspond  to 

maximum 

distances 

of 

£033  and 

1 2  A  3  £  xm. 

respectively.  Over  these  distances  the  errors  and  variances 
due  to  interference ,  fading,  tiltirg  a^d  scattering  increase 
to  a  point  that  even  wide  aperture  antennas  cannot  produce 
useful  fix  information. 

An  important  consideration  for  narrow  aperture  antennas  is 
that  the  effective  aperture  of  one-quarter  wavelengtn  at  c-? 
MHz,  e  typical  lcnghaul  daytime  frequency,  oecomes  a 
cne-si xteenth  wavelength  at  5  MHz.  a  nighttime  frequency.  Tne 


loss  of  effective  aoerture  further  exaceroates  the  problem  cf 
deterning  a  bearing  and  its  variance.  The  effective  height  of 
the  arter.r.a  is  also  a  function  of  frequency;  therefore,  ore 
can  expect  the  array  pattern  to  change  with  tne  change  in 
operating  frequencies. 

3.  NARROW  A? ??.TrJR?  D?  AM  TUN NAS 

The  knowledge  of  the  ionospnere  has  grown  extensively  in 
the  past  forty  years.  Investigators  can  now  feel  reascnaoly 
comfortable  with  the  developed  models  and  the  improved 
sensors.  especially  the  ex t ra-ter~es  t riaily  sensing 
satellites.  General  predictions  are  possible  and  a  new 
favorite  computer  aid  is  the  software  that  predicts 
propagation  and  displays  ray  tracings  (see  ?ef.  13  and 
Appendix  C ' .  An  IF??  engineer  can  review  the  general 
propagation  scheme  with  an  assurance  that  he  umerstanis 
sufficiently  tne  problems  presented  by  the  ionosphere.  But  in 
the  case  of  narrow  aperture  HFDF  antennas  one  must  guard 
against  the  feeling  of  confidence  induced  by  a  knowledge  cf 
the  general  situation.  One  is  reminded  of  the  situation  where 
a  blind  man  feeling  the  t runs  of  an  elepna.it  attempts  a 
general  description  of  the  elephant.  Ir.  tne  case  of  a 
1. 5-meter  aperture  antenna  sampling  a  wavefront  in  tee 
60-meter  band,  the  dimensional  comparison  with  a  nani  and  an 
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elephant  is  accurate. 

There  are  two  commonly  used  tyoes  cf  narrow  aperture 

antennas.  One  type  relies  on  amplitude  comparison  to 
determine  direction  of  arrival  and  the  second  type  relies  on 
phase  corpar'ison.  An  example  of  the  former  is  tne  simple  loop 

and  of  the  latter  is  the  Adcocs  .  (leferer.ee  5  points  out  tnat 

the  phase  and  amplitude  distinction  is  not  clear  cut  in  the 

ase  of  the  simple  loop  is 

E  v 


V 

Figure  2 

Simule  loop  sensing  direction  of  arrival 

The  direction  of  arrival  of  the  signal  is  determined  tv  tne 
relative  orientation  of  the  loop  and  the  horizontal  component 
of  the  magnetic  field. 
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The  case  of  the  Adcock  (actual lv  one  half  of  a  U-Adcock) 


is 


illustrated  as: 


v  s  (zttcA I'k)  sm  e 
9  s\nK  (vk') 

K=  ^  / 2iTci 


Figure  3 


Simple  Adcock  sensing  direction  of  arrival 


The  direction  of  a^ival  of  the  signal  is  determined  'ey  tne 
phase  difference  between  the  two  elements. 

The  two  examples  above  only  serve  to  illustrate  how 
direction  of  arrival  information  is  determined.  Real  systems 
emplov  more  elements  to  resolve  ambiguities,  improve  accuracy 
and  enhance  resistance  to  noise.  The  point  is  that  tne 
fundamental  process  relies  on  an  element  sensing  amplitude  or 
phase.  This  fundamental  process  is  in.  tu~n  the  fundamental 
difficulty  for  narrow  aperture  HFD7  antennas. 

In  Ref.  5  Ge thing  uses  computer  simulation  to  plot  wave 
interference  of  multimode  signals  in  terms  of  surfaces  of 
constant  phase  (C?)  and  constant  amplitude  { CA).  Tor  the 
ideal  case  of  a  single  specular  component  with  no  scattering, 
the  surface  of  C?  and  CA  is  a  plane  whose  normal  is  the 
direction  of  oropagation.  In  the  case  of  two  rays.  tne 


interference  patterns  repress  n  ted  by  the  surfaces  of  0?  ana 
CA  vary  with  the  anrular  separation  of  the  rays  in  elevation 
and  azimuth.  In  all  of  the  patterns  presented  in  Fef.  5  in 
which  the  amplitudes  of  the  two  component  rays  differ  by  only 
ten  percent,  major  distortions  to  the  ideal  planes  of  CA  arc 
C?  occur.  Approximately  planar  portions  of  the  surfaces  of  CA 
and  C?  extending  to  several  wavelengths  in  length  are  up  tc 
sixty  degrees  different  from  the  true  angle  of  arrival.  Tne"e 
are  also  '■cinks  in  the  ohase  fronts  that  vary  tne  phase  up  tc 
ninety  degrees  in  less  than  the  space  of  one  wavelength.  In 
the  cases  where  mo~e  than  two  rays  are  present,  the 
interference  patterns  become  much  mo^e  complex. 

It  is  obvious  that  the  spatial  extent  of  wide  aperture 
antennas  is  needed  to  resolve  suc.n  interference  patterns  in  a 
short  period  of  time.  Falser  and  Smith  in  Ref.  7  explained 
that  when  the  outputs  of  two  antennas  were  correlated,  the 
antennas  had  to  be  spaced  forty  wavelengths  in  tne  case  of 
single  hop  and  ten  wavelengths  in  the  case  cf  multi  nop  to 
lower  the  correlation  coefficient  to  C.5  .  For  a  narrow 
aperture  antenna  to  detect  phase  or  amplitude  distortions  of 
this  magnitude,  the  time  of  observation  must  be  relatively 
long.  However,  it  is  necessary  to  detect  such  distortions  to 
permit  an  assessment  of  reliabilty  to  be  assigned  to  bearings 
measured  in  distorted  fields. 

The  description  above  of  interference  patterns  was  for  two 
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sources  nearly  equal  in  amplitude.  The  condition  of 
comparable  amplitudes  is  one  that  results  in  severe 
distortion.  As  the  amplitude  ox"  one  of  tne  rays  becomes 
substantially  less  than  the  othe7'  ray .  the  interference 
pattern  approaches  the  ideal,  undistorted  planar  pattern 
(implicit  is  tne  assumption  that  single  mode  scatter  is  also 
very  weak).  Assuming  that  at  least  one  of  the  rays  of  c  two 
ray  interference  pattern  is  fading,  it  can  be  expected  that 
for  short  periods  of  time  ideal  planar  CA  and  C?  wavefronts 
can  be  observed.  There  is  no  ready  means  of  identifying  these 
moments;  however,  if  the  fading  is  random,  the  planar  CA  and 
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rotation  induced  elliptical  polarization  the  fade  rates  are 
measured  in  seconds  per  cycle.  In  Ref.  12  polarization  fading 
with  periods  of  1?  seconds  and  23  db  fade  depths  were 
reported  as  common.  If  two  or  more  rays  are  present,  fading 
is  measured  in  cycles  per  second.  This  indicates  that  a 
narrow  aperture  HFDF  system  will  require  approximately  a 
second  to  recognize  the  fading  condition  if  strong  multipath 
interference  exists.  The  time  required  to  average  the 
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interference  pattern  is  related  to  the  polarization  fading  of 
t.ne  dominate  mode.  The  amplitude  of  a  polarization  faded 
signal  is  a  stocnastic  process;  therefore,  tnere  is  no 
deterministic  functional  relationship  between  time  and 
fading.  A  measure  of  the  rapidity  with  wnich  facing  is 
fluctuating  can  be  expressed  in  terms  of  a  fading  power 
spectrum.  (Section  5.4.3  of  Ref.  1  discusses  tne  concept  of 
fading  power  spectrum.)  If  tnere  is  a  large  portion  of  tne 
"fading  power"  in  the  higher  frequencies  (130  to  1 i d k  Hz  , 
the  facing  is  fast.  If  the  "fading  power”  is  primarily  in  tne 
the  ?.l  to  1.0  Hz  region,  the  fading  is  slow.  In  tne  case  of 
polarization  fading  it  nas  already  been  noted  t.nat  fading  is 
typically  in  the  seconds  per  cycle  range.  Therefore,  an 
antenna  which  does  not  have  sufficient  spatial  aperture  to 
average  interference  patterns  must  rely  on  fading  to  permit 
time  averaging.  The  time  required  for  averaging  is  a  function 
of  the  interference  pattern  and  appears  tc  be  on  the  oraer  of 
five  to  ten  seconds. 

A  measurement  experiment  reported  by  Bain  in  Ref.  3 
demonstrated  how  time  averaging  of  bearings  reduced  the 
variance  associated  with  the  mean  bearing.  Using  a  U-Adccci 
with  buried  feeders,  bearings  or.  skywave  signals  were 
recorded  at  five  bearings  per  second.  An  autocorrelation  <'>f 
the  bearings  was  computed  and  the  resulting  curve  was 
approximated  by  the  exponential  expression: 


R(  T  )  =  exp(  -  T/t.  ) 


wnere  is  a  parameter  associated  with  fitting  an 
exponential  curve  to  the  measured  hearings.  The  formula 

relating  the  variance  of  the  mean  hearing  (  <5 1  and  tne 

x 

variance  of  a  single  observation  ((5)  is: 


where  T  is  the  time  i nt e^va 1  over  which  the  hearings  were 
averaged.  Bain  reported  that  for  (T^ =  0.55  (corresponding  to 
considerable  hearing  fluctuation),  the  variance  was  reduced 
by  a  factor  of  10  in  12  seconds.  The  12  second  duration 
roughly  co",'esponds  in  order  of  magnitude  to  the  reciprocal 
of  an  average  fade  ’'ate. 

G.  SUMMARY  0?  THEORETICAL  CONSIDERATIONS 

The  narrow  aperture  HE'D?  antenna  is  physically  limited  to 
tire  averaging  operation  against  skywave  signals.  In  tne  case 
of  ideal  ionospheric  propagation,  the  antenna  system  can 
perform  within  equipment  and  site  limitations.  If  tne  site 
errors  are  known,  the  equipment  and  array  calibrated  and 
tnere  is  a  good  SNR,  average  bearing  errors  of  0.5  to  1 
degree  and  variance  of  5  degrees  squared  should  be  uossible. 
If  multipath  oropagation  exists  with  fading  up  to  20  db.  time 
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averaging  over  at  least  ten  seconds  with  sampling  at  about 
five  per  second  should  reduce  most  of  the  variance  due  to  tae 
complex  inter:’ e^ence  patterns. 

The  difficulty  of  obtaining  accurate  HfD?  against  short 
duration  signals  using  a  narrow  aperture  array  is 
considerable.  In  the  case  of  multipath  interference  in  which 
at  least  two  rays  are  comparable  in  amplitude,  the  D F  error 
on  a  short  duration  signal  with  only  one  sample  bearing  could 
be  up  to  ninety  degrees.  This  is  tne  extreme  of  tearing  error 
due  to  phase  and  amplitude  front  distortion.  In  a  more 
hospitable  multipath  environment  the  system  performance 
should  be  much  better,  but  there  is  little  experimental 
evidence  by  which  one  can  assign  typical  tearing  errors  ard 
variances.  The  analvsis  of  the  narrow  aperture  antenna  system 
in  the  follovirg  sections  provides  performance  data  on  an 
experimental,  state  of  the  art  system. 
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The  Electr omagnet ics  Division  of  Southwest  research 
Institute  :'5'vRI),  located  in  San  Antonio.  Texas.  has 
developed  and  tested  a  new  design  for  a  narrow  aperture  RED? 
antenna  system  to  operate  against  hoth  ground  wave  anc. 
shyvave  signals.  The  significance  of  tnis  nev;  cesign  is  tnat 
it  is  a  nast  noun  table  narrow  aperture  ar.  tenna  that  is  a 
fixed  array.  There  are  no  moving  parts;  therefore,  it  is 
ideal  for  the  shipboard  environment.  The  elements  of  tae 
array  are  simple  loops  and  soaced  loops.  The  latter  will  be 
shown  to  have  polarization  independent  qualities  and. 
therefore,  to  be  ideal  for  exploiting  s’sywaves .  The  array  and 
the  associated  instrumentation  of  the  system  are  high  speed 
and  computer  controllable. 

The  primary  ^efe-’ence  for  the  analysis  that  follows  is  an 
in-house  report  prepared  by  the  system  architects  [Ref.  D <  . 
The  system  herein  described  has  beer  patented.  The  author  cf 
this  thesis  has  visited  the  San  Antonio  site  and  has  operated 
the  HFDE  system  with  the  assistance  of  the  3WRI  personnel. 
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3.  THEORY 

To  understand  the  operation  of  the  spaced  loop  array  one 
roust  review  the  theory  of  the  simple  loon.  The  i  igure  or.  the 
next  page  depicts  a  simple  loom  set  in  a  coordinate  system 
with  an  incoming  signal  ray  (3  field  components  lace  led  2v 
and  Eh).  The  angle  phi  ( 4>)  is  the  azimuth  measured  in  tae  XY 
plane  of  the  incoming  ray.  The  plane  of  the  loop  is  aligned 
with  the  XZ  plane.  The  incidence  angle  tneta  O'1  is  measured 
in  the  plane  defined  ’ey  the  Z  axis  and  the  signal  ray.  The 
signal  is  considered  to  have  a  vertical  and  horizontal 
electric  field  component  (2v  and  Eh).  The  expression  for  the 
output  voltage  (sinusoidal  input,  time  variation  suppressed; 
is  : 

VI  =  -  (  T  v  ' )  C  C  5  4>  -  f  Eh')  (  cose  3  IN  $  )exTH  j  <J>  h  (i) 
where: 

Zv ' :  relative  amplitude  of  the  vertical  component 
Fh ' :  relative  amplitude  of  the  horizontal  component 
4>  :  a  zimuth 

0  :  angle  of  inciderce 

^  h :  phase  of  horizontal  component  relative  to  t ne 
vertical  component 
VI:  simple  loop  output  voltage 

With  respect  to  sky-wave  signals,  tne  significance  of  equation 
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Simple  Loop  and  Coordinate  System 


(1)  is  that  it  is  polarization  dependent.  Hi'PR  systems 
(generally  ~ely  on  isolating  a  null  in  the  array  pattern  that 
can  be  related  to  the  azimuth  of  the  incoming  signal.  Tne 
null  used  should  or.l1'’  be  a  function  of  the  target's  bearing. 
The  simple  loop  works  well  with  ground  waves  fo”  which  case 
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VI  =  -Ev'  COS  <£> 

The  simple  loop  does  not  function  acceptably  against 
skyvaves.  In  the  case  of  skyvaves  the  loop  voltage  is  a 
function  of  the  two  spatial  variables,  theta  and  phi.  and  tne 
relative  phase.  The  nulls  c”°ated  by  these  three  variaoles 
are  too  numerous  and  the  available  measurements  too  sparse  to 
resolve  all  the  ambiguities. 

A  solution  to  tne  polarization  dependence  limitation  of 
the  simple  loop  is  to  combine  two  simple  loops  into  a  two 
element  interferometer  as  illustrated  in  figure  2.  The  loops 
are  corrected  in  parallel  with  con  os ing  phase.  Tne  output 
voltage  of  this  array,  Known  as  coaxial  spaced  loops,  can  be 
determined  by  nattern  mul t ini ic a t ion  [Ref.  13  ].  The  pattern 
of  the  spaced  loop  array  is  equal  to  the  product  of  tne  group 
pattern  and  the  element  pattern.  The  g~oup  pattern  cf  the 
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Figure  5 

Coaxial  Spaced  Loops  and  Coordinate  System 
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array  is: 


G-r  =  j  fi  a  S  IN  (  9  ) 3 1 N  (  0>  ) 
d:  separation  cf  the  tvc  loops 

/3  :  2  tt.'Ts 

7m  :  signal  wavelength 

New  let: 

Ev  =  jEv'  d/2 
Eh  =  jEn'  d/2 

This  permits  the  spaced  loop  output  voltage  to  he  written  as: 
V  a  =  -Ev  SIM  SSI  M  2$J-  Eh  (SIM  2©S  I M  4*)  exp  ■  j$h  )  (2) 


The  significance  of  equation  ( 2  ’>  is  the  existence  of 
polarization  independent  ri:ll5 .  The  output  voltage  equals 
zero  whenever  the  azimuth  ar.ule  equals  2  or  ISC  degrees.  The 
incidence  angle,  the  ~e].  at,  i.  ve  phase  and  the  relative 
amplitudes  of  the  electric  field  components  do  net  affect 
these  nulls.  It  is  due  to  the  interferometer  structure  teat 
these  nulls  exist;  they  are  therefore  called  interferometer 
nulls  to  distinguish  taem  from  the  simple  loop  nulls,  figure 
3 .  taken  from  reference  11,  graphically  displays  the 
polarization  independent  nulls  different  conditions  of 


incidence  and  polarization. 

Equation  (2N  is  an  expression  for  a  fixed 
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Figure  6 

Coaxial  Spaced  Loop  Patterns  as  a  Function  of 
Signal  Polarization  and  Angle  of  Incidence 
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tne  spaced  loop  array  in  the  coordinate  system.  To  make  the 


orientation  arbitrary  the  variable  alpha  is  introduced  into 

(2) . 

VaCX  =  -S’'  SIN  S I  M  2  (  <$)  -  <X)  +  Eh  S  I  N2  6  S  I N’2  (  $  -  *  (  3  - 

.'<$>-  DO:  relative  azimuth  angle 


The 

alpha 

variable 

perm  its 

the  expr 

°  s  s  i  o  r. 

o  f  t  he 

0  u  t  pi.  t 

volt 

age  of 

a  spaced 

loop  array 

oriented 

a  1  pha 

d egrees 

from 

the 

x-axi s 

to  be 

written  as 

s  h  o  v  r.  i  r 

(3)  . 

This  v  ill 

later 

allow  equation  (3)  to  express  the  output  voltages  of  more 
than  one  pair  of  spaced  loops  set  at  different  angles  in  the 
coordinate  system. 

3y  defining: 

C  =  E v  SIN© 

A3  =  Eh  S I N ( 2  ©  / ex p ( j  $  h '  V  2 
A 2  =  -C  SIN  2<$>  -  A3  CCS  2 

32  =  -A3  SIN  2$  x  C  COS  2  <*> 

Equation  (3)  can  be  written  as: 

Vaoc  =  A0  +  A 2  COS  2 <*.  +  P2  SIN  2^.  ( 4' 

This  form  permits  a  Fourier  series  interpretation  of  the 
spacea  loop  output  voltage: 

A 25  =  dc  term  of  the  output  voltage 
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A2  and  32  are  coefficients  of  tne  second  harmonic 

The  significance  of  (4)  is  that  for  a  pixel  value  of 
target  azimuth  and  elevation,  the  spaced  loop  vclta-e  as  a 
function  of  relative  azimuth  is  limited  to  a  second  carbonic 
of  the  relative  aziruth.  The  application  of  the  Nyquist 
sampling  criterion  ~eveals  that  the  voltage  pattern  can  be 
duplicated  by  four  sample  values.  Therefore.  a  spinnir.,; 
spaced  loop  can  be  synthesized  hy  a  rinimum  of  four  samples 
taker,  equally  spaced  through  35?  degrees  of  azirutn. 

The  solution  for  the  tearing  (the  azimuth  angle  pnij  is 
derived  from  equation  (4)  and  the  definitions  given  atove  for 


C  .  A 

rsj 

(i-* 

i\j 

and  32. 

Py  algebraic 

:  man i pula ti 

on  it 

i  s 

ie 

tt 

:rm  i  n 

that . 

Vi 

C  - 

*/-  ( 

A  2  *  B  2 

-  b  ?> 

)  (5 

) 

A 2 ,  A2  and 

32  will 

be  s 

h  c  w  r. 

to 

be 

measurea 

c  le  c 

ua  n  t  i 

t  i 

.  es  . 

is 

determ 

i  ned 

f  r  on 

ear 

a  t: 

.  on 

( 5  'i 

above . 

Us 

ir  g  t 

relationships  , 

SIN 
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AC 

)] 

f 

( C  )  ( A2  )  + 
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?2) 

n 

) 

(  5 

CCS 
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-[  l/( 

c  + 

A3 

)] 

[ 
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-  (CM 

32) 

] 

‘7 
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-  f 
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( e .  5  • 

AP.CTAN 

'  SIN' 

2  (J> 

!  r 

:os 

2$  ) 

n  •  13  2 

n  = 

2, 

,  1 

Inherent  in  equation  (4)  and  made  obvious  in  equation 
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are  four  null  ambiguities.  There  a^e  two  nulls  162  decrees 
apart  that  can  be  attributed  to  the  simple  loops  and  two 
nulls  163  decrees  apart  that  are  the  interferometer  nulls. 
The  SWF  I  analysis  shows  that  by  aiding  tne  simple  loop 
phasors  into  the  analysis.  tne  simple  loop  nulls  can  be 
determined  and  then  discarded.  By  comparing  the  spaced  loop 
output  to  the  simple  loop  output,  the  c^r^ect  interferometer 
null  which  represents  the  desired  bea^ir*-  can  be  identified 
The  engineers  at  S>,DI  used  the  ideas  they  developed  above 
to  design  a  fixed  spaced  loop  array.  Tne  undesirable 
mechanical  feature  of  the  rotating  spaced  loop  was  eliminated 
by  using  four  spaced  loops  fixed  in  an  array  to  synthesize 
rotation  as  shown  below. 

<x 


2  45 

9? 

135 

*  o 

© 

130  225 

2?  ? 

315 

Figure  7 

Spaced  loot)  array  geometry 
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The  Nyquist  criterion  requires  a  minimum  of  four  samples  to 
syntnesize  equation  (4).  This  could  reliably  be  accomplished 
by  three  nairs  of  spaced  loops,  but  to  provide  for  additional 
reliability  in  the  presence  of  noise,  a  four  pair  spaced  loop 
array  was  constructed. 

Assuming  the  orientation  given  in  the  diagram  above,  one 
can  determine  A '2,  A2  and  32  in  terms  of  tne  individual  spaced 
locus.  Solving 


Vat*  = 

AT 

+  A 2  CCS 

2  <X  + 

32  S 

t  erns 

of 

aluha  yields. 

= 

2 

Va2  = 

A  2  + 

A  2 

*  = 

45 

Va45 

=  A  2 

+  32 

OC  = 

52 

V  a  92 

=  A0 

-  A  2 

CA  = 

135 

V  a  135 

=  A  2 

-  32 

wnere  VaT.  7 a 45,  Va90  and  Val35  are  the  phasors  of  tne  spaced 
looms  in  the  array  shewn  above. 

This  provides  four  equations  to  solve  for  three  unknowns.  One 


solves 

for  A  0 . 

A 2  and 

32  by  the  following  equations. 

A0 

=  (2.25) 

( V  a  0 

'/  a  4  5  V a 9 2  x  7 a  1 3 5  : 

A  2 

II 

CJl 

( V  a  '2  - 

7a92 ) 

32 

=  (2.5) 

( Va45 

-  7  a  1 3  5  ) 

After  A 7 ,  A 2  and  32  a~e  determined  f^om  the  phasor  equations 
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above,  they  are  substituted  into  equations  (S'  and  (7,  vhi oh 
in  turn  are  used  to  solve  equation  (S)  for  the  tour  possible 
bearings . 

Rurtner  algebraic  and  trigonometric  analysis  detailed  in 
Reference  shows  that  the  simple  loco  nulls  can  be 

determined  by. 

(X-  ARC  TAM  (  VL0/-VL90  )  =  <$  -  A  1C  TAN  (  C  '  *Z  ) 

wnere  VL2  and  VLS0  are  the  ohascrs  of  the  two  simple  loop 
pairs  . 


Once  the  simple  loop  nulls  are  Known.  tne  sign  of  C  in¬ 
equation  (5)  can  be  determined.  This  in  turn  leads  to  the 
unambiguous  selection  of  the  proper  interferometer  null. 

4>  =  AFCTAM  (  VI /V2  ) 


where , 


vi  =  [j/(-a/  -  cz)1 

V 2  =  f  .1  / ( -A 01  -  CX)] 

it  was  noted  above 
solve  for  three 


f  (  C  )  (  VL9<3 )  -  ( Ad  )  t  YLd  )  ] 

(  (  C  )  ( V  L0 )  *  ( A3  )  ( VL90 )  l 

that  four  equations  are  available  to 
unknown  coefficients.  The  additional 


information  permits  two  separate  solutions  fo~  the  AO  term. 


f  A01  =  ( ? .  5  '  (  Va?  -  Vase  '• 


r A0l  '  =  (7.5-  (  Va45  +  Val25  1 

The  difference  between  these  two  A0  terms  snould  ideally  oe 
zero.  I?  the  iiffe"ence  is  rot  zero  thene  is  an  inconsistency 
within  the  system.  This  difference  is  called  tne  A4  term 
because  it  corresponds  exactly  with  the  coefficient  of  tne 
fourth  harmonic  of  a  Fourier  series  expansion  of  tne  spaced 
loop  pattern  in  azimuth.  The  A4  term  is  therefore  an 
important  parameter  in  determining  bearing  quality. 

C.  SYSTEM  DESI1N  AMD  INSTRUMENTATION 

The  array  of  spaced  loops  and  simple  loops  suitaole  for 
mast  mounting  shown  in  figur°  8  was  built  by  S.vRI.  There  are 
four  pai^s  of  spaced  loops  in  the  lower  bay.  Each  pair 
consists  of  40  inch  high  tv  2?  inch  wide  simple  loops 
separated  by  50  inches.  The  output  of  these  diametrically 
opposite  loops  are  connected  in  parallel  opposition.  The 
simple  loops  in  the  upper  bay  are  used  to  resolve  ambiguities 
in  the  tearing  algorithm.  These  diametrically  opposite  simple 
loops  are  connected  in  parallel  assistance.  The  reference 
antenna  is  synthesized  by  quadrature  addition  of  the  simple 
loops . 

A  block  diagram  of  the  equipment  suite  is  drawn  in  figure 
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Coaxial  Spaced  Loop  Mast-Mountable  Array 


n  n 


QH: 

Quad¬ 

rature 

Hybrid 


PD: 


Power  Divider 


Figure  9 

Block  Diagram  of  SWRI  Spaced  Loop  Antenna 
System  Instrumentation 
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S.  Tne  RJ  sequencer  is  a  computer  controlled  switch.  that  is 
necessary  to  provide  nigh  speed  switching  of  the  different 
elements  of  the  array.  A  dual  channel  receiver  is  used  to 
provide  a  receiver  channel  for  the  reference  signal  and  a 
channel  for  tne  loop  voltages.  Tne  predetected  output  of  botn 
channels  of  the  receiver  is  monitored  by  a  pnasemeter  that 
provides  a  digital  measurement  of  tne  onase  of  the  antenna 
elements  with  ~es  pee  t  to  the  ~ef e^ence .  The  detectors  are  a 
pair  of  precision  peak  detectors.  The  detector  output  is 
sampled  and  digitiz°d  by  the  analog  to  digital  converter.  The 
digital  data  is  "oute‘J  to  a  minicomputer  for  processing,  iacn 
data  frame  is  approximately  20  ms  ir.  duration.  Tne  data  frame 
consists  of  the  six  comple-r  numbers  representing  tne  six 
voltage  phase's  ' V a ?  .  Va 45,  Ya90.  Val35.  VL0  and  VL90  ■  .  Net 
all  of  tne  data  frames  are  acceptable.  Tne  voltage  values 
must  te  within  the  linea~  ~ange  of  the  detectc1”  and  receiver 
circuits.  The  acceptable  data  frames  are  tne  input  to  tne 
algorithms  that  determine  the  hearings  and  resolve  the 
ambiguities . 
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IV.  DATA 


A.  DATA  PILES 


The  SWPI  equipment  suite  is  a^ranmed  so  that  data  measured 
from  the  spaced  loop  antenna  instrumentation  is  stored  on 
magnetic  disk.  This  permits  the  DP  operator  to  postprocess 
the  data  using  statistical  techniques  to  derive  a  more 
accurate  tearing.  Also  available  is  the  capability  of  mass 
storage  on  magnetic  tape.  It  was  on  magnetic  tape  that  3  V.  R I 
provided  the  Naval  Postgraduate  Sc.ocol  with  nine  files  of 
data  in  1S7S  and  four  files  in  I960.  The  1979  data  consists 
of  the  following  files.  The  test  source  was  a  transmitter 
placed  close  to  the  array  to  provide  a  ground  wave  ir. 
approximately  the  same  direction  as  'V'.vV. 


File  Number 

Frequency 
(MHz  ) 

S  ource 

T  i  n  6 
( CST  ) 

Date 

i 

10 

fc'WV 

3E  : 33 

2/13/79 

2 

23 

vwv 

12:15 

2/13/79 

3 

15 

j/'.v  V 

12:33 

2/13/79 

4 

5 

VW  V 

37:33 

2/14/79 

5 

15 

wv 

39:63 

2/14/  79 

f 

5.71 

Test 

7 

13.31 

Test 

8 

15.31 

Test 

9 

23.01 

Test 
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1983 

files 

are  as  fol 

lows  : 

Rile 

Number 

F  requenc y 

Source 

Time 

Date 

(MHz  ) 

(  CST  ) 

10 

5 

!i  ,v  V 

22*40 

2/6/80 

11 

10 

ww  v 

07  : 40 

2  /5  /3  0 

12 

15 

VW  V 

09:00 

2/5/80 

13 

8.566 

KLC 

08  : 15 

2 -7 /c  2 

It  is  known  that  file  1?  is  a  low  SNR  data  set.  Files  11  and 

12  a  "e  data  sets  with  SNR's  in  excess  of  22  dS. 

B.  a  tfV  AND  KLC 

V*  V*  V  is  an  ideal  target  because  it  is  an  amplitude 

modulated  signal  with  no  carrier  suppression.  For  the 
majority  of  the  hourly  information  duty  cycle,  trie 
information  modulated  on  the  carrier  is  simple  4  40,  5 2  2  and 

600  Hz  tones’*  ticks;  and  occassional  voice  announcements.  The 
ThV  signal  is  stable  to  +/-  two  marts  in  t 1<  ,  .  and  it  is 

available  on  5,  10,  15  and  20  MHz  for  24  hours  a  day.  Tne  vva// 
signal  is  transmitted  from  Boulder,  Colorado,  which  is 

geographically  fixed  at  43.8  M  and  105.1  *.  The  true  bearing 
of  the  greet  circle  arc  passing  througn  San  Antonio.  Texas, 
and  Boulder  is  336.7  degrees,  and  the  length  of  tie  arc  is 
1367  km. 

KLC  is  a  manual  morse  ship-t o-snore  station  transmitting 
from  a  platform  in  the  Gulf  of  Mexico.  This  signal  was  cnosen 


as  a  target  because  of  the  on-off  tewing  (OOK)  modulation  and 


because  its  relatively  short  distance  from  San  Antoni 


results  in  a  shywave  with  a  high  angle  of  arrival  at  tne  S'a'PI 
antenna  array.  The  COS  modulation  is  important  because  it  is 
a  favorite  mode  of  tactical  communications;  it  is  brief  and 
reliable.  The  true  bearing  f^om  San  Antonio  to  EEC  is  3S9  S 
degrees  and  the  distance  is  372  -cm. 

C.  DATA  RECORDS 

Each  file  consists  of  13,220  records.  Each  record  consists 
of  the  A3,  phase.  A4  and  tearing  terms  calculated  from  frames 
of  data  (six  voltage  phasors)  wnich  were  generated  every  33 
ms.  The  A3  in  the  data  is  a  normalized  version  of  tne  A3 
explained  in  the  previous  section.  In  that  section, 

A3  =  (2.5'  3n(5IN2 ©  )exp( j $ n) 

The  data  on  the  tape  is  A3  normalized  by  the  factor  ZvSI'i  , 
which  yields: 

A0n  =  (  Eh/Ev  ]  COS  0  exp(j4>h' 

This  is  a  complex  number  of  wnicn  only  the  magnitude  is  use’. 
All  further  references  to  AO  in  the  data  and  analysis  section 
will  be  the  magnitude  of  the  normalized  value: 

A0  =  f  Eh/Ev  ]  CCS  0 
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It  should  he  noted  hens  that  the  A3  term  is  a  measure  of  tne 
amount  c?  horizontal  polarization  present.  If  Eh  is  greater 
than  Ev.  the  ~atio  Vh/Ev  will  tend  tc  ra-ce  A0  a  number 
greater  than  one.  If  the  ar_gl°  f"om  the  perpendicular,  tneta. 
is  large,  there  is  a  greater  effective  a~ray  aperture  in  the 


plane  of  Ev,  and 

k?  is  small0". 

I  f 

Ad 

is 

large . 

t,  ne 

horizontal  com pone 

n  t  is  dominant  . 

If 

A  T 

i  s 

sma 1 1 . 

t  ne 

vertical  component 

The  second  term 

dominates. 

of  the  record  is 

the 

phase . 

It  is 

t  ne 

phase  of  the  horizontal  electric  field  relative  to  the 
vertical  component  of  tne  field.  It  is  a  calculated  value 
betweer  -lc0  and  +  130  degrees.  If  this  phase  angle  is  a 
constant  zerc.  the  polarization  is  linear.  If  it  is  a  nonzero 
constant,  the  polarization  is  elliptical.  The  phase  value  is 
typically  noted  to  vary  randomly  within  a  limited  range  over 
short  tine  durations.  Over  durations  of  several  minutes.  it 
will  varv  over  the  entire  -130  to  ->-130  degree  range  due  to 
changes  in  Eh  and  Ev  oath  lengths  and  multimode  interference. 

The  third  term  ir.  a  record  is  the  amplitude  of  the  A4  terr. 
which  was  discussed  in  the  previous  section.  It  is  a  measure 
of  the  inconsistency  within  the  spaced  loop  HEDF  system.  It 
is  in  large  part  due  to  noise,  hut  it  can  also  be  tc  a 
limited  degree  a  measure  of  circuit  imbalance,  measurement 
e^ro ~ ,  componen t  failure ,  software  failure  and  site  error. 
Its  value  is  that  i+  is  a  measure  of  performance?  however,  it 
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is  not  a  system  diagnostic  tool. 

The  fourth  member  cf  the  record  is  the  calculated  bearing. 
It  is  an  integer  value  from  2  to  359  degrees.  Tor  the  .v.vV 
signals  this  beaming  is  the  system's  estimate  cf  the  angle  cf 
arrival  of  the  signal  wavefront  which  should  not  vary  far 
from  the  value  335.7  degrees.  7or  5LC  the  true  bearing  is 
239.9  degrees. 

D.  IONOSPHERIC  DATA 

No  ionospheric  sounding  information  was  available  for  tne 
time  periods  during  which  tne  data  was  recorded.  However, 
pronag a  t i on  information  was  nrovided  by  tne  Naval  Ocean 
System  Center  (NOSC'1.  Using  the  grown  sun  spot  number  for  tne 
data  recording  dates,  they  employed  a  computer  propagation 
prediction  program  known  as  PROPHET  to  orovide  ray  trace 
diagrams;  MU?.  LU?  and  power  predictions;  and  24  hour  line  of 
bearing  variance  curves.  This  data  represents  a  good  estimate 
of  propagation  conditor.s  between  Boulder  and  San  Antonio  for 
the  times  of  interest.  Examples  of  the  program  output 
provided  bv  NOSC  are  reviewed  in  Appendix  C.  Using  the  NC5C 
data,  propagation  information  for  the  *,  vV  files  is  taoulated 
in  TABLE  I . 
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*r 
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5 

10 

11 

12 


TABLE  I 

Summary  of  Ionospheric  Data 

Freq  Time  Date  Variance  !iUF  LUF  Ionospheric 
MHz  C-MT  degrees  MHz  MHz  Mode 

squared  Prediction 


10 

14-30 

2/13/79 

1 

17 

2 

.5  probable  multimode 

1 ,  2 ,  83  nops 

20 

IS:  15 

2/13/7S 

1 

22 

5 

possible  multimode 

1  hop 

15 

18  :30 

2/13/79 

1 

22 

K 

highly  probable 

multimode  182  a ops 

5 

13-00 

2/14/79 

3 

12 

2 

possible  multimode 

terminator  1.2.3  hops 

15 

17  : 3? 

2/14/79 

1 

22 

5 

probable  single 

rode 

2 :  40 

2/*-*/?0 

•7  m  4 

12 

2 

possible  multimode 

1  . 283  hops 

10 

13:40 

2/5/80 

3 

17 

2 

probable  multimode 

teT-rinato1-  1,2  hops 

1  5 

17:00 

2/5/90 

1 

23 

5 

single  mode 
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V.  A  MALTS  IS  OF  SWF  I  DATA 


A.  INTRODUCTION 

The  purpose  of  this  analysis  is  to  e?ar  the  capability  of 
the  SWP.I  space!  loop  RFj?  antenna  sister  to  determine  the 
angle  of  arrival  of  a  short  duration  signal.  A  short  duration 
signal  is  considered  to  be  from  l'?d  to  1000  ms  in  duration. 
It  is  important  to  note  that  this  is  not  a  gene ral  analysis 
of  the  performance  of  the  antenna  system,  the  bearing  a.id 
sense  algorithms  or  the  post-processing  algorithms  developed 
by  SWF  I .  It  is  also  important  to  recognize  that  system 
development  is  not  complete,  but  is  the  subject  of  on-going 
research.  The  data  provided  to  this  investigator  was  provided 
from  a  system  configuration  not  optimized  for  short  duration 
signals  or  for  some  of  the  target  frequencies  reco”iec  . 
Because  the  subject  of  this  analysis  is  short  duration 
skywave  signals.  the  following  analysis  takes  into 
consideration  the  need  to  make  maximum  use  of  trie  available 
data.  Whereas  SV/RI  algorithms  stress  a  bearing  selection 
process  that  eliminates  a  large  percentage  of  trie  data 
records  tc  enhance  the  reliability  cf  the  estimated  bearing, 
this  analysis  recognizes  that  a  signal  of  £00  ms  duration  is 
represented  by  cr.lv  1?  data  frames  and  that  some  compromise 
to  T’eliabi  li  ty  must  be  made.  The  term  reliability  in  t.nis 


report  is  used  as  a  measure  of  confidence  in  the  validity 
the  data.  If  one  associates  a  standard  deviation  of  53 
decrees  with  a  data  reco~d  and  2?  decrees  vita  a  second  data 
record,  the  latter  would  he  considered  mere  reliable. 

The  analysis  fi^st  concentrated  cn  examining  the  data 
primarily  tv  filtering  or.  the  A4  te^m  .  the  indicator  r f 
system  inconsis tencv.  If  the  Ai  term  is  small,  the  bearing  in 
that  record  should  be  considered  more  reliable  tnan  a  hearing 
associated  with  a  high  A 4  value.  Using  this  approach  a 
.FORTRAN  program  called  DR  ERR  (DR  EFRor)  was  developed  to 
examine  each  file  ana  report  the  average  bearing  error, 
standard  deviation  of  bearing  error  and  two  other  statistics 
of  short  duration  signals  with  respect  to  the  Ai  term.  A  full 
discussion  of  the  program  is  detailed  in  the  next  s act  ion.  It 
was  discovered  that  tae  A4  term  is  a  useful  parameter  for 
determining  bearing  reliability  in  tae  majority  of  cases. 
However,  when  the  A 4  threshold  is  set  to  only  allow  t.ne  data 
record  with  A4  approaching  close  tc  zero  (the  theoretical 
ideal).  there  is  not  a  large  probability  of  determining  a 
bearing  on  a  short  duration  signal.  In  an  attempt  to  improve 
on  the  sole  use  of  the  A4  term  as  a  r el i ab i 1 ty  indicator,  a 
probabilistic  likelihood  ratio  matrix  based  or.  all  of  the 
available  signal  parameters  was  employed  in  the  analysis.  The 
details  of  this  approach  are  in  the  section  titled  Lrir-.T 
(Likelihood  MATrix). 
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A  closer  examine!  nation  was  made  of  the  problem  of  bearirg 
ambiguity.  The  technique  employed  in  this  po"tion  of  tne 
analysis  is  given  in  the  section  titled  AMBIGUITY  RESOLUTION. 
This  analysis  gives  some  useful  insight  into  possible 
difficulties  within  the  antenna  system  that  may  prove  to  be 
the  most  tractable. 

5  .  DFFRR 

The  purpose  of  this  analysis  was  to  determine  new 
accurately  a  OP  bearing  could  be  calculated  fron  the  given 
data.  The  data  consists  of  23  2  seconds  of  'dtIV  per  file  (one 
file  of  KLC).  To  study  short  duration  signals  it  is  only 
necessary  to  consider  the  230  seconds  of  data  to  be  a 
continuous  concatenation  of  short  duration  signals.  Tne 
10.030  records  in  each  file  contain  the  data  for  tne  23  ms 
sampling  oeriods»  therefore,  integer  multiples  of  records 
correspond  to  different  signal  durations.  To  examine  system 
performance  against  a  230  ms  signal,  one  need  only  examine  a 
file  10  records  at  a  time.  A  203  second  file  may  be  theugnt 
of  as  containing  1003  signals  of  233  ms  duration.  Similarly, 
for  a  signal  duration  of  1  second.  53  records  may  be  used  to 
synthesize  the  signal,  and  the  file  is  made  up  of  li'o 
signals . 

A  FORTRAN  urogram  named  DFFRR  was  written  to  examine  tne 
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The  general  purpose  of  DFFPF  is  to  examine  system 
performance  as  two  parameters  are  varied.  The  first  parameter 
is  signal  duration  5  the  second  is  the  A4  term.  Tne  A4  term, 
explained  in  section  III.  is  the  measure  of  inconsistency 
within  the  D?  system.  If  the  A4  term  is  large,  tne  bearing 
value  in  a  ”ecord  is  not  considered  reliable.  If  a  4  is 
small. more  confidence  is  Placed  in  the  bearing.  The  relative 
descriptor  large  and  small  have  vet  to  be  evaluated.  In 
order  to  evaluate  the  pertinent  "ange  of  A 4  values,  the  A4 
threshold  (A4MAX )  is  varied  between  a  small  value.  3.1,  and  a 
large  value.  1.0.  in  increments  of  0.1. 
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records  a  re  being  processed. 
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immediately  rejected  because  the  value  of  A4  is  greater  than 
A4MAX.  The  remaining  four  "eco^ds  a^e  called  ”A4  admissible" 
and  are  used  to  determine  a  tearing  mean  and  standard 
deviation:  (The  mean  and  standard  deviation  formulas  used  in 
DTEP.R  are  derived  in  Appendix  A.) 

id  TAN  =  339  deg  STD  =  44  deg 

These  statistics  are  used  to  f  onn  a  window  centered  at  309 
degrees  extending  44  degrees  on  ei  the”  side  of  the  mean.  The 
A4  admissible  hearings  must  pass  through  this  window  f c ^ 
further  consideration.  Record  number  354  is  not  window 
admissible"  and  is  discarded.  The  remaining  tnree  records, 
being  both  A4  and  window  admissible,  a”0  used  to  compute  a 
second  mean  ani  standard  deviation: 

MFAN  =  33-5  *  eg  STD  -  4  deg 
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This  mean  is  DFFRR's  best  estimate  cf  the  bearing  her  this 
one  102  rs  signal  (~ecords  252  through  254  . 

This  reported  bearing  is  compared  with  the  true  Gearing . 
337  oeg,  and  the  bearing  error  is  computed  as  335-337=-2  deg. 
Additionally,  the  valid  signal  counter  is  augmented  by  one. 
The  number  of  valid  signals  will  b°  used  later  to  determine 
the  probability  of  obtaining  a  bearing  (?0:). 

If  there  are  less  than  three  A  A  admissible  records  or  less 
than  two  window  admissible  records  in  a  given  signal 
duration,  the  signal  is  considered  invalid  aue  to 
insufficient  ^  a  ta  and  is  counted  in  an  invalid  signal 
counter.  To  avoid  the  loss  of  reliable  data  in  the  case  of  a 
small  standard  deviation  of  the  A4  admissible  beamings,  the 
screening  window  is  not  permitted  to  be  narrower  than  12 
degrees . 

The  DFZRR  urogram  processes  ten  separate  signal  durations 
at  ten  different  A4MAX  settings.  Program  output  consists  cf 
four  tables  on  separate  pages.  Each  table  is  a  ten  by  ten 
matrix;  the  ’’ovs  correspond  to  the  A 4 MAX  vlaues  and  the 
columns  correspond  to  the  signal  duration,  see  tables  II 
through  '7.  The  first  table  (table  II)  is  the  average  bearing 
e'TO".  This  is  the  average  of  all  the  separate  means 
reported.  In  the  example  above  fo1’  120  ms,  the  rear.  335  would 
be  one  of  a  maximum  possible  220?  values  that  would  be 
averaged  and  then  reported  in  toe  row  A4MAX=0.2  and  the 
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column  signal  duration  equals  100  ms. 

The  second  table  (table  III)  is  the  standard  deviation  of 
the  bearing  errors  reported  in  the  first  table,  for  tne  100 
ms  signal  duration  category  the  standard  deviation  would  be 
of  a  maximum  possible  2300  mean  values  (the  actual  sample 
size  is  equal  to  the  number  of  valid  signals).  As  tne  signal 
duration  increases,  the  sample  size  decreases.  For  tne  1J03 
ms  column,  the  maximum  possible  number  of  signals  is  200.  If 
medium  duration  signals  are  examined  with  TFiP.R.  one  must  be 
attentive  to  the  sample  size.  For  a  10  second  duration  signal 
category,  the  sample  size  has  diminished  to  twenty.  For  103 
second  signal  durations ,  there  are  only  two  samples  and  tne 
validity  of  a  standard  deviation  is  highly  questionable. 

The  third  table  of  output  (table  17)  is  tne  average 
intra-signal  standard  deviation.  This  is  tne  average  of  tne 
standard  deviations  reported  for  individual  signals.  In  the 
previous  example.  the  3TD=4  would  be  one  of  a  maximum 
possible  2300  standard  deviations  to  be  averaged.  If  tne 
average  intra-signal  standard  deviation  were  equal  to  four, 
the  interpretation  would  be  that  for  all  signals  of  100  ms 
duration,  after  the  unreliable  bearings  are  discarded.  the 
expected  standard  deviation  of  the  remaining  cluster  of 
bearings  is  four  degrees.  The  term  intra-signal  is  used  ic 
distinguish  it  from  the  standard  deviation  of  bearing  errors. 

The  fourth  table  (table  V)  is  the  compilation  of  valid 
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Table  II 

DFERR  Output  Page  1 ,  Average  Bearing  Error 


OSTANDARD 

DEVIATION 

OF 

BEARING 

ERROR 

AS  A 

FUNCTION  OF 

SYSTEM 

NOISE 

AND 

SIGNAL 

.1.  .  0 

33 

30 

26 

24 

24 

2  p 

:i  9 

20 

15 

19 

0.9 

33 

30 

26 

24 

24 

22 

19 

20 

15 

19 

M 

0 . 8 

33 

30 

26 

24 

24 

22 

19 

20 

1 5 

1 9 

A 

0.7 

33 

30 

26 

24 

25 

22 

19 

20 

15 

18 

X 

0.6 

33 

30 

26 

23 

25 

23 

18 

0  7 

A#..  \J 

17 

18 

0 . 5 

31 

29 

27 

24 

25 

23 

20 

23 

16 

17 

A 

0.4 

?9 

28 

26 

24 

23 

21 

20 

21 

17 

17 

4 

0.3 

29 

27 

26 

25 

23 

1 9 

20 

21 

18 

17 

0.2 

29 

31 

■ , ,;) 

W  A... 

31 

28 

24 

27 

28 

2  2 

22 

o .  l 

20 

22 

25 

25 

28 

23 

24 

pp 

23 

22 

100 

200 

300 

400 

500 

600 

700 

800 

900 

1000 

SIGNAL  DURATION  ( MILL I SEC ) 

Table  III 

DFERR  Output  Page  2 ,  Standard  Deviation  of  Bearing  Error 
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Table  IV 

DFERR  Output  Page  3 ,  Intra-Signal  Standard  Deviation 
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DFERR  Output  Page  4,  Number  of  Valid  Signals 


signals  for  each  of  the  A4MAX  and  signal  duration  categories. 
This  is  an  important  statistic  needed  to  compute  the 
probability  of  obtaining  a  bearing.  In  the  case  of  the  102  ms 
signal  duration,  2200  are  possible,  if  A4MAX=v5.2  and  1200 
signals  are  valid,  the  FOB  for  100  ms  is  equal  to  the  ratio 
of  valid  signals  to  possible  signals.  In  this  case.  P0S=3.c. 
In  general,  as  either  A4MAX  or  tae  signal  duration  increases, 
the  number  of  valid  signals  approaches  the  number  of  possible 
signals.  For  A4MAX  above  0.4  or  the  signal  duration  above  one 
second,  the  ?0B  is  approximately  one. 

DFEFP  was  used  to  process  all  of  the  files,  including  tne 
OCX  modulated  XLC  signal.  The  XLC  file  required  a  slightly 
modified  version  of  PFFRR  because  each  record  of  data  was 
searched  for  a  flag  that  indicates  that  the  signal  is  indeed 
present.  Search  was  also  made  fo^  flags  that  indicate  that 
tne  system  is  saturated.  XLC  is  considered  a  file  of  signals 
separated  by  noise.  It  is  not  a  continuous  signal  like  the 
WWV  signal?  in  fact,  the  duty  cycle  is  less  than  twenty-five 
percent . 

The  tabular  output  produced  by  DFSRR  can  be  considered  a 
second  level  data  base.  To  simplify  follow-on  discussions, 
the  statistics  in  the  DFSRR  output  will  be  referred  to  as  the 
FOUR'  statistics.  Examination  of  the  DFSRR  data  base 
revealed  that  the  FOUR  statistics  are  a  strong  function  of 
the  A4  term.  As  the  A4  term  increases,  corresponding  to  tie 
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acceptance  of  more  inconsistent  data,  the  standard  deviation, 
intra-signal  standard  deviation  and  number  of  valid  signals 
increases.  Based  on  the  large  standard  deviations  recorded 
above  A4=0.4  and  the  insignificant  increase  in  P0£  above 
A4=0.4,  the  range  of  interest  was  restricted  to  A4  values  in 
the  0.1  to  0.4  region.  Even  with  this  restriction,  the  amount 
of  data  is  too  large  to  present  in  this  report;  however,  it 
should  be  noted  that  the  changes  in  the  FOUR  statistics  are 

typically  monotonic  as  A4  varies.  The  data  presented  herein 

is  for  A4MAX=0.2  and  A4MAX=2.4.  The  best  nerformance 

category.  A4MAX=0.1.  is  not  presented  because  of  its  low  PCB 
and  because  it  is  not  significantly  different  from  A4rtAX=3.2 
The  tabular  data  of  DFERR  does  not  permit  easy  visual 
perception  of  the  characteristics  of  the  FOUR  statistics.  A 
plot  program  was  written  to  display  the  data.  Tne  figures  at 

the  end  of  this  section  are  of  the  V  5.10  and  20  MHz  files 

of  1979  (Fig.  10-24)  and  1980  (Fig.  25-34).  For  eacn 
frequency  there  are  five  graphs.  In  eacn  set  of  five,  tne 

first  two  and  the  last  are  the  most  important  as  tney  are 

concerned  with  the  short  duration  signals.  The  second  two 

graphs  are  for  signal  durations  of  one  to  ten  seconds.  They 

are  included  to  illustrate  how  the  FOUR  statistics  tend  to 
reach  steady  state  values.  The  fifth  graph  is  the  histogram 
of  the  tearing  errors  for  A4MAX=0.2  ar.d  signal  duration  equal 
to  200  ms.  Annotated  in  the  upper  right  hand  corner  of  all 
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the  compute”  drawn  graphs  are  the  average  hearing  (3  . 
standard  deviation  (-)  and  average  hearing  error  (  TJ  for  the 
entire  file  portrayed  by  the  graph.  These  are  tne  statistics 
for  a  signal  duration  of  200  seconds  subject  to  tne  A4MAX 
constraint  labelled  on  the  graph.  If  one  assumes  tnat  the 
true  angle  of  arrival  was  fixed,  that  no  multipath  existed 
and  that  there  was  no  slow  tern  variance  due  to  tne 
ionosphere,  the  curves  plotted  should  converge  to  tne  200 
second  average  bearing  error.  Both  the  standard  deviation  and 
the  average  intra-signal  standard  deviation  should  converge 
to  the  20?  second  standard  deviation.  The  P03  should  converge 
to  one.  However,  the  PROPHET  program  graphs  show  that  in  many 
cases  there  probably  existed  multimode  conditions  (actual 
amount  of  interference  unknown )  and  that,  in  all  cases, 
varying  degrees  of  variance  existed  due  to  polarization 
fading  and  interference.  Despite  this  difficulty.  the  200 
second  statistics  can  be  considered  an  approximate 
convergence  point. 

The  first  curve  (X)  of  each  graph  is  the  average  bearing 
error.  This  is  the  difference  between  the  calculated  bearing 
and  the  true  geogranhic  bearing  to  the  signal  transmitter 
(337  or  090).  The  most  prominent  feature  of  all  the  average 
bearing  error  curves  is  that  tney  are  not  a  strong  function 
of  signal  duration.  This  could  have  been  anticipated 
realizing  that  the  average  of  the  majority  of  small  subsets 
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of  a  large  set  will  closely  approximate  the  average  of  tne 
large  set.  In  this  respect  the  average  tearing  error  as  a 
function  of  signal  duration  is  not  particularly  useful.  A 
much  more  significant  view  of  the  average  bearing  error  is  a 
histogram  of  the  bearing  errors  from  each  valid  signal.  A 
histogram  for  A4MAX=0.2  and  a  signal  duration  equal  to  200  ms 
is  the  fiftn  figure  in  each  set  of  five  (Fig. 
14,19,24,29,34).  These  histograms  show  that  the  distribution 
of  the  bearing  error  for  most  of  the  files  is  only  rougnly 
approximate  to  a  normal  distribution,  .^ajor  deviations  from 
the  normal  curve  are  the  accumulation  cf  bearings  in  tne  120 
degree  ambiguity  region  and  the  distributions  where  multimode 
propagation  was  highly  probable.  The  major  difficulty  with 
the  average  bearing  e^ror  data  is  the  absence  of  calibration 
data.  Tae  amount  of  correctable  bias  error  is  unknown. 

The  standard  deviation  curve  (*)  is  tne  most  important 
information  on.  the  graphs.  If  a  bearing  is  to  be  used  with 
other  bearings  to  compute  a  fix,  the  standard  deviation  is 
used  to  compute  the  fix  area  for  a  given  probability  that  the 
target  will  be  within  the  fix  area.  The  standard  deviation  is 
also  a  measure  of  confidence  in  a  single  bearing.  If  the 
bearing  distribution  is  normal  and  the  standard  deviation  is 
20  degrees,  one  can  expect  that  the  bearing  calculated  is 
within  20  degrees  of  the  true  bearing  (assuming  average 
bearing=true  bearing)  about  57  percent  of  tne  time.  Because 
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the  distributions  for  the  data  files  are  only  roughly  normal, 
amplification  of  the  significance  of  the  standard  deviations 
is  pertinent.  The  following  data  are  the  approximate 


percentages  (within  five  percent)  of  bearings  falling  within 

‘or  a  2O0  ms  signal 


the  standard  deviation 

of  each 

file 

duration  and  A4MAX='0.2 

• 

File  ft 

STD 

0/ 

o 

1  (Fig.  19) 

31.1 

90 

2  (Fig  .  24  ) 

25.6 

3  5 

■7 

31.1 

75 

4  (Fig.  14) 

15.5 

55 

5 

37.4 

90 

10  ( Fig  .29 ) 

59.6 

75 

11  (Fig.  34 ) 

29.5 

95 

12  (Fig.  36 ' 

16.2 

88 

13 

19.  3 

83 

This  data  indicates 

that  the  < 

i is  t  i 

the  normal  and  that  system  performance  is  better  than  one. 
thinking  in  terms  of  a  normal  distribution,  would  believe. 
For  example,  in  file  11  (Fig.  34)  the  standard  deviation  is 
about  30  degrees,  but  95  percent  of  tne  data  are  within  tne 
standard  deviation.  The  expression  of  standard  deviation 
cannot  be  separated  from  its  distribution  and  still  retain 
meaning . 

Observing  the  curves  (Fig.  10-13.  15-18,  cd- 23,  25-28, 
30-33)  the  reader  will  note  that  the  standard  deviation  is  in 
almost  all  cases  a  monotoni cally  decreasing  function  of 
signal  duration.  An  interesting  comparison  is  the  standard 
deviation  at  100  ms  (STD1)  and  at  10  seconds  (STD2),  again 
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with  A4MAX=P .2 : 


File 

*  5TD1 

5TD2 

STD1/STD2 

1 

29 

5 

5.3 

2 

25 

11 

2.3 

3 

31 

13 

2.1 

4 

16 

5 

3.6 

5 

42 

10 

4.2 

13 

57 

16 

3.2 

11 

30 

A 

7.5 

12 

16 

9 

1  .8 

13 

17 

6 

2.9 

These  results  are  comparable  to  those  reported  by  Eair. 
(Ref.6j  and  consistent  in  order  of  magnitude  to  the  time 
required  to  average  the  fluctuating  surfaces  of  constant 
ampli tude . 

The  third  curve  (.)  is  the  average  intra-signal  standard 
deviation.  It  is  typically  a  mono t cni cally  increasing 
function  of  signal  duration  from  12C  ms  to  12  seconds.  The 
fact  that  the  intra-signal  standard  deviation  is  the  smallest 
at  130  ms  indicates  that  the  major  factors  affecting  tne 
variance  are  not  rapidly  fluctuating,  i.e.  the  period  is 
larger  than  102  ms.  The  intra-signal  standard  deviation  at 
122  ms  is  small,  four  to  eight  degrees,  except  for  file  12 
(Fig.  25-25)  which  was  recorded  with  a  low  SNR.  In  tnis  case, 
a  major  source  of  variance  was  noise  and  its  fluctuations 
were  rapid.  As  the  signal  duration  increases,  the  pnenomena 
causing  the  majority  of  the  variance  have  mere  effect  within 
individual  signals. 

The  fourth  category  of  data  on  the  graphs  is  the 
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probability  of  obtaining  a  line  of  tearing.  The  derivation  of 
this  number  has  previously  been  explained.  Its  usefulness  is 
a  proper  subject  for  operations  research;  it  does  not  provide 
the  engineer  with  any  useful  information  about  t ne  process  of 
HFDF  on  short  duration  signals.  However,  one  does  not  need  a 
specific  operational  context  to  know  that  hign  POT  is  good 
and  low  FOB  is  bad. 

One  can  summarize  this  section  by  stating:  (1)  A4  is  an 
effective  measure  of  reliability  of  tne  data,  (2)  the 
variance  is  high  for  snort  duration  signals,  (3;  tne  average 
bearing  error  is  not  useful  without  calibration  data  and  (4) 
variance  improvement  with  time  averaging  corresponds  in  order 
of  magnitude  with  that  predicted  bv  fading1  phenomena. 
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Figure  10 

WWV  5  MHz  2/79  Short  Signal  Duration  A4MAX=0.2 
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Figure  11 

WWV  5  MHz  2/79  Short  Signal  Duration  A4MAX=0.4 


SOURCE:  WWV  5M  7 : C  C  2/7S  A4MAX=0.4 
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Figure  12 

WWV  5  MHz  2/79  Medium  Signal  Duration  A4MAX=0.2 


B=  3  3  5  .  5 

SOLRCE  :  UV4  V  5N  7:00  2/  79  A4I^AX  =  0.2 

a  =7 . 3 

AVE  BE/RING  ERROR  (DEGREES)  :  X  B 

ST  C  OF  BEARING  ERPCR  (DEGREES)  :  *  e  =  1.5 

AVE  CF  INTRA-SIGNAL  STC  (DEGREES)  :  . 


11  .CCC 


3.  5CCG 


♦ 

I 

I 

I 

I 

I 

I 

I 

I 

1 

<• 

I 

I 

I 

I 

I 

I 

1 

I 

I 


-4.C000  ♦ 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X  + 

- +  - 

0.0 

SIGNAL  CLRATICN: 

1CCC  . 

2000. 

3000. 

40  00  . 

5000. 

5000. 

60  OC  • 

7CCC. 

8000. 

5000 

10000  . 

.10GCC  . 

(  MIL  ISECCNCS  ) 

PROBABILITY  CE  : 

1  .000 

1  .000 

1  .0  00 

l.OOC  l.CCO  1 

.000 

L.000 

1.000  ] 

L  .  CCC 

1. 000 

OBTAINING  LOB 

X  : 

-4  . 

-4. 

-4. 

-4. 

-  4. 

-  4. 

-4. 

-4. 

-4. 

-4  . 

4  : 

7. 

7. 

t. 

6. 

6. 

6. 

6. 

5  . 

5. 

5. 

•  • 

9. 

10. 

1  1. 

11. 

11. 

1  1  . 

11. 

11. 

1  1. 

11. 

Figure  13 

WWV  5  MHz  2/79  Medium  Signal  Duration  A4MAX=0.4 
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WWV  5  MHz  2/79  Bearing  Error  Histogram 


Figure  15 

WWV  10  MHz  2/79  Short  Signal  Duration  A4MAX=0.2 
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Figure  16 

WWV  10  MHz  2/79  Short  Signal  Duration  A4MAX=0.4 
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Figure  17 

WWV  10  MHz  2/79  Medium  Signal  Duration  A4MAX=0.2 


SOURCE :  M>V  10P  8:30  2/79  A4RAX=0.2 
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Figure  18 

WWV  10  MHz  2/79  Medium  Signal  Duration  A4MAX=0.4 
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WWV  10  MHz  2/79  Bearing  Error  Histogram 


Figure  20 

WWV  20  MHz  2/79  Short  Signal  Duration  A4MAX=0 . 2 
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Figure  21 

WWV  20  MHz  2/79  Short  Signal  Duration  A4MAX=0.4 
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C.  LMAT 


The  previous  section  detailed  the  results  of  accepting 
data  record  hearings  on  the  tasis  of  the  A4  term.  The 
calculated  variances  are  large  and  in  some  cases  the 
corresponding  POB  is  small.  There  is  a  strong  need  to  lower 
the  variances  and  increase  the  POB  to  approximately  one.  To 
do  this  requires  a  mere  complete  use  of  the  data;  this 
includes  the  A?  and  Phase  terms. 

The  theoretical  equations  of  the  coaxial  spaced  loons 
prove  the  existence  of  polarization  independent  nulls  in  the 
ideal  case.  However,  the  construction  of  perfect  loops  and 
free  space  siting  of  these  loops  is  clearly  impossible. 
Recognizing  that  the  spaced  loop  array  does  not  perform 
ideally,  one  attempts  to  identify  the  sources  of  fixed  error 
and  correct  them  with  a  calibration  data  set.  Otner  sources 
of  error  are  random  and  are  described  in  probabilistic  terms. 
It  is  usually  the  case  that  more  than  one  of  the  major  errors 
is  random.  This  results  in  joint  probability  density 
functions  that  are  impossible  to  derive  analytically  and 
impossible  to  isolate  in  order  to  measure. 

In  an  attempt  to  maximize  the  use  of  information  available 
without  a  precise  knowledge  of  either  the  fixed  or  random 
errors,  it  was  decided  to  take  a  decision  theory  approach.  An 
assessment  of  a  bearing's  reliability  is  what  is  most  needed. 
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The  decision  to  he  made  is  binary;  a  hearing  is  reliable  or 
unreliable.  Once  the  bearing  reliability  is  determined,  tne 
statistical  procedure  developed  for  DFERR  can  be  used  to 
determine  a  mean  bearing.  The  decision  to  determine  if  a 
bearing  is  reliable  will  be  based  on  the  likelihood  that  it 
i s  rel iable . 

As  an  example,  assume  the  following  data  records: 


Record  # 

A0 

A4 

Bearing 

1 

0  .113 

0.214 

342 

2 

0.612 

0  .208 

358 

3 

0.514 

0.421 

288 

4 

3.815 

•  *  • 

0.113 

*  •  » 

305 

•  •  « 

bearing  is 

known 

to  be  337. 

Assume  that  a  se 

is  recorded  on  an 

unloca  ted 

target  transmitter: 

Record  # 

A0 

A4 

Bearing 

1 

0.056 

0.298 

105 

2 

0.822 

0.173 

042 

3 

0.109 

0  .202 

093 

4 

0 . 666 

0.432 

087 

ccrd  is 

most  re 

liable?  It 

i s  known  in  data  set 

that  the  first  record  is  tne  most  reliable  because  its 
bearing  is  the  most  accurate.  Using  the  maximum  lik'el  incod 
(ML)  criterion,  one  would  estimate  that  the  record  in  data 
set  two  that  most  closely  resembles  record  one  of  data  set 
one  is  the  most  reliable.  By  ML  criterion  record  tnree  would 
be  selected  because  its  A0  and  A4  parameters  are  the  closest 
match  to  record  one  of  the  first  data  set.  In  this  scheme 
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there  is  no  filtering  on  preset  A0  or  A4  limits.  If  tne 
lowest  level  of  A0  and  A4  were  the  reliabilty  criterion, 
record  one  of  the  second  data  set  would  have  been  selected  as 
the  best  estimate.  The  advantage  of  the  ML  criterion  is  tnat 
it  is  a  "use  what  works"  technique.  The  deductive  approach  of 
analysis  of  the  system  and  errors  is  discarded  because  it  is 
too  complicated.  Instead,  the  inductive  technique  of 
observing  and  classifying  provides  the  more  attractive 
approach  in  this  case. 

To  use  the  ML  criterion  with  the  StoRI  data  files  one  must 
expand  upon  the  ideas  presented  in  the  example.  Define  two 
matrices,  a  bearing  acceptable  and  a  bearing  unacceptable 
matrix.  3ach  matrix  is  three  dimensional;  a  dimension  is 
allotted  to  k?> ,  one  to  Phase  and  the  third  to  A4.  This  spans 
all  the  information  in  a  data  record.  A0  is  a  measure  of  the 
horizontal  field  component;  Phase  is  a  measure  of 
polarization;  and,  A4  is  a  measure  of  system  inconsistency. 
These  three  measures  are  not  sufficient  to  completely  specify 
system  performance?  in  fact,  they  are  not  sufficient  to 
completely  specify  the  horizontal  field  component, 
polarization  or  system  inconsistency .  Tney  are  what  is 
available.  Pach  dimension  is  divided  into  ten  increments.  The 
A?  and  A4  dimension  values  are  between  0  and  1.0  in 
increments  of  0.1  .  The  Phase  dimension  is  between  -183  and 
+1S0  in  increments  of  36  degrees.  The  structure  of  both 
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matrices  can  be  sketched  as: 


A0 


^  j - 1 - 1 - 1 - r- 


"t — r 


PHf^SE 


M 


I 

O 


0.5 


.  r  i  '  '  1 
~\80  O 


0.5 


1.0 


' I 
\80 


1.0 


The  matrices  are  defined  to  be  expressions  of  tne 
conditional  probability  mass.  They  are  created  by  using  a 
file  or  files  of  data  on  a  target  of  known  position.  Knowing 
the  true  bearing  permits  one  to  address  either  the  bearing 
acceptable  or  bearing  unacceptable  matrix.  Further 
explanation  is  best  presented  using  an  example.  The  following 


record  s 

are 

available  from 

a  data  set 

with  a 

true  bearing 

known  to 

be 

337  deg: 

Re 

cord  #  A0 

Pha  se 

A4 

3ear ing 

( deg ) 

(deg; 

i 

0.911 

-43 

0.621 

266 

2 

0.813 

—  64: 

0  .518 

273 

3 

0.517 

120 

0.231 

357 

4 

0.342 

72 

0.165 

340 

All 

the 

elements 

of  the 

bearing 

acceptable  and 

unacceptable 

matrices  are 

initialized 

at  zero  . 

The  criterion 

99 


for  placement  in  one  matrix  or  tne  other  is  the  value  of  the 
bearing.  A  window  about  the  true  bearing  is  defined.  If  tne 
window  is  defined  as  1?  degrees,  a  bearing  is  acceptable  if 
it  is  within  5  degrees  either  side  of  the  true  bearing.  Eacn. 
record  in  a  file  is  examined.  Record  number  one  is 
unacceptable  because  the  beaming  value  is  not  in  the  232  o  2 
deg  window.  Therefore,  its  A 0,  Phase  and  A4  values  are  mapped 
into  increments  along  the  respective  dimensions  of  the 
unacceptable  matrix,  thereby  addressing  one  of  the  elements 
of  the  matrix.  A  one  is  added  to  the  contents  of  tnis 
element.  Each  record  in  turn  is  examined  .and  a  one  is  added 
either  to  an  element  in  the  bearing  acceptable  or  the  bearing 
unacceptable  matrix.  (Record  number  four  is  an  example  of  a 
record  that  would  apply  to  the  bearing  acceptable  matrix./ 
After  all  records  have  been  processed,  each  element  in  the 
bearing  acceptable  matrix  is  divided  by  the  number  of  records 
that  contributed  to  that  matrix.  Similarly,  the  elements  of 
the  bearing  unacceptable  matrix  are  divided  by  the  number  of 
records  that  contributed  to  it.  This  produces  an  expression 
of  conditional  probability.  If  a  bearing  is  acceptable,  the 
probability  of  it  having  a  particular  A0,  Phase  and  A4  value 
is  equal  to  the  value  in  the  acceptable  matrix  addressed  by 
the  given  A0 ,  Phase  and  A4  values.  Likewise,  if  a  bearing  is 
unacceptable,  the  probabilities  of  its  A0.  Phase  and  A4 
values  can  be  read  from  the  element  addressed  by  those 
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values  . 


Tne  acceptable  and  unacceptable  natrices  are  constructed 
from  data  on  a  located  target.  If  one  processes  a  file  of 
data  on  an  unknown  target,  the  natrices  are  used  as  follows. 
Using  the  A 2,  Phase  and  A4  values  to  address  bcth  of  the 
natrices.  the  element  values  of  the  natrices  are  compared. 
Suppose  the  value  from  the  acceptable  matrix  is  9.211.  and 
from  the  unacceptable  matrix,  it  is  0.997.  The  ML  criterion 
makes  the  decision  for  the  matrix  with  the  nignest 
probability  (maximum  likelihood).  3.211  is  the  highest  value; 
therefore,  the  decision  is  that  the  bearing  is  acceptable. 

The  requirement  to  address  two  matrices  and  compare  the 
returned  values  can  be  eliminated  by  forming  a  single 
likelihood  ratio  matrix.  This  matrix  is  constructed  by 
dividing  each  element  of  the  bearing  acceptable  matrix  by  the 
corresponding  element  of  the  bearing  unacceptable  matrix.  The 
decision  can  be  made  on  a  single  element  of  tne  likelihood 
ratio  matrix  (also  simply  called  likelinood  matrix  or  denoted 
as  [L]  ) .  In  the  example  above,  the  ratio  of  the  two  values  is 
2.175.  The  decision  is  that  a  bearing  is  acceptable  if  the 
addressed  element  of  the  likelihood  matrix  is  greater  than  or 
equal  to  one.  If  less  than  one,  the  bearing  is  unacceptable. 

The  above  ideas  car  be  set  into  mathematical  notation  as 
follows.  Matrices  are  denoted  with  brackets.  The  elements  of 
a  matrix  are  indexed  with  i,j,.-c  in  the  general  case. 


101 


[A]  :  bearing  acceptable  matrix 

fTJ]  :  bearing  unacceptable  matnix 
[A  (i,,i,k)  ]  =  (A  (  A0  ,Pha  se  .  A4 )  ] 

fTJ  (  i  ,  j  ,  k )  1  =  fu  ( A0  ,  Phase  ,  A4 )  ] 

Na  =  number  of  records  tabulated  in  l A ] 

Mu  =  number  of  records  tabulated  in  [U] 

Na  +  Mu  =  (10000)  X  (number  of  files  used) 

The  probability  mass  matrices  are: 


[Pal  =  ( 1  /Na  )  T Al 
(Pul  =  ( 1  /Nu  )  f U 1 

The  maximum  likelihood  (ML)  criterion  can  be  used  to 
decide  the  acceptability  of  a  beaming, 
let  the  record  be  fA3=<*, ?hase=.^,A4='£.Bo } .  then, 

DL 

A  <  U 

D\  “'PS* 


where,  the  decisions  31  and  32  are. 

DP:  Pearing  is  acceptable 
31:  Searing  is  unacceptable 

If  A  >  U,  the  decision  is  DP?  if  A  <  U.  the  decision  is  Dl. 

This  comparison  process  can  be  simplified  and  at  tne  same 
time  be  made  more  flexible  by  defining  a  livelihood  ratio 
matrix  (Li  : 
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(LI  =  f  Pa]  ©  (Pul 


where  the  symbol  (©)  means  to  divide  each  element  of  the 
matrix  [Pa]  only  by  its  i.j.k  counterpart  in  (Pul.  (Define 


x/0  =  infinity.  ) 

The  ML  criterion  can  be  rewritten  as: 

D2 

[L]  <  [1]  [1]  =  12x12x12  matrix  of  1  's 

Dl 

o 

OZ 


The  above  idea  can  be  extended  to  the  Minimum  Probability  of 
Error  (MPE)  and  the  3ayes  Cost  (3C)  criteria 
s  traig’ntf  orwardly : 


m 


D.2 

> 

< 

Dl 


P{3  acceptable} 


P(B  unacceptable} 


Dl 


<<n,  -  c„  > 

P f B  acc } 

?{B  unacc} 

C 

C 

c 

c 


21  * 

• 

1 1  • 
12  *• 
22.  : 


cost  of 
cost  of 
cost  of 
cost  of 


deciding 

deciding 

deciding 

deciding 


acceptable  when  unacceptable 
unacceptable  when  acceptable 
unacceptable  when  unacceptable 
acceptable  when  acceptable 


The  terms  on  the  right  side  of  the  decision  symbols  are 
matrices.  The  ratio  of  P{3acceptable}  to  P( Bunacceptable} 
must  be  determined  separately  for  each  element  of  the  matrix. 
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The  sane  will  probably  be  true  for  the  Bayesian  costs.  The 
MPE  and  EC  criteria  are  nore  sophisticated  and  will  probably 
yield  better  results,  but  data  limitation  prevented  further 
investigation.  Very  large  amounts  of  data  would  be  needed  to 
determine  probability  ratios  for  each  element  of  a  10x10x10 
matrix.  It  would  also  require  a  operational  input  for  a 
realistic  assessment  of  Bayesian  costs. 

Using  the  above  concepts.  the  FORTRAN  program  LMAT  was 
written  to  produce  the  li-celinood  ^atio  matrices  iLj  f^om 
individual  files  and  from  groups  of  files.  Testing  of  tine 
method  centered  on  the  15  MHz  files  because  they  are  the  most 
numerous  for  a  single  frequency  (three  files).  An  L  matrix 
was  created  from  file  3  for  two  separate  definitions  of 
"bearing  acceptable".  The  first,  denoted  L[3.S0j.  considered 
the  bearing  acceptable  window  to  be  ninety  degrees  wide.  The 
second,  Lr3.20).  used  a  twenty  degree  window  to  construct  tne 
matrix. 

Using  the  L  matrices,  several  data  files  were  processed 
employing  the  program  LFILE  (L  matrix  modified  FILE).  This 
program  reads  the  A0 ,  Phase  and  A  4  term  of  each  record  and 
uses  them  to  address  the  L  matrix.  If  the  element  addressed 
is  greater  than  one.  the  bearing  in  that  record  is  considered 
acceptable,  and  the  record  is  written  into  a  new  data  file 
with  a  new  A4  term  equal  to  0.01.  Only  the  A4  term  is 
changed.  The  new  data  file  is  subsequently  processed  by  BFERR 
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which  will  always  treat  an  acceptable  bearing  (A4=0.01)  as 
valid  for  statistical  processing.  If  the  L  matrix  element 
addressed  is  less  that  one.  the  bearing  is  unacceptable  and 
assigned  the  A4  value  10. £  in  the  new  data  file.  DFiRP.  will 
reject  any  bearing  with  this  large  A4  value. 

The  file  of  irterest  is  the  1990  15  MHz  file,  figure  35  is 
the  result  of  simple  DFEP.R  processing  on  this  file;  figure  36 
is  the  histogram  of  the  bearing  error.  Several  examples  of  L 
matrix  processing  on  this  file  a -e  graphed  in  figures  37,  59 
and  41.  The  explanation  of  these  graphs  is  identical  to  the 
explanation  given  in  the  previous  section.  It  will  be  noted 
.however,  that  the  graphs  are  labelled  with  the  L  matrix 
notation.  The  A4!^AX  term  is  not  applicable.  Each  grapn  has  an 
accompaning  histogram  of  the  bearing  erno".  figure  37  is  tie 
file  processed  by  an  L  matrix  made  from  the  file  itself. 
Figure  59  is  the  file  processed  bv  an  L  matrix  created  from  a 
1973  file  at  the  same  frequency.  Figure  41  is  the  file 
processed  by  an  L  matrix  created  f”om  f oim  1979  4 s V  files  at 
different  frequencies.  Mot  unexpectedly,  the  best  performance 
is  by  the  L  matrix  formed  from  the  file  itself  (Fig.  37).  But 
it  is  interesting  to  note  that  the  L  matrix  created  from  the 
single  file  of  identical  freauency  (Fig.  39)  performed  better 
than  the  matrix  made  up  from  the  four  files  (Fig.  41).  This 
tends  to  confirm  frequency  sensitivity  in  the  L  matrix.  It 
may  prove  valuable  to  add  a  frequency  dimension,  making  the  L 
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matrix  four  dimensional. 


A  considerable  amount  of  processing  was  devoted  to 
determine  an  optimum  window  for  the  likelihood  matrix.  No  one 
window  width  outperformed  other  widths  in  all  categories.  The 
central  difficult?'-  is  that  one  thousand  elements  in  the  L 
matrix  must  be  determined.  Accuracy  in  terms  of  ML  for  each 
element  is  a  limiting  process;  as  the  number  of  bearings  used 
to  determine  eacn  element  apprcac.oes  infinity,  the  true  ML 
ratio  is  determined.  To  approximate  the  infinite  sample  ML 
ratio  within  1?  percent  would  require  very  rougnly  120 
records  per  element.  But  some  elements  are  rarely  addressed; 
therefore,  large  amounts  of  data  are  needed  to  fill  tne  L 
matrix.  The  amount  of  data  can  be  estimated  in  a  very  rougn 
vay.  Using  four  Mies  and  a  ninety  degree  window,  it  was 
observed  that  about  40  percent  of  the  ratrix  elements  are 
addressed  zero  or  one  times.  If  the  requirement  is  that  the 
probability  of  this  happening  should  be  less  than  0.201.  tne 
expression  (0.4)  =(€.001)  determines  the  number  (n)  of  sets 
of  four  files  required.  In  tnis  case  about  thirt?/-  files  would 
be  required.  This  corresponds  to  300.000  records  or  100 
minutes  of  data.  This  is  an  easily  achieved  number  at  the 
antenna  system  development  site. 

The  most  overall  successful  runs  were  made  with  a  ninety 
degree  window.  In  table  VI  are  the  results  of  processing  the 
I960  data  files  with  a  likelihood  matrix  constructed  from  the 
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1979  data.  The  first  four  1S79  files  were  used  ;  the  fiftn,  a 
second  15  MHz  file,  was  not  included  to  avoid  a  possible 
dominant  15  MHz  bias.  The  data  in  table  VI  is  the  variance 
and  POB  from  simple  DFEP?  processing  ( A4MAX=0 .2 )  and  f^on  L 
matrix  processing  (L (1-1,90) ).  For  the  193d  5  MHz  file  the 
variance  is  high  in  either  case.  Starting  at  333  ms  signal 
duration,  the  L  matrix  processing  has  a  lower  variance,  and 
as  the  signal  duration  increases.  L  matrix  processing  is 
better  and  tetter  compared  to  DFEP.R  processing.  Peraaps  mucn 
more  significant  is  that  the  L  processing  has  a  much  higher 
PGB.  .Vith  the  1980  19  MHz  file.  L  matrix  processing  is 
superior  at  every  signal  duration.  This  is  also  true  for  the 
15  MHz  file.  This  is  a  significant  result.  L  matrix 
processing  is  generally  superior  to  the  A4  filter  process  of 
DFERR.  In  some  cases  the  variance  is  halved  and  the  ?C3  is 
more  than  doubled.  However,  performance  against  2LC  is  poorer 
as  can  te  seen  in  the  table.  For  tne  shortest  signal 
durations.  the  L  matrix  variance  is  higher  than  DFEFR 
processing.  but  tne  two  approach  each  other  past  503  ms.  Tne 
P03  is  slightly  higher  for  the  L  matrix  process. 

The  results  in  table  VI  are  very  encouraging.  A  matrix 
created  from  1979  data  significantly  outperformed  DFE  ?.?. 
processing  on  the  1989  data.  Tne  poorer  results  for  tiLC  are 
not  surprising.  a.v'V  bears  337  deg  from  San  Antonio,  and  K^C 
bears  999  deg.  The  antenna  array  which  is  mounted  atop  a 
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hcrizontallv  girded  building  suffers  pattern  distortion  due 
to  antenna  element  coupling  with  the  building.  The  distortion 
to  the  pattern  is  dependent  both  on  the  azimuth  and  elevation 
angle  at  which  the  signal  arrives.  In  fact,  the  XLC  results 
are  encouraging.  The  fact  that  L  matrix  processing  ana  DF2RR 
processing  results  are  so  close  suggests  tnat  azimuthal 
dependence  is  not  an  overiding  factor.  It  may  be  possible  to 
achieve  acceptable  results  using  L  matrices  created  from 
azimuth  sector  data.  If  a  first  guess  is  tnat  a  bearing  is  in 
the  sector  213  to  27'7  ,  post  processing  could  use  an  L  matrix 
specifically  cheated  for  that  50  degree  sector.  In  all  six  L 
matrices  would  be  required  tc  cove’’  the  full  azimuth. 

The  optimism  expressed  above  must  be  tempered  by  the  fact 
that  the  technique  and  the  results  were  derived  with  an 
inductive  approach.  A  theoretical  basis  with  which  tne  L 
matrix  success  could  have  been  predicted  was  not  derived.  L 
matrix  success  was  observed,  not  predicted.  Confidence  that 
the  technique  is  functional  in  the  general  case  will  require 
observation  in  situations  whe^e  frequency,  azimuth.  SNR. 
vertical  angle  of  arrival,  polarization,  shipboard  siting  and 
propagation  modes  a^e  varied. 
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TABLE  71 


Comparison  of  DFEP.B  and  L  'Matrix  Processing 
For  all  DEEP.?  processing,  A4MAX=0.2.  For  all  L  matrix 
processing,  the  L  matrix  is  L (1-4 ,92). 
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Figure  35 

WWV  15  MHz  2/80  Short  Signal  Duration  A4MAX=0.2 
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Figure  36 

WWV  15  MHz  2/80  Bearing  Error  Histogram  A4MAX=0.2 


Figure  37 

WWV  15  MHz  2/80  L(12 ,90) 
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Figure  38 

WWV  15  MHz  2/80  L(12,90)  Bearing  Error  Histogram 


Figure  39 

WWV  15  MHz  2/80  L(3}90) 
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Figure  40 

WWV  15  MHz  2/80  L(3,90)  Bearing  Error  Histogram 
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Figure  41 

WWV  15  MHz  2/80  LCl-4,90) 
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Figure  42 

WWV  15  MHz  2/80  L(l-4,90)  Bearing  Error  Histogram 


D.  AMBIGUITY  RESOLUTION 


The  histograms  of  the  average  bearing  e"ror  demonstrate 
one  of  the  difficulties  with  toe  SWFI  spaced  loop  antenna 
system  in  trying  to  DF  shoT't.  duration  signals.  It  can  be  seen 
by  looking  at  the  histograms  that  the~e  are  bearing  error 
numbers  that  cluster  in  relatively  large  values  at  points  of 
1S0  and  90  degree  ambiguity.  This  corresponds  to  the  antenna 
system  deciding  on  the  wrong  null  in  t.ie  antenna  array 
pattern.  The  severity  with  which  this  can  affect  the 
calculated  variance  is  graphed  in  figure  43.  This  graph, 
plotted  from  an  equation  developed  in  appendix  B. 
demonstrates  that  a  small  number  of  ambiguities  have  a  very 
significant  effect  on  the  variance.  The  graph  assumes  1020 
bearings  constitute  the  total  sample.  All  the  ambiguities  a~e 
at  133  degrees.  The  abscissa  is  the  number  of  130  degree 
ambiguities.  The  ordinate  is  the  ratio  of  the  new  variance  to 
the  old  variance  which  is  set  at  twenty  degrees  squared.  If. 
for  example,  there  are  ten  193  degree  ambiguities  in  the 
original  sample  and  they  are  removed,  the  new  standard 
deviation  will  be  less  than  ten  degrees  squared.  In  the  case 
of  many  small  sets  representing  short  duration  signals, 
evaluated  separately,  the  reduction  in  variance  will  be 
larger  than  that  predicted  in  figure  43. 

To  examine  the  effect  of  removing  the  ambiguities  from  the 
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SV’RI  data  files,  a  FORTRAN  program  AM3I3-  was  written  to  set 
the  A 4  value  to  13. 6  for  every  record  wnich  is  within  a 
narrow  window  about  the  93  and  1S0  degree  ambiguities. 
Setting  the  A4  to  such  a  large  value  permitted  the  use  of  the 
DFERR  program  which  automatically  rejects  records  witn  such  a 
high  A4  term.  This  test  was  run  on  several  data  files  and  the 
results  are  displayed  in  figures  44  and  45.  It  can  be  seen 
that  sharp  decreases  in  variance  occur,  serving  to  ma.-ce  tne 
bearing  estimate  more  reliable. 

The  resolution  of  ambiguities  with  narrow  aperture  systems 
is  a  difficult  Droller.  The  algorithms  are  especially 
sensitive  to  vertical  angle  of  arrival,  low  SNR  and  multimode 
propagation.  Time  averaging,  if  the  signal  is  sufficiently 
long,  can  overcome  multimode  propagation,  but  vertical  angle 
of  arrival  is  very  dependent  on  array  geometry,  and  SNR  is 
often  totally  uncontrollable.  The  spaced  loop  array  ambiguity 
problem  may  be  due  to  equipment  errors,  but  it  is  more  likely 
that  the  physical  limitation  of  the  aperture  is  the  primary 
difficulty.  The  aperture  is  insufficient  to  sense  wavefront 
distortions  or  to  sharply  isolate  pattern  nulls.  In  tne  case 
of  burst  communications,  time  averaging  is  impossible.  The 
narrow  aperture  array  must  be  assisted  by  wide  aperture 
elements  to  work  against  short  duration  signals.  If  a  spaced 
loop  array  is  mounted  aboard  a  naval  vessel.  it  may  be 
possible  to  add  a  simple  loop  element  at  the  bow  and  stern 
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and  two  amidships.  This  simple  crass  shaped  interferometer 
array  would  he  medium  aperture.  By  itself  it  would  not 
provide  reliable  but  it  could  improve  the  reliability 
of  the  narrow  aperture  system  by  sensing  wavefront 
distortions  and  providing  deep,  well  defined  nulls  to  resolve 
amcigu i ti es . 
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Figure  43 

Ambiguity  Suppression  Curve 
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Figure  44 

Ambiguity  Suppression  for  WWV  5  MHz  2/80 
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Figure  45 

Ambiguity  Suppression  For  WWV  15  MHz  2/80 
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VI.  CONCLUSIONS  ANL  PSCOMMSNDATI ONS 

A.  CONCLUSIONS 

i  •  Performance  and  the  Ionos phere  . 

The  narrow  anerture  antenna  is  physically  incapable  of 
spatiallv  sensing  phase  or  amplitude  wavefront  distortions 
that  have  spatial  periods  of  manv  wavelengths .  This  inherent 
disadvantage  must  be  overcome  by  tine  averaging.  Under  the 
assumption  that  wavefront  distortions  fluctuate  about  tne 
undistorted  wavefront  (the  mean  wavefront),  the  minimum  time 
necessary  to  average  out  the  distortions  is  equal  to  the 

longest  period  of  the  primary  phenomena  causing  the 
distortion.  The  most  severe  distortions  occur  with  multimode 
interference  when  two  o^  more  rays  are  comparable  in 

amplitude.  In  this  case  the  time  required  to  average  out  the 
fluctuations  is  the  fading  period  of  the  major  spectral 
component  of  the  fading  power  spectrum,  roughly  on  the  order 
of  ten  seconds  for  an  HF  signal.  This  means  that  a  narrow 
aperture  antenna  cannot  reliably  determine  the  direction  of 
arrival  of  a  burst  transmission  in  the  presence  of  severe 
multimode  interference.  As  the  severity  of  distortion  is 

reduced,  the  reliability  of  DF  is  improved.  In  the  case  of 

the  SWF. I  spaced  loop  tiFDF  system.  the  range  of  standard 
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deviation  for  the  200  ms  signals  studied  nere  is  from  15 
degrees  for  signals  sysnthesized  from  data  in  file  4  to  59 
degrees  for  signals  synthesized  from  data  in  file  12.  Tne 
average  standard  deviation  over  the  nine  data  files  is  29  1 
degrees.  These  data  files  are  believed  to  include  both 
multimode  and  single  mode  signals  and  signals  with  nigh  and 
low  SN ? 's . 

2.  Performance  Spec if i ca Hons 

The  performance  of  a  narrow  aperture  HFDP  system 
against  short  duration  signals  is  very  dependent  on 
ionospheric  propagation.  Therefore.  the  performance  of  tne 
system  cannot  be  specified  separately  from  specifications  on 
the  state  of  the  ionosphere.  To  state  that  a  s vs  tern  must 
perform  with  a  given  standard  deviation  of  bearing  error 
during  quiet  ionospheric  conditions  with  single  mode 
propagation  is  much  different  from  requiring  the  same 
standard  deviation  of  bearing  error  during  disturbed 
ionospheric  conditions  or  with  multimode  propagation. 

3.  api  bearing  Reliability 

There  are  several  sources  of  fading  in  ionospheric 
propagation?  however,  fading  can  be  considered  a  good  measure 
of  the  distortion  of  local  constant  phase  and  constant 
amplitude  wavefronts.  If  fading  is  severe,  on  the  order  of  12 
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db ,  one  can  be  confident  tnat  significant  distortion  is 
occurring  to  the  wavefront  and  that  narrow  aperture  derived. 
D?  bearings  on  sho^t  duration  signals  will  be  generally 
unreliable . 

4.  Livelihood  Patio  Matrix 

The  concept  of  using  a  likelinood  matrix  developed 
from  given,  reliable  data  to  rake  a  binary  (acceptable  or  not 
acceptable)  decision  on  a  random  data  set  has  not  been  fully 
explored.  However,  the  results  from  the  limited  processing 
are  very  encouraging.  It  should  prove  possible  to  develop 
appropriately  sized  and  numbered  L  matrices  capable  of 
significantly  outperforming  algorithms  based  only  on 
filtering  bv  parameter  limits. 

5 .  Ambigui tv  Elimination 

The  filtering  to  eliminate  Qh  and  130  degree 
ambiguities  demonstrated  that  a  significant  reduction  in 
bearing  variance  can  be  accomplished  if  the  source  of  the 
ambiguities  can  be  corrected.  It  is  doubtful,  however,  that 
the  correction  can  be  accomplished  without  the  introduction 
of  supportive  medium  aperture  interferometer  elements. 


B.  RECOMMENDATIONS 


1.  Use  of  All  Ionospheric  Data 

HFDF  with  narrow  aperture  antennas  on  short  duration 
signals  rust  maximize  the  use  of  as  much  real  time 
information  as  possible  on  the  ionosphere.  Evaluation  of  the 
vertical  angle  of  arrival  and  polarization  is  needed  to  help 
determine  if  a  ground  wave  or  a  sky wave  burst  transmission  is 
being  received.  The  vertical  angle  of  arrival  will  further 
permit  an  estimate  of  range.  The  range  and  vertical  angle  of 
arrival  information  should  be  used  with  a  propagation 
prediction  program  to  make  an  overall  evaluation  of  tne 
reliability  of  a  calculated  bearing  based  on  predicted 
ionospheric  induced  beaming  variance.  To  improve  the 
propagation  prediction  program,  real  time  measurements  of 
ionospheric  parameters  are  needed.  This  could  be  accomplished 
by  updated  inputs  of  geophysical  data,  shipboard  ionospheric 
sounders,  solar  observations,  data  derived  from  satellite 
beacons  and  the  use  of  geographi cal ly  fixed  transmitters  as 
beacons.  The  provision  of  this  information  will  probably  not 
permit  the  calculation  of  a  more  accurate  bearing,  but  it 
will  help  to  determine  the  reliability  of  the  bearing  (the 
variance  of  the  bearing). 
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2.  Use  of  All  Signal  Parameter  Data 

Various  receiver  parameters.  especially  the  AGC 
voltage,  should  be  monitored  to  estimate  the  fading  of  a 
signal.  A  large  variation  of  the  AGC  voltage  could  signify 
fading  and  somewhat  reduce  the  confidence  in  the  calculated 
bearing.  A  small  variation  in  the  AGC  voltage  would  be 
inconc lusive . 

3.  Detailed  Performance  Specifications 


The  specifications  for  the  required  performance  of  a 
narrow  aperture  antenna  system  must  include  the  ionospheric 
conditions  under  which  performance  is  to  be  measured.  An 
example  specification  is  that  a  system  must  have  an  average 
tearing  error  of  five  degrees  and  a  standard  deviation  that 
encompasses  at  least  67  percent  of  the  data  when  tne  signal 
is  received  in  the  20  meter  band  via  a  predominantly  one  nop 
propagation  path.  The  signal  power  of  other  propagation  modes 
should  be  30  db  below  the  primary  mode,  and  tne  SNR  against 
background  noise  should  be  at  least  12  db .  To  conduct  such 
measurements  aboard  a  ship  may  require  that  the  ship  be 
positioned  close  to  a  wide  aperture  array  that  would  be  able 
to  resolve  the  local  mode  structure  of  a  known  target 
transmitter's  signal. 
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4.  Likel inood  Ratio  Matrix  Technique 

Though  the  results  of  using  the  probabalistic 
likelihood  matrix  were  were  not  totally  conclusive,  it  is 
recommended  that  this  approach  ce  further  investigated.  The 
strong  appeal  of  the  likelihood  matrix  approach  is  the 
maximum  use  of  information  in  ar  imprecisely  known  and 
imprecisely  knowaole  environment.  Turtner  investigation 
should  include  the  following  areas.  Optimize  the  dimensions 
of  the  likelihood  matrix,  perhaps  eliminating  tne  phase 
dimension.  Define  the  likelihood  matrix  in  terms  of  tne 
minimum  probability  of  error  criterion  instead  of  the  maximum 
likelihood  criterion.  Define  the  likelihood  matrix  in  terms 
of  Bayesian  costs  and  compare  the  results  with  those  of 
maximum  likelihood  and  minimum  probability  of  error. 
Determine  tne  frequency,  SNR  and  fading  ser.sitivies  of  the 
likelihood  matrix.  Of  greatest  value  and  greatest  difficulty 
would  be  a  general  theoretical  development  of  the  lixelinood 
matrix  ir  te”ms  of  t.ne  spaced  loon  system  narameters.  It  may 
prove  to  be  the  case  that  the  likelinood  matrix  could 
concisely  store  calibration  data.  This  could  be  calibration 
in  azimuth,  elevation,  polarization,  frequency  and  SNR,  ail 
in  a  finite  number  of  likelihood  matrices.  Using  the 
measurements  of  frequency  and  SNR  and  tne  estimates  of 
azimuth,  elevation  and  polarization,  a  stored  likelihood 
matrix  would  postnrocess  system  data  to  refine  tne  oearir.g 
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es  tima  te . 


5.  Ambiguity  Diagnostic  Algorithm 

There  is  a  need  to  develoo  a  diagnostic  algoritnm  that 
will  use  a  known .  fired  transmitter  to  tabulate  ambiguity 
errors  as  a  check:  o'  system  sensitivity  and  cf  the  phase 
sensing  elements,  pa r ticula ri ly  the  nhasemete~. 

€  •  ^§<Liur  Apertupe  Aid 

The  spaced  loop  array  does  not  provide  reliable  data 
for  various  important  tactical  situations.  The  variance 
associated  with  a  burst  signal  bearing  is  too  large  for  most 
fix  and  targeting  algorithms.  To  improve  performance.  the 
feasibility  of  adding  a  simole  medium  anerture  interferometer 
should  be  investigated. 
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APP3NDIX 


AVZBAG3  AND  STANDARD  DEVIATION  CALUL AT ICNS  ON  BEARINGS 

The  c omputer  programs  used  in  the  analysis  portion  of  this 
report  ccrputes  average  and  stania^d  deviation  of  hearings. 
These  statistical  calculations  are  s t raightf o rwa^d ,  cut  not 
in  the  forn  generally  recognized.  In  this  ~eport.  bearings 
have  been  treated  as  integers  in  the  set  d  to  359.  When 
calculating  averages  of  bearings  that  overlap  the  D-359 
boundary,  it  must  be  ta^en  into  account  that  359  differs  from 
0  by  only  one  degree.  For  this  reason  the  following  formulas 
were  developed. 

Assume  that  n  bearings  are  available  for  averaging: 

51,  B  2 . 3  i . Bn 

Normally  the  average  would  be  calculated  by: 

3  =  (1  'n)  £  3i 

rewrite  the  bearings  in  terms  of  a  reference  bearing.  B  ref 
and  a  difference  Ai. 

Bi  =  B  ref  *  Ai 

Then  , 

B  =  (1/n)  £  (3  ref  +  Ai)  =  B  pe?  ■*  (l/n)£  Ai=  3  ref  -A 
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Average  bearing  =  Reference  Bearing  +  Average  Difference 

It  is  most  convenient  to  let  B  ref  =  2‘,  then, 

3  =  A 

However,  the  greatest  difference  allowed  in  a  closed  set  2  to 
35S  is  183  J  therefore,  pax  =  180.  and  the  conditional 

definition  is: 

Ai  =  Bi  if  Bi  <  190 
Ai  =  Bi  -  353  if  Bi  >  130 

Now  the  average  formula  can  be  written: 

B  =  ( 1  -'n)  £  A  i 

subject  to  the  conditional  definition  of  ,  which  is  easily 
ippleper.ted  in  FORTRAN  . 

For  the  standard  deviation,  a  typical  formula  is: 

=  ( 1  /n  )  £  3  i  -  ( 3  )*’ 

Again  let : 

Ai  =  Bi  3  <  Bi  <  183 

Ai  =  Bi  -  353  193  <  Bi  £  353 

Then,  in  a  lixe  panr.er', 

<5^  =  ( 1  /  n )  21  A  i  -  ((lhlE  Ai/ 
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APPENDIX  3 


AMBIGUITY  SUPPRESSION  VARIANCE  EQUATION 

In  the  analysis  portion  of  this  report,  the  results  of 
suppressing  93  and  1SU  degree  ambiguities  from  a  rearing  data 
set  were  discussed.  The  equation  used  to  describe  tne 
significance  of  ambiguity  suppression  is  derived  in  tnis 
appendix  . 

Assume  that  a  histogram  of  bearing  errors  closely 
approximates  the  probabiltiy  density  function  p(b).  where 
p(b)  is  sketched  as: 


The  delta  functions  represent  that  portion  of  tne  histogram 
due  to  ambiguities.  The  function  g(b)  is  general;  the  only 
restriction  associated  'with  it  a"e  that  the  origin  of  the  b 
axis  coincide  with  the  mean  of  gib''  and  that  gib)  be  vide 
compared  with  the  distribution  of  the  ambiguities. 

One  writes  p ( b )  as: 
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p(b)  =  b)  +  ^  ?/?)$(  b-1-! 80  )  +  L  S(  b-*-9? )  +  ?.S(b-90 

+  (?/?)&( b-l£C i  (1 

Using:  J' p ( b  )  db  =  1 
and  f g  (  b )  db  +  P+L  +  P.  =  1 

then  /*(»  db  =  1  -  (P  +  L  +  R)  (?) 

The  variance  of  x»(b'  is: 

VAR f 3]  =y*ban(b)  db  -  |/b 

/b  p(b )  db  =  /b 

/  p(b)  db  =  2  +  ( L )  ( -S0 )  +  ( P. )  (92)  =  (R-LH90) 

VAR  [51  =  / b%  (  b  )  db  +  (P  )  (190  )+(L+R)  (93  >-(?-L  *  (90  *  (3) 

If  the  ambiguities  a^e  eliminated,  a  new  pdf  (b)  results 


b  x>  ( b  )  d  b 


rb  t>(b)  db  =  /b  «(o)  db  *  (P)(18d)2  +  (  L-E  '(9U;2 


g'(b)  is  a  scaled  version  of  g(b)>  therefore. 
K  g'(b)  =  g(b)  where  K  <  1 
/-(b)  db  =  1  ==>  f b!  db  -  K 

combining:  this  with  (?)  yields: 

K  =  1  -  (P  +  L  +  R)  (4) 


Tne  variance  of  the  new  pdf  is: 
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V  A  R  [ 


3']  =  /' b  g'(b)  db  -  j^b  g'(b)  a' 

J  b%'(b>  db  =  (1/X>  ybS:' b)  db 

yb  g'(b)  db  =  0  (by  assumption,  centered  at  the  origin 
VAP.  {?']  =  (1/X)  /b%(b)  db  ( 51 


Substitute  (5)  into  (3)  and  solve  fc"  the  new  variance: 


VAR  [3]  =  X  VAR(3'j  *  (?)(1S0)  +  (L+R)(90)  -  (R-L) (92) 

z  x 

VAR  [B ]  +  ( ( R-L )  -  ( L+R )  -  4?  U 92 ) 

VAP  !  p'l  - -  (6,1 

1  -  ( ?  +  L  +  R ) 

Equation  (o'  is  the  expression  that  approximates  the  new 
variance  of  a  general  distribution  waen  the  ambiguities  are 
eliminated.  Tnere  is  one  condition  that  must  be  observed. 
Because  variance  is  always  a  positive  quantity,  the  numerator 
of  (6^  must  be  a  positive  quantity.  (The  denominator  is 
quaranteed  to  be  positive.)  Therefore, 

2.  "i. 

VARfB]  +  ((R-L)  -  (L+R)  -  4P)(90)  >  0  wnich  can  be 

re-expressed  as, 

^  >  90  I”  4?  +  (  L  +  R  ^  -  (R-L)1’!  ( 7  ) 
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APPENDIX 


n 


PROPHET 

The  data  used  to  compile  the  information  of  Table  I  was 
provided  by  'i1'.  Bob  Pose  of  Naval  Ocean  Systems  Center 
(NOSC).  The  diagrams  included  i  r.  this  appendix  are  sore 
examples  of  the  graphical  output  available  from  tne 
ionospheric  prediction  program  PROPHET. 

Figures  46  and  49  are  ray  trace  plots  of  tne  signal  path 
from  V V  at  Boulder,  Colorado,  to  San  Antonio  Texas.  Tne 
necessary  incuts  to  tne  program  are  the  date,  time  and 
geophysical  data.  The  date"  and  time  were  selected  to 
correspond  with  the  5WRI  soaced  loco  antenna  data.  Tne 
sunspot  number  and  x-ray  flux  were  determined  from  published 
geophysical  data.  The  12.7  cm  flux  is  a  number  that  can  be 
determined  from  the  sunspot  number.  PROPHET  also  requires  the 
transmitter's  location,  power  and  gain  which  in  the  case  of 
'.vVV  was  determined  from  published  sources.  The  traces 
represent  tne  signal  at  launch  angles  from  2  to  50  degrees  in 
5  degree  increments?  the  traces  are  framed  in  a  spatial 
coordinate  system,  altitude  versus  range.  The  receiver's 
location  is  denoted  bv  an  asterisk  at  the  correct  distance 
along  the  range  axis.  This  distance  is  the  great  circle  arc 
connecting  the  transmitter  and  receive". 
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Refraction  causes  the  traces  to  tend  tack  to  the  Earth 
unless  the  launch  angle  is  sufficiently  high  to  permit  tne 
ray  to  escape  the  Earth.  This  is  the  case  for  the  50  degree 
launch  ray  in  figure  46.  Multimode  occurs  when  ray  traces 
cross  each  other  and  return  to  the  surface  at  approximately 
the  same  position.  In  figure  49  there  is  a  potential 
multimode  condition  at  the  receiver.  The  severity  of 
multimode  interference  depends  on  the  strengths  of  tne 
multihop  traces  after  reflection  at  the  surface  of  the  Earth 
and  D  layer  adsorption. 

Figures  4?  and  50  are  relative  power  diagrams  in  a 
frequency  versus  24  hour  coordinate  system.  The  curves  are 
relative  rover  contours.  The  tor  cuve  represents  tne  MUF,  and 
its  relative  power  value  is  -30  do.  The  bottom  contour  is  tne 
LUF ,  and  it  too  has  a  relative  rower  value  of  -30  db.  As  one 
works  inward  from  the  outer  two  contours,  each  succeeding 
contour  is  a  +1?  db  higher  than  its  outside  neighbor.  Tne 
inside  contours  may  reach  values  of  +10  and  +2 l  db.  Tne 
actual  power  received  is  dependent  upon  the  transmitter  power 
and  the  propagation  conditions. 

Plots  of  ionospheric  induced  variance  are  shown  ir.  figures 
4£  and  51.  Tne  values  of  variance  are  especially  useful  for 
direction  finding  work  in  that  they  are  a  measure  of  error 
induced  by  the  ionosphere  independent  of  any  HFEF  system.  Tne 
values  a +e  calculate'1  from  empirically  derived  formulas.  The 


13? 


fi^st  order  effect  is  frequency  dependent,  and  toe  second 
order  effect  is  cased  on  the  Poss  curve  of  variance  as  a 
function  of  range.  Reference  13  is  the  M OS G  documentation  on 
the  development  of  the  emT>i"ical  formulas.  Figures  46  and  51 
show  that  variance  is  typically  oetveen  1  and  2  degrees 
squared  with  occassional  sharp  peaks  reaching  3  degrees 
squared  for  about  one  hour  duration.  The  major  peatc  at  1333Z 
is  due  to  sunrise  effects  (terminate’'.;. 

Additional  information  on  the  capabilities  of  PR0P.-1ET  car. 
be  obtained  from  M03C . 
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5  hhO  1900  15002 

SUNSPOT  HUMBER  *  181 
10.7  CM  FLUX  =  225 
X-RAY  FLUX  =  1.000E-004 


LAUNCH  ANGLES: 
MIN  -  0.0  DEG 

MAX  »  50.0  DEG 

INC  *  5.0  DEG 


TPra/iMiM^  m 


GO 

CD 


transmitter:  boulder 

40. 8N,  105. IN 
POWER  =  10000  WATTS 
ANTENNA  GAIN  =  0  dB 
FREQUENCY  =  15000.0  kHz 


receiuer:  sa 

29. 5N,  98. 6W 

RANGE  »  1387  KM 

BEARING  =  152.9  DEGREES 


> 


Figure  46 

WWV  15  MHz  2/5/80  Ray  Trace 


SUNSPOT  NUMBER  *  181 
10.7  CM  FLUX  ®  225 
X-RAY  FLUX  1.000E-004 


tfihi/irphk  m 


•P 

O 


TRANSMITTER:  BOULDER 
40. 8N,  105. 1W 
POWER  =  10000  WATTS 
ANTENNA  GAIN  =  0  dB 
FREQUENCY  =  10000.0  kHz 


RECEIUER: 
29. 5N, 
RANGE  * 
BEARING 


SA 

98.614 
1387  KM 

«  336.6  DEGREES 


MUF 


FREQ 


LUF 


:  > 


Figure  47 

WWV  10  MHz  2/5/80  Relative  Power  Diagram 
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5  FEB  1980  BQHBB  FRQ  »  15008  KHZ 
X-RAY  FLUX  =  1 . 00E-004  SUNSPOT  *  «  181 

tgt:  target  lat  «  40.80;  lon  «  105.10 

PWR  =  10000  WATTS;  ANT  GAIN  =  0  DB 


u 

A 

R 

I 

A 

N 

C 

E 


TIME,  ZULU 

Figure  48 

WWV  15  MHz  2/5/80  24  Hour  Variance  Diagram 


13  FEB  1979  18302 

SUNSPOT  NUMBER  *  151 
10.7  CM  FLUX  *  194 
X-RAY  FLUX  =  1 . 000E-004 


LAUNCH  ANGLES: 
MIN  »  0.0  DEG 

MAX  «  50.0  DEG 

INC  «  5.0  DEG 


TPIWIFSas  m 


M 

-P 

fNj 


TRANSMITTER:  BOULDER 
40. 8N,  105. 1W 
POWER  =  100O0  WATTS 
ANTENNA  GAIN  =  0  dB 
FREQUENCY  =  15000.0  kHz 


receiver:  sa 

29. 5N,  98. 6W 

RANGE  *  1387  KM 

BEARING  *  152.9  DEGREES 


:  > 


Figure  49 

WWV  15  MHz  2/13/79  Ray  Trace 


13  Ft B  1979  BBBII 
SUN-1  POT  HUMBER  =  151 
10.7  CM  FLUX,  =  194 
X-RAY  FLUX  =  1.080E-004 


-P 

CO 


TRANSMITTER:  boulder 
40. SN,  105. 1W 
POWER  =  10000  WATTS 
ANTENNA  GAIN  =  0  dB 
FREQUENCY  =  10000.0  kHz 


receiuer:  sa 

29. 5N,  98. 6W 

RANGE  *  1387  KM 

BEARING  =  336.6  DEGREES 


F 

R 

E 

Q 

U 

E 

N 

C 

Y 


MUF 


FREQ 


LUF 


: >00000000000 


Figure  50 

WWV  10  MHz  2/13/79  Relative  Power  Diagram 
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13  FEB  1979  BBBSB  FRQ  *  10000  KHZ 
X-RAY  FLUX  =  1 . 00E-004  SUNSPOT  #  *  151 

tgt:  target  lat  *  40.80;  lon  =  105.10 

PWR  =  10000  watts;  ant  gain  =  0  DB 


[fcgp’/lfcSHK  g] 


-p 

-p 
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Figure  51 

WWV  10  MHz  2/13/79  24  Hour  Variance  Diagram 


non  ooo  ooo 


DFERR 


-p 

cn 


C  THIS  IS  A  PROGRAM  TO  CALCULATE  OF  ERROR  STATISTICS  OF  THE  SWRI 
C  SPACED  LOOP  AMT  EMMA  SYSTEM.  A  5IVEN  SET  OF  BEARINGS  IS  ANALYZED  TO 
C  FIND  THE  MEAM  AND  STANDARD  DEVIATION.  THESE  STATISTICS  ARE  USED  TO 
C  RE-EVALUATE  THE  DATA  SET  BY  TRIMMING  OFF  BEARINGS  WITH  LARGE 
C  DEVIATIONS  FROM  THE  MEAN.  A  NEW  MEAN  IS  C  AL  CUL  AT  ED  AND  COMPARED 
C  TO  THE  TRJE  BEARING  TO  p  I  NO  THE  FR.ROR.  STATISTICS  ARE  THEN 
C  COMPUTFD  FOR  GIVEN  SIGNAL  DJRATION  ACROSS  THE  ENTIRE  DATA  SET 
C  AMD  FOR  SE.ECTED  VALUES  OF  A4  WHICH  IS  A  MEASURE  OF  SYSTEM  NOISE. 

C 

C 

LQGICA.n  SOURCE!  30J /30*'  '/ 

I NTEGER  *2  N9RG2 

data  :lr/,lr,',:sr/,sr'/,csa/,sa*/ 

DIMENSION  NBVEC!  5  000) , A4VEC!  5  000) ,ND AT  A ! 1 0 , 10 , 4 ) 

C  FOR  CP/CMS,  MJST  REFORMAT  DSRN  01 

CALL  DSDSET (I, 1500,2,16) 

WRITE!  I  0,  600) 

WRITE! 10,601) 

READ! 5, 500)  SOJRCE 
WRITE! 10,602) 

READ! 5, 501 )  NREC 
WRITE! 10, 603) 

RE  AD ! 5 , 502)  MINS 
1M I N S  =  M I  NS720 
WRITE! 10,604) 

R  EAD! 5 , 501 )  MAYS 
M AXS=MAXS/20 
WRITE! 1 0, 605) 

READ!  5 , 50  2)  NINE 
NINC=NI NC/20 
WRITE! 10,606) 

READ ( 5 , 503 )  NT3RG 
WRITE! 10, 607) 

R  EAD ( 5 , 50  4 )  NS  ID 

CALCULATE  ROTATION  SHIFT  TO  SIMPLIFY  DETERMINATION  OF  BRG  ERROR 
NSHIFT=360-NT}RG 

SET  UP  LOOP  FOR  A4MAX  INCREMENTING 

DO  100  11=1,10 
A4MAX=3.1*FL0AT! 11-11) 

SET  UP  LOOP  FOR  SIGNAL  DJRATION  DECLARATION 


BRG0621 0 
BRG06220 
3RG06230 
BRG06240 
3RG06250 
BRG06260 
BRG06270 
BRG06280 
BRG06290 
BRG06300 
BR306310 
BRG06320 
BRG06330 
R  RG06340 
BRG06350 
BRG06360 
BRG06370 
BRG06380 
BRG06390 
BRG06400 
BRG0641 0 
RRG06420 
BRG06430 
9RG06440 
BRG0645  0 
BRG06460 
RRG0647  0 
BRG06480 
BRG06490 
BRG06500 
BRG06510 
BRG06520 
BRG0653  0 
BRG06540 
BRG06550 
BRG06560 
BRG06570 
BRG06580 
BRG06590 
BRG06600 
BRG066 1 0 
BRG06620 
RRG06630 
BRG06640 
BRG06650 
BRG06660 
BRG06670 
BRG06680 


noon  OOO  OOO  OOOOr—  o  oo  ooo  ooo 


DO  10L 
NC  1  =  0 
NC2  =  0 
M2  =  0 
M  3=0 
M4=0 
M  10  =  0 
M  12  =  0 


ji=mins  ,maxs,nin; 


UP 

LOOP  TO 

READ 

THRU 

ENTI  RE 

F 

IL 

E 

M  5= 

NR  E 

C/Jl 

00 

102 

K1  =  1 

» M5 

UP 

LOO 

?  TO 

READ 

NUMB 

E 

R  OF  REC 

S 

CORRES 

PON 

DIN 

G 

TO 

DO 

103 

Ll  =  l 

t  J  L 

READ!  1 

,  505 1 

END  = 

20  1 » 

E 

R 

R=103 ) 

A 

4, 

SStNBR 

32 

N8R 

G  =  N 

BRG2 

S  = ' 

, • SR*  <  OR 

’  S  A 

•  IF* 

» 

THE  RE 

FE 

RE 

NCE 

SI 

GMAL  I 

S 

LOW 

LO 

BE 

0ISCA 

ROE  0 

.  THI 

S 

IS  00 

NE 

BY  SE 

TT 

!  M3 

THE 

A  4 

IF  ( 

SS. 

EQ.CLR.0R 

.SS. 

E 

Q 

.CSR.O 

R. 

SS 

.EQ. 

CS 

A) 

A4  = 

10 

.0 

COUNT 

fOTAL 

NJM 

BER 

0 

F 

RECOR 

OS 

R 

EAD 

HI 

TH 

NC2 

03 


NC2=NC2+1 
A4VEC { L 1 ) = A4 
M  B  V  E  C  (  L 


(  L 

1 J  =N 

BR  3 

ONTA 

INS 

BRGS  FOR  3 

I  YEN 

A 

SSOC 

14  TE 

0  WITH  THE 

SE  3 

TR 

I M  (  M 

bve: 

♦ A4VE C  » J 1 » 

NS  TO 

IC 

I  ENT 

DAT 

A  TO  DETERMINE 

SJ 

F.EQ 

.1  > 

GO  TO 

200 

p 

RROR 

C/> 

a? 

in 

TWEEN 

TRUE 

AND 

4N+-NSH 

IF  T 

AM 

.GT. 

180) 

ME4N  = 

=  ME  AN 

-360 

I F ( Ml . 3  E. 36  0)  Ml =Ml-360 
I F( Ml «GE. 180)  Ml=Ml-360 

TOTAL  THE  ERRORS  IM  M2  AMO  ERRORS  SQUARED  IN  M3 
TOTAL  MEAN  IM  M12  AND  MEAN  S3J4RED  IN  M10 


M2  =  M2*-M  1 
M  3=M3  -frM 1 *M1 


BRG06690 
BRG06700 
BRG0671 0 
BRG06720 
BRG06730 
BRG06740 
BRG06750 
BRG06760 
BRG06770 
BRG06780 
BRG06790 
BRG06800 
BR  GO 68 10 
BRG06820 
BRG06830 
BRG06840 
BRG06850 
BR306860 
BRG06870 
BRG06880 
BRG06890 
BRG06900 
BRG0691 0 
BRG06920 
BRG06930 
3RG06940 
BRG06950 
RRG06960 
BRG06970 
BRG06980 
8RG06990 
B  RG07000 
BRG07010 
BRG07020 
BRG07030 
BRG07040 
3RG07050 
RRG07060 
BRG07070 
BRG07080 
BRG07090 
BRG071 00 
BRG07110 
BRG07120 
BRG07130 
BRG07140 
BRG07150 
BRG07160 


M10=M1 D+MSTD 


-p 


C 

C  COUNT  NUMBER  OF  SIS 
C 

nci=n:h-i 

GO  TO  L 02 
200  IM4=M4  +  I  N5UF 

102  CONTI N  JE 

C 

C  READ  TO  EM D  OF  DATA 
C 

if<nc2.ge.nre: ) 

M8=NREC-NC2 
DO  104  K2*lfM3 
104  READC 1, 505, ENO= 
NBRG=NBRG2 
206  REWIND  1 
C 

C  FIND  AVE  BEARING  ER 
C 

IF(NCl.EQ.O)  GO 
GO  TO  203 
C 

C  IF  NO  SIGNALS  FOR  W 
C  THE  OUTPUT  WI T H  999 
C 

202  M6=99999 
M  7=Mfe 
NC1=M6 
M4  =  M6 

M 1 3=  M  6 
M  14  =  M6 
GO  TO  204 
C 

C  CALCULATE  MEAN  BEAR 
C  CALCULATE  MEAN  BRG 
C 

203  M6=M2/NCi 
X1=FL0AT( M2 ) / F  _ 
X2=FL0AT(M3)/FL 
M7= I NTC SORT ( X2- 
M 14=M  1 O/NC  1 

Mll  =  Jl /MI  NS 
C 

C  STORE  VALJE  TO  BE  3 
C 

204  NDATAC I  1, Ml  1, L ) 
NDATAC  I  1 1  Ml  1 «  2 ) 


MALS  FOR  WHICH  STATS  ARE  OBTAINED 


FILE  IF  NECESSARY  INORDER  TO  REWINO  FOR  NEXT  RUN 
GO  TO  206 

2  0  61  A4»  SS  »  N3RG2 

R  AND  ST J  OF  BEARING  ERR 
TO  202 

HICH  STATS  OBTAINED,  THIS  SHOULD  BE  FLAGGED  IN 
99. 


ING  ERROR  ( M 6 )  AND  STD  OF  BEARING  ERROR  (M7J. 
( M 1 4 )  AND  STD  OF  BRGS  ABOUT  MEAN  C  M 13 ) 


OAT(NCl) 
OAT(NCl) 
X1*X1) ) 


RINTED  IN  THE  ARRAY  *  NDATA '  . 

=  M6 
=  M  7 


BRG07170 
BRG07180 
BRG07190 
BRG07200 
BRG0721 0 
BRG07220 
BRG07230 
BRG07240 
BRG07250 
BRG07260 
BRG07270 
BRG07280 
BRG07290 
BRG07300 
BRG0731D 
BRG0732U 
BRG07330 
B  RG07340 
BRGD7350 
BRG07360 
BRG07370 
BRG07360 
BRG07390 
RRG07400 
9RG074I 0 
RRG07420 
BRG07430 
BRG07440 
BRG07450 
BRG07460 
BRG07470 
BRG07480 
3 R GO  749  0 
BRG07500 
BRG07510 
PRG07520 
BRG07530 
BRG07540 
BRG07550 
BRG07560 
BRG07570 
BPG07580 
BRG07590 
BRG07600 
BRG076 1 0 
BRG07620 
BRG07630 
BRG07640 


101 

100 

201 

205 

105 

106 

C 

C 

500 

501 

502 

503 

504 

505 

C 

C 

600 

601 

602 

603 

604 

605 

606 

607 

608 

609 

610 

611 

612 

613 

614 


GO  TO 


DO 


(  1 1 , Ml  1, 3 ) =NC1 

PRG07650 

(  I  1 ,  Ml  1 , 4 ) =  M14 

BRG07660 

NUE 

BRG07670 

NJE 

BRG07680 

205 

BRG07690 

(6,609) 

SRG07700 

NJE 

BRG0771 0 

6  11=1 ,3 

BRG0772D 

(6,614)  SDJRCE, NSTD 

BRG07730 

5  12=1,4 

BRG07740 

. E  Q.  1 )  WRI  TE  (6, 6  10) 

BRG07750 

. EQ. 2 )  WRITE(6»611) 

RRG07760 

.EQ.3)  WRI  T  E  (  6 , 6 12  ) 

RRG07770 

. E  0. 4 )  WRI TE ( 6, 6  13) 

BRG07780 

PRINT! NDATA, 12, MINS, N I NC  ) 

BRG07790 

DO 


CALL 
CONTI M  J  E 
STOP 


FORMAT ( 30 A 1 ) 

FORMAT! 15 ) 

FORMAT! 14) 

FORMAT! 13 ) 

FORMAT!  ID 
FORMAT! 8X,44,42,  A2) 


FORMAT! 
A '  0  OF  S 
FORMAT! 
FORMAT! 
FORMAT! 
FORMAT! 
FORMAT! 
FORMAT! 
FORMAT! 
FORMAT! 
FORMAT! 
FORMAT! 
A*  NOISE 
FORMAT! 
A ' OF  STS 
FORMAT! 
A* AS  A  F 
FORMAT! 
FORMAT! 
A*  USED  T 
END 


DTHIS  P 
GNAL  DJ 


OENTER 
•OENTER 
•OENTER 
' CENTER 
• OENTER 

•  OEM  TER 
• OENTER 
1017  ) 

• OEMD  0= 

•  OAVFRAG 
( A4  TERM 
• OSTAMDA 
TEM  NO  I S 
• ONUMBER 
UNCTION 

'  0  AV  E  ST 
' 1S0URIE 
0  DETERM 


FINDS  AV E  BEARING 
AND  SYSTEM  NOISE 
AND  T IME! 30A1) •  ) 

OF  RECORDS!  15) • ) 

_  _ N AL  JURAT  ION, XPLE 

MAX  SIGNAL  DURATION, XPLE 
SIGNAL  INCREMENT, XPLE  OF 
TRUE  BEARI NG  !  13  )  •  1 
STD  XPLE  FOR  TRIMMING  (ID* 


ROGRAM 
RATION 
SOURCE 
NUMBER 
MIN  SI 


ERROR  AS 
!A4  TERM) 


FUNCTION'  ,  /, 


OF 

OF 

20 


20  MSEC!  14) ' ) 
20  MSEC!  15)'  ) 
MSEC! 14)  '  ) 


) 


DATA  FILE  ENCOUNTERED.') 

E  REARING  ERROR  AS  A  PJNCTION  OF  SYSTEM 
)  AND  SIGNAL  DURATION',//) 

RD  DEVIATION  OF  BEARING  ERROR 
E  AND  SIGNAL  DJRATION* , //) 

OF  VALID  SIGNALS  OF  A  GIVEN  DURATION  ', 

OF  SYSTEM  NOISE  AND  SIGNAL  DURATION',//) 
ANDARD  DEVIATION  OF  BRGS  ABOUT  AVE  BRG* , // ) 


AS  A  FUNCTION  • 


:  ',30A1,/,' 

I NE  BEARING 


OSTANDARD 
WINDOW  =  • 


DEVIATION 

11,//) 


MULITPLE  « 


BRG07800 
BRG0781P 
BRG07820 
BRG07830 
BRG07840 
BRG07850 
BRG0786O 
BRG07870 
BRG0788D 
BRG07890 
BRG07900 
BRG0791 0 
BRG07920 
8RG07930 
BRG07940 
BRG0795  0 
BRG07960 
BRG07970 
BRG07980 
3RG079VD 
BRG08000 
BRG08010 
BRG08020 
BRG08030 
BRG0804D 
BRG08050 
BRG08060 
BRG08070 
BRG08080 
BRG08090 
BRG08100 
BRG08110 
BRG08120 


ooo  ooooo  oooo  ooo  ooooooooooooo 


SUBROUTINE  TRIM 

THIS  SUBROUTINE  CALCULATES  MEAN  AND  STD  OF  •  N  •  ELEMENTS 
(0-360  BEARINGS)  0  =  VECTOR  NVEC  AFTER  FILTERING  OUT  NVEC  TERMS 
ASSOCIATED  WITH  A*  VALUES  IN  XVEC  THAT  ARE  TOO  LARGE.  FILTERED 
ELEMENTS  ARE  STORED  IN  IVEC.  USING  MEAN  AND  AN  INTEGER  MULTIPLE 
OF  STD,  ELEMENTS  0=  IVEC  ARE  FILTERED.  VA.JES  WITHIN  LIMITS  OF 
MEAN  ARE  USED  TO  COMPUTE  A  NEW  MEAN  (MEAN).  THERE  MUST  BE  AT  LEAST 
VALUES  IN  NVEC  AND  AT  LEAST  2  VALUES  AFTER  FILTERING,  IF  NOT,  THE 
VALUE  OF  ' I NSUF '  IS  SET  TO  1 


SUBROUTINE  TR  [  M(  NVi 


XV  E! 


DIMENSION  NVEC<  N)  ,XVEC(  N) 
DATA  IVEC/500D*0/ 


, N , NSTD, XMAX , M S T 0 , ME AN , I N SU F ) 
, I VECI 5000) 


INITIALIZE  COUNTERS  AND  SUMMERS 

N  C I  =  0 
M  1  =  0 
M2  =  0 
I NSUF=D 
ME  AN=  0 
M  =0 
M  11  =  0 


IF  LESS  THAN  3 
I  NS JF  =1 . 


BEARINGS  IN  NVE 


t 


RETURN  TO  MAIN  PROGRAM  WITH 


IF(N.LE.Z)  GO  TO  200 


FILL  IVEC  WITH  ELEMENTS  OF  NVEC  THAT  PASS  A4MAX  FILTERING 

COUNT  NUMBER  OF  VALID  NVEC  ELEMENTS  WITH  NCI 

Ml  AND  M2  SUM  TERMS  FOR  FINDING  MEAN  AND  STANDARD  DEVIATION 


100 

C 


BRG08130 
8RG08140 
BRG08150 
BRG08160 
3RG08170 
BRG08180 
BRG08190 
BRG08200 
BRG08210 
BRG08220 
BRG08230 
RRG08240 
BRG08250 
BRG08260 
BRG08270 
BRG08280 
BRG08290 
BRG08300 
BRG08310 
BRG08320 
BRG08330 
BRG08340 
BRG08350 
BRG0836D 
BRG0837  0 
BRG08380 
BRG0839  0 
BRG08400 
BRG0841 0 
BRG0842  0 
BRG08430 
BRG08440 
BRG08450 
BRG08460 
BRG0847  0 
BRG08480 


DO  100  11=1 ,N 

3RG08490 

IFIXVECdl)  .GT.XMAX)  GO  TO  100 

BRG08500 

BRG08510 

OUNT  NUMBER  OF  BR3S  THAT  PASS  XMAX  FILTERING  WITH  NCI 

BRG08520 

BRG0853  0 

NC1=NC 1 +1 

BRG0854D 

IF(NVEC  (  I  1)  .3T.  180)  NVEC(  I  1)  =  NVEC(  ID-360 

BRG08550 

IVEC(NCl) =NVEC(  I 1  ) 

BRG08560 

M1  =  NVEC  (  ID  +M1 

BRG08570 

M2  =  NVEC(I1)*NVEC(  Il)+M2 

BRG08580 

CONTI NUE 

BRG08  59  D 

BRG08600 

C  IF  LESS  THAN  2  BRSS  IN  IVEC,  STATS  WILL  bE  MEANINGLESS,  RETURN  TO 
C  MAIN  PROGRAM  WITH  INSUF=1. 

C 

I F<  NCI . LE. 1  )  j 0  TO  200 
C 

C  CALCULATE  MEAN  (M5)  AND  STD  (M3)  OF  XMAX  FILTERED  ELEMENTS. 

C  FLOATING  POINT  IS  JSED  IN  SOME  OPERATIONS  FOR  BETTER  ACCURACY. 

C 

X1=FL0AT(M1 )/FLOAT(NCI) 

X2=FL0AT(M2)/:_DAT(MC1) 

M5  =  I NT( XI  ) 

IF(M5.LT.O)  M5=M 5+360 
M  3  =  I  N  T  ( SORT ( X2 -  X 1 *X  1 )  ) 

IFIM3.LT.6)  M3  =  6 
C 

C  APPLY  STD  MULTIPLIER  NST3,  CALCULATE  ACCE3TABLF  WINDOW  ( M6CM7 ) 

C  ABOUT  THE  MEAN  (M>)  AND  ZEROIZE  COJNTER  M  AND  SUMMER  M9. 

C 

M4=NSTD*M3 
M6=M5+M4 
M  7  =  M  5  -  M  4 
N  W I  N  =  M  5  — M  7 
M  =  0 
M9=0 
M  10  =  0 


C 

C  FIND  BEARINGS  WITHI 
C 

DO  101  12  =  1, NCI 
M8=IVEC (12) +3SD 
I F( ( I  YE  C (  I2J.LT 
M=M  +  1 

I F( I VEC ( I  2 )  .LT. 
M9  =  M9  +  l VEC(  12) 
M10=M10+I VEC( 12 
101  CONTI N J  E 

C 

C  COMPUTE  MEAN  OF  BE  A 
C 

IF(M.EO.O)  GO  T 
X4=FL0AT( M9 )/FL 
X5=FL0AT(  M10)  f- 
MST0=INT(SQRT( X 
ME AN= I  NT ( X4 ) 

GO  TO  201 
C 

C  IF  NUMBER  OF  VALID 
C 


N  WINDOW  ABOUT  THE  MEAN 

. M7.0R.IVECI 12)  .GT  •  M6 ) .AND.M8.GT .M6) 
M7)  I  V  EC ( 1 2 ) = I  V  E  C ( I  2 )  +  3  S  0 
)*IVEC(  12) 

RI  NGS 

0  200 
OATI  M ) 

L  0  AT (  M ) 

5-X4*X4)J 

ELEMENTS  NOT  SUFFICIENT,  SET  INSUF=1 


GO 


T0  101 


BRG08610 
BRG08620 
BRG08630 
BRG08640 
BRG03650 
BRG08660 
3  R  GO  86  70 
BRG08680 
BRG08690 
BRG08700 
BRG06710 
BRG08720 
B  RG  0873  0 
BRG08749 
BRG08750 
BRG08760 
BRG08770 
B RG0878  0 
Bp  G08790 
BRG08800 
BRG088 10 
BRG08820 
BRG08830 
BRG08840 
BRG08850 
BRG08860 
BPG08370 
BRG08380 
BRG08890 
BRG0890  0 
3RG08910 
BRG08920 
BRG08930 
3RG08940 
BRG06950 
BRG08960 
RRG0897  0 
BRG08980 
3RG08990 
BRG09000 
BRG09010 
BRG09020 
BRG09030 
BRG09040 
BRG090S0 
BRG09060 
BRG09070 
BRG09080 


200  I  N  S  U  F  =  1 

201  RETURN 
END 

C 

C 

C  THIS  SUBROUTINE  IS  USED  TG  PR I M  T  THE  DAT  A  CALCULATED  IN  THE 
C  MAIN  PROGRAM 
C 

SUBROUTINE  PR I NT(NDATA, N, MINS, NINCJ 
L0GICAL*1  TITLE! 10) 

DATA  TITLE/*  •,*  •  , * M • ? • A • ,  1 X  •  »  •  «,'A*,'4*,‘  * , •  •/ 

DIMENSION  A4LB.  (  10) ,NDATA(  10,  10,4) ,NSD-BL<10) 

DO  100  11=1,10 
A4LBL  (I  1)  =FLDAH  11-1  1  )*0.1 
100  NSDLBL!  II 1=MINS*20+! ( 1 1- 1 i *NI  NC*20) 

WRITEtb ,600)  ( TITLE! I  ) , A4L3  H  I I , ( NDATA! T , J , N ) , J  =  1 , 10 ) , I = 1 , 10) 

600  FORMAT! 5X, 41,5 X,F 3.1, 2X, 101  6, /I 
WRITE!S,601)  NSDLBL 

601  FORMAT! //,16X, 1016) 

WRITE! 5,602 ) 

602  FORMAT!//, 25X,» SIGNAL  DJRATION  t MILL  I  SEC ) ' / // ) 

RETURN 

END 


BRG09090 

BRG091 00 

BRG09110 

BRG09120 

BRG09130 

BR  G091 40' 

BRG09150 

BRG09160 

3RG09170 

BRG09180 

BRG09190 

BPG09200 

BR GO 92 1 0 

BRG09220 

BRG09230 

BRG09240 

BRG09250 

BRG09260 

BRG0927  0 

BRG092B0 

BRG09290 

BRG09300 

BRG09310 


non  oooo  ooo 


LMAT 

C  THIS  PROGRAM  READS  4  FILE  OF  SWRI  SPACED-L 
C  A  THREE  DIMENSIONAL  C  QND  I T  I OM  AL  PROBABI  LTY 
C  MATRIX  (U 
C 

LOGICAL*!  SOURCE!  30) /30**  •/ 

INTEGER*?  NB2 

DIME NS I  ON  A(1D,10,10),U!10,10,10),XL(1 

SET  UP  FOR  Co/CMS  READ  OF  DATA  FILE 

CALL  DSDSET(1 ,1600,2,16) 

INITIALISE  THE  ACCEPTABLE  (A),  THE  UNACCEP 
LIKELIHOOD  RATIO  (XL)  MATRICES  WITH  ALL  ZE 

DATA  A/ 10  00*0.  / ,U/1000*0.  /,  XL/ 1000*0./ 

SET  IN  PARAMETERS: 

DATA  CLR/'LR'/, CSR/*  SR* /,  CSA/'  SA«/ 

WRI  TE  ( 1  0,601) 

601  FORMAT! 'DENTER  SOURCE  (3DA1)') 

READ ( 5, 501 )  SDJRCE 

501  FORMAT ( 30 A1 ) 

WRITE!3,602)  SOURCE 

602  FORMAT! 30A1) 

C  NUMBER  OF  RECORDS  AND  NJMBER  0=  ^ILES 
NREC= 1 0000 
NFI  LE=4 

C  SET  UP  TWO  ACCEPTABLE  BEARING  WINDOWS  TO  A 
C  CROSSOVER  INHERENT  IN  BEARING  PROBLEMS.  E. 
C  WINDOW  AROUND  B=350,  SET  UP  TWO  WINDOWS  33 
C  THE  WINDOW  DOES  NOT  SPAN  THE  CROSSOVER,  SE 
C 

C  LOWER  BOUND  OF  ACCEPTABLE  BEARING  WINDOWS 
NBMI N 1 =0D0 
NBMI N2=  292 

C  UPPER  BOUND  DF  ACCEPTABLE  BEARING  WINDOWS 
NRMAX1  =  02  2 
N  BM  AX2  =  36  0 

C  INITIALIZE  ACCEPTABLE  AND  UNACCEPTABLE  BEA 
N  A  =  0 
NU=0 

C  INITIALIZE  CONSTANTS  NEEDED  FOR  THE  SUBROJ 
C  AO,  PHASE,  AND  A4  INTO  APPROPRIATE  MATRIX 
AOMI N  =  D  . 


OOP  DATA  AND  CONSTRUCTS 
LIKELIHOOD  RATIO 

D, 10,  10) 


TABLE  (U)  AND  THE 
ROS. 


CCOUNT  FOR  360  TO  000 
G.  IF  YOU  WANT  A  20  DEG 
0-360  AND  000-010.  IF 
T  WINDOW  2  =  WIMOOW  1. 

I  AND  2 


1  AND  2 


RING  COUNTERS 


TINE  THAT  MAPS 
ADDRESS  . 


BRG03990 
BRG04000 
BRG04010 
BRG04020 
BRG04030 
BRG04040 
BRG04050 
BRG04060 
BRG04070 
BRG04080 
BRG04090 
BRG04100 
BRG04110 
BRG04120 
B  R  GO  4130 
9RGD4140 
BRG041 50 
RRG04160 
BRG04170 
BRG04180 
BRG04190 
BRG04200 
BR GO 42 10 
3RG04220 
BRS04230 
BRG04240 
BRG04250 
BRG04260 
BRG04270 
BRG04280 
BRG04290 
B  RG04300 
BRG04B10 
BRG04320 
BRG04330 
PRG04340 
BRG04350 
BRG04360 
BRG04370 
BRG04380 
BPG04390 
BRG04400 
BRG04410 
BRG04420 
BR GO 4430 
BRG04440 
BRG04450 
RRG04460 


A  OMAX  =  L  . 

NAODI ^=13 
PMI N=-l 80. 

PMAX=1 80. 

NPDIM=1 0 
A  4MI N  =  0  . 

A4MAX=I . 

N  A4DI M= 1 0 

USE  A  DO  LOOP  TO  RE  AD  THRU  THc  ENTIRE  FILE  AND  PLACE 
APPROPRIATE  ENTRIES  INTO  THE  A  AMO  U  MATRICES. 

C 

DO  103  18=1 tNFILE 
DO  109  11=1,  NREC 

READ(  1,500,  E  N  0  =  109)  AO,  NP  ,  A4,  SS  ,NB2 
NB=NB2 

500  FORMAT! 3A4, A2,A2> 

IF(SS.EQ.CLR)  GO  TO  109 
IF( SS.EQ.CSR)  GO  TO  109 
IF(SS.EQ.CSA)  GO  TO  109 
P  =  FLOAT ( N  PI 


C 

C 

r 


c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


CALL  SUBROUTINE  M A 3  F  TO  DETERMINE  MATRIX  ELEMENT  ADDRESS  (I,J,K). 
AO  DETERMINES  I,  3HASE  DETERMINES  J,  AND  A 4  DETERMINES  X. 


DECIDE  IF  A 


C 

C  Tt 
C 

200 


NU=NU+1 


N  A=NA  + 1 


PF 

(AO, AO 

MIN,  A 

PF 

<  P  .PMI 

N,  PMA 

PF 

<  A4,  A 4 

MIN.  A 

1 

SH3JL 

D  BE 

e 

NBMI  Nl 

.AND. 

E. 

NBMI N2 

.  AND. 

AU 

GMENT 

UN  ACC 

)  = 

U<  I, J, 

0+1. 

09 

AUGMENT 

ACCEP 

)  = 

A(  I,  J, 

<)  +1. 

109 
108 
C 

C  CONSTRUCT  LIKELIHOOD  RATIO  MATRIX 


CONTI NJC 
CONTI NUE 


BRG04470 
RRG04480 
BRG0449  0 
BRG04500 
BRG045 1 0 
BRG04520 
PRG04530 
BRG04540 
BRG04550 
RRG04560 
BRG04570 
BRG04580 
BRG04590 
BR  GO  4600 
n  RG0461 0 
BRG04620 
n  R  G0463  0 
BRG04640 
BRG04650 
BRG04660 
BRG04670 
BRG04680 
BRG04690 
BRG04700 
9RG04710 
BRG04720 


,  J  ) 

BR GO 4730 

40  I  M »  K  1 

BRG04740 

BRG0475  0 

0  A  OR  U. 

BRG04760 

BRG04770 

BMAXl)  GO  TO 

200 

BRG04780 

BMAX2 )  GO  TO 

200 

BRG04790 

BRG04000 

COUNTER  AND 

MATRIX. 

BR  G  04810 

BRG04820 
BR  G  04830 
BRG04840 
RRG04850 
BRG04860 
3RG04870 
BRG04880 
BRG04890 
BRG04900 
BRG0491 0 
BRG04920 
3RG04930 
BRG04940 


cn 

-P 


c 

BRG 04950 

DO  101  12  =  1  , NA  0  D  I  M 

BRG04960 

DO  102  J2  =  1  ,  NPO I M 

BRG04970 

DO  103  K2=l ,  NA  4D  I  M 

BRG04930 

C  AVOID  DIVISION  BY  ZERO 

BRG04990 

I  F<  U  (  I2,J2,K2I.EQ.O.)  GO  TO  201 

BRG05000 

I F(A< 12 » J2, K2) . EQ.O. I  GO  TO  103 

BRG0501 0 

C 

BRG05020 

C  CONSTRUCT  XL  FROM  THE  RATIO  OF  THE  PROBABIlTY  MASS  MATRIX  ELEMENTS 

BRG05030 

C  OF 

A  AND  U. 

BRG05U40 

C 

3RG05050 

XL( 12, J2, K2 )=(  A( I 2, J2,K2) /FLOAT ( NA)  ) /( Jl 12, J2  ,K2 )  /FLOAT  (  N'J )  ) 

RRG05060 

GO  TO  103 

BRG05070 

C 

BRG05080 

C  IF 

THE  U  MATRIX  ELEMENT  IS  ZERO, SET  XL  MATRIX  ELEMENT  TO  9999.99 

BRG05090 

C 

BRG05100 

201 

XL ( 12, J2, K2 )=9999.99 

BRG 05110 

103 

CONTI N  JE 

BRG05120 

102 

CONTI N  JE 

BRG05130 

101 

CONTINJE 

BRG05140 

C 

BRG05150 

C  WRITE  XL  MATRIX  INTO  A  STORAGE  FILE 

BRG05160 

C 

BRG05170 

WRI  TE(  8,600)  X'. 

BRG051 80 

600 

FORMAT ( FI  2. 3) 

BRG05190 

STOP 

BPG05200 

END 

BRG05210 

C 

«RG05220 

C 

BRG05230 

C 

BRG05240 

C 

SUBROUTINE  MAPF 

BRG05250 

c 

BRG05260 

SUBRO JT  INE  MAPF(VAR,VARMIN»VARMAX,M1*M2) 

BRG05270 

c 

THIS  SUBROUTINE  MAPS  A  VARIABLE  (VAR)  WHICH  LIES  ON  OR  BETWEEN 

BRG05280 

c 

A  MINIMUM  VALJ5  (VARMIN)  AND  A  MAXIMJM  VALUE  (VARMAX)  TO  AN 

3  RG05290 

c 

ADDRESS  (M2)  OF  AN  ARRAY  0=  LENGTH  (Ml). 

BRG05300 

c 

BRG05310 

I F( VAR . LT . VARMI N  )  V AR=VARMIN 

BRG05320 

IF(VAR. GT. VARMAX )  VAR  =  VARMA  X 

BRG05330 

M2=INT  (  FL0AT(MI  ) *  ( VAR- VAR  MI  Nl / < V ARM AX  - V ARMI N) ) +1 

BRG05340 

IF(M2.3T.M1)  H  2  =  M 1 

BRG05350 

RETURN 

BRG05360 

END 

BRG05370 

non  \j\  oooooo  ooo 


LFILE 

C  THIS  PROGRAM  USES  A  LIKELIHOOD  RATIO  MATRI <  (XL)  to  CREATE  AND 
C  OUTPUT  FILE  THAT  IS  THE  SAME  AS  THE  INPUT  =ILE  EXCEPT  FOR  TH£ 

C.  A4TERM. 

C 

INT£GER*2  NB2 

LOGICAL*!  SOURCE ( 30) /30**  •/ 

DATA  CLR/'LR'^  »  C  S  R / ' SR'/tCSA/'SA'/ 

DIMENSION  XL ( 1 0  » 1 0 , 1 0 ) 

CALL  DSDSET  (1,1600,2,15) 

CALL  DSDSET  (2,1600,2,16) 

READ (9,502)  SOJRCE 
502  F  ORM AT ( 30  A1 ) 

READ! 8, 500)  X. 

500  F0RMATIF12.3) 

NREC  =  1 0  000 

REQUIRED  MAPPING  CONSTANTS 

A  OM I N  =  0  . 

AOMAX  =  1  * 

NAODI M=  10 
PMIN=-l 80. 

PMAX=180. 

NPDI M= 1 0 
A4MI N=0 . 

A  4MA  X  =  1  . 

NA4DI M= 10 

READ  THRU  THE  INPJT  FILE  AND  USE  MAPF  TO  FIND  ADDRESS  OF  XL. 

IF  XL  IS  LT  1  OR  IF  RCD  TAGGED  LR,  THE  BEARING  IS  MADE 
UNACCEPTABLE  BY  SETTING  A 4= 10 . 0  .  IF  GREATER  THAN  1,  THE  BEARING 
IS  ACCEPTABLE  AMD  A4  IS  SET  EQJAL  TO  0.01  . 

DO  100  11  =  1,  NREC 
READ (1,501)  AO, NP,A4,SS,NB2 
N  B  =  NB2 

1  FORMAT! 3  A  4, A2,A2  ) 

P=FLOAT(NP) 

FIND  THE  APPROPRIATE  I,J,K  ADDRESS 

CALL  M4PF(A0,  AOMIN, AOMAX, NAODI M, I ) 

CALL  MA PF( P ,PMI N, PM4X, NPDIH , J ) 

CALL  MAPF (A4, A4M I N, A4MAX , NA4D I M , K ) 

C 

C  COMPARE  THE  XL  MATRIX  ELEMENT  r(ITH  1. 


BRG05400 
BRG05410 
BRG05420 
BRG05430 
BR  G  05440 
BRG05450 
BRG0546  0 
BRG0547  0 
BRG05480 
RRG05490 
BRG05500 
3RG055 10 
BRG0552  0 
BRG05530 
BRG05540 
BR  GO  5550 
BRG05560 
BRG05570 
BRG05580 
BRG05590 
BRG05600 
BRG05610 
BRG05620 
BRG05630 
3  R  GO  5640 
BRG05650 
BRG05660 
BRG0567  0 
BRG05680 
BRG05690 
BRG05700 
BRG05710 
3  RG05720 
BRG05730 
BRG05740 
BRG05750 
BRG05760 
BRG05770 
BRG0578Q 
BRG05790 
3RG05800 
RRG05810 
BRG05820 
BRG05830 
BRG05840 
BRG05850 
BRG05860 
BPG05870 


c 


I F (XL ( I  ,J,K).LT.1.0)  GO  T3  200 
IF(SS.EQ.CLR)  53  TO  200 
C 

C  THIS  PATH,  BEARING  ACCEPTABLE,  SET  A4=0.0l  . 

C 

A4=  0.  01 
GO  T3  201 
C 

C  THIS  PATH,  BEARING  UNACCEPTABLE,  SET  A4=10.0  . 

C 

200  A  4= 1 0 . 0 
C 

C  WRITE  THE  NEW  RECDRD 
C 

201  WRI TE { 2  * 600 )  A0,NP,A4,NB 

603  FORMAT { 4A  4 ) 

100  CONTINJE 
STOP 
END 
C 
C 

C  SUBROUTINE  MAPF 
C 

SUBROUTINE  MA=> F ( V AR , V AR Ml N , VA R MAX , M 1, M2 ) 

I F( VAR . LT.VARMI N )  VAR  =V  AR  MI N 
IF(VAR.GT.VARMAX)  VA  R  =  VARMAX 

M2  =  INT(  FLOAT!  Ml)  *  (  VAR-VARMI  N)  /  (  VARM  AX- J  AP  M I  N  )  )+l 

IF(M2.GT.M1)  M  2  =  Ml 

RETURN 

END 


BRG05880 
BRG05890 
BRG05900 
PRG0591 0 
BRG05920 
BRG05930 
BRG05940 
BRG05950 
RRG05960 
BRG05970 
BRG05980 
BRG05990 
BRG06000 
BRG06010 
BRG0602  0 
BRG06030 
3  RG06040 
BRG06050 
BRG06060 
BRG0607  0 
BRG06080 
SRG06090 
BRG06100 
PRG06110 
RRG06120 
BRG06130 
BRG06140 
BRG06150 
BRG06160 
BRG06170 
BRG06180 


PL0T3 

C  THIS  PROGRAM  PLOTS  THE  RESULTS  OF  ANALYSIS  OF  SWRI  OF  SYSTEM  DATA. 

C 

C 

LOGICA'_*l  SOURCE  (  30)  /30*'  •/ 

DIMENSION  A  ( 1 0 ) ,  S  ( 1 0  ) ,  C  (  i  0 )  ,V(10)»P(10)»RANGE(4),X(6)» 

ADKlOJt  D2(10)«SC1(10)  *S:2(10) 

DATA  01  /100. , 200. ,300. ,400. , 500. t600. , 700. , 800. , 900. ,1000./ , 

A  02  /1000.  ,  2000. ,3000.,  4000. ,5000. ,  6000  .,7000  .  , 

B  8000. ,9000. ,10000. /, 

C  SSI  /2000.,  1000. ,666. ,  500. ,400. ,  333. ,285. ,250. ,222., 

0200./, 

E  SC2  /200. ,  100. ,66. ,50.  ,40. ,33. ,28. ,25.  ,22. ,20./ 

C 

C  IF  MORE  THAN  ONE  DATA  FILE  IS  TO  BE  PLOTTED,  SET  NFILE=1 
C 

NFILE=1 

204  CONTINUE 
C 

C  NUMBER  OF  PLOTS  TO  BE  MADE 
C 

203  RE AD (10, 502, END=2 02, ERR =203)  NPLOT 

502  F  ORMAT ( I  2 ) 

C 

C  FOR  SHORT  (1)  OR  LONG  (10)  SIGNAL  DURATIONS 
C 

N  1  =  1 
C 

C  FIX  THE  SCALE  PARVMETERS 
C 

RANGE! 1 )=1000frFL0AT(Nl) 

RANGE ( 2 )  =  0 . 

C 

C 

DO  100  11=1  , NPLOT 
R  EAD(  10,500)  SOURCE 

500  FORMAT! 30A1) 

READ! 10 ,501 )  A,S,C,Y 

501  FGRMAT(10F7.0) 

I  F ( N 1 .E0.10)  GO  TO  200 
C 

C  DETERMINE  THE  SCA.ING  PARAMETERS  FOR  THE  Y  AXIS  OF  THE:  PLOT 
C 

CALL  MI  NMAXIA,  X<  1  )  ,X(  4)  ,10) 

CALL  MINMAX(S,  XI  2  )  ,X(  5)  ,  10) 

CALL  MINMAX(V»X(3)»X(6),10) 

CALL  MINMAXIX,  RANGE(4),RANGE(3),6) 


BRG09340 
BRG09350 
BRG09360 
BRG09370 
BRG09380 
BRG09390 
BRG09400 
BRG09410 
BRG09420 
BRG09430 
BRG09440 
BRG09450 
BRG09460 
BRG0947  0 
BRG09480 
BRG09490 
BRG09500 
BRG09510 
BRG09520 
B  RG09530 
BRG09540 
BRG09550 
BRG09560 
BRG0957D 
BRG09580 
BRG09590 
BRG09600 
Bp  G0961 0 
BRG09620 
3RG09630 
BRG09640 
BRG09650 
RRG09660 
BRG09670 
9RG09680 
8  RG  09690 
BRG09700 
BRG09710 
BRG0972  0 
BRG09730 
B  RG097  40 
BRG09750 
BRG09760 
BRG09770 
BRG09780 
BRG09790 
BRG09800 
BRG0981 0 


IF(RANGE(4) .Gr .3. )  RANGED  4 ) =0 . 

DO  101  12=1,10 
P(I2)=:< 12)/s:i( 12) 

101  CONTI  NUE 
GO  TO  201 

200  DO  102  13=1 ,10 

102  P(I3)=C( 13) /SC2(  13) 

201  CONTINUE 

W  R I  T  E  (  6 ,600)  SOURCE 

600  FORMAT! «1*, 33(f) , IX, 'SOURCE :  ',30A1,/) 

WRITE! 6 ,601  ) 

601  F  ORMAT ( 9X , ' AV  E  REARING  ERROR  ( OEGREES )  :  X  •  ,/ , 

A9X, *  STO  OF  REARING  ERROR  (DEGREES)  : 

B9  X , • A VE  OF  INTRA-SIGNAL  STO  (DEGREES)  :  .*,//) 

CALL  UTPLTT (01, S, 10, RANGE, l ,1  ) 

CALL  UTPLTT(D1,V, 10, RANGE, 1,2) 

CALL  UTPLTT  (D1,A,  10, RANGE,!.  ,3  ) 

WRI TE(  6 , 603  )  01  , P 

60  3  F  ORMAT (//,1X,' SIGNAL  DURATION:  '  ,3X, 10r  S . 0, / , 2X, ' (MILL  I  SECOND  S  )  •  , 

A/ , IX,* PROBABILI TY  OF  :* ,2X, 10F6.3,/,2X, 'OBTAINING  LOB*) 

WRITE!  6 ,604 )  A,S,V 

604  FORMAT ( / , 13  X,  ' X  : • , 3 X , 1 OF 6 . 0 , / , 1 3 X, *  *  : * , 3X , 1 0F6. 0 , / , 

A13X,*.  :• ,3X,10F6.0) 

REWIND  5 
100  CONTI NJE 

IF(NFILE.EQ.I)  GO  TO  203 

202  CONTINJE 
STOP 
END 

C 

C  THIS  SUBROUTINE  FINDS  THE  MIN  (AMIN)  AND  MAX  (AMAX)  VALUES 
C  OF  AN  ARRAY  OF  LENGTH  N. 

C 

SUBROUTINE  MI  NM  AX  (  ARR  AY  ,  AMIN, AMAX, M) 

DIMENSION  APRAY(N) 

AMAX  =  ARRAY (  1 ) 

AMIN=ARRAY( 1) 

DO  10  D  11=2, N 

IF(ARRAY(I1).GT .AMAX)  A M A X= ARR A Y < 1 1 ) 

IF(ARRAYdl).LT.AMIN)  AMI  N  =  ARR  AY  (II) 

100  CONTINUE 
RETURN 
END 

AMB I  GUI TY 

C  THIS  PROGRAM  EDITS  A  SWRI  FILE  TO  MAKE  THE  A4  TERM  =  10.0  FOR  ALl 
C  RECORDS  THAT  ARE  NITHIN  A  DEFINED  WINDOW  A30JT  THE  90  AND  1B0  DEGREE 
C  AMBIGUITIES  ASSOCIATED  WITH  THE  TRUE  BEARING. 


BRG09820 
BRG09830 
BRG09840 
BRG09850 
RRG09860 
9RG09870 
BRG098P0 
BRG09890 
BRG09900 
BRG09910 
BRG09920 
BRG09930 
BRG09940 
BRG09O50 
BRG09960 
PRG0997  0 
BRG09980 
BRG09990 
RRG10000 
RRG10010 
BRG10020 
BRG10030 
BRG10040 
BRG10050 
PRG10060 
BRG10070 
BRG10080 
BRG10090 
BRG10100 
BRG10110 
BRG10120 
BRG10130 
R  RG 1 0140 
BRG10150 
BRG10I60 
RRGI0170 
BRG101R0 
BRG10190 
B^G 10200 
BRG1021 0 
BRG10220 
BRG10230 
BRG1D240 
BRG10250 
BRG10260 
BBG10270 
BRG10280 
BRG10290 


c 

c 

INTEGER  N  BRG  2 
C 

C  NEED  DSDSET  FOR  CVCMS 
C 

CALL  OSDSETI1, 1600,2, 16) 

CALL  DS0SEK2,  1600,2, 16) 

C 

C  SET  NTBRG  WITH  THE  VALUE  OF  THE  TRUE  BEARING. 

C 

N  T8RG  =  3  37 
NREC=10000 
r 

C  SET  THE  WIDTH  OF  THE  AMBIGUITY  WINDOW.  PROGRAM  ACCOUNTS  FOR 
C  FOLD  AROUND  AT  36D  DEGREES. 

C 

NW!  DTH  =  1 0 
N  A1 80  =  N  TBRG-1 3  D 

IF(NA130.LT.O)  N A 180=NA19 0- 360 
NAP90=NTBRG+90 

I  FCNAP90.GT  .360)  NAP90=NA P90-360 
NAM90=NTBRG-90 

I F ( NAM9  0. LT  .0)  NAM90=360*NAM90 
C 

C  USE  WINDOW  SUBROUTINE  TO  DETERMINE  WINDOW  LIMITS  ABOUT  THE 
C  90  AND  180  DEGREE  AMBIGUITIES. 

C 

CALL  WINDOW(NA180»NW!DTH,Ni,N2 ) 

CALL  WI NDOW (NAD90»NWI  DTH,  N3,N4) 

CALL  WI ND0W(NAM90,NWIDTH, N5 ,N6 ) 

C 

C  LOOP  THRU  ALL  OF  THE  RECORDS  0=  THE  FILE. 

C 

DO  100  11=1  ,NRE  C 
READ! 1,500)  AO , NP  ,  A4 , SS , NBRG2 
500  FORMAT ( 3A4 , 2A2 ) 

NB=NBRS2 

C 

C  IF  BRG  IN  A  RECORD  IS  INSIDE  AMBIGUITY  WINDOW,  SET  A4  VALUE  TO 
C  10.0  AND  WRITE  WH0..5  RECORD  TO  FILE  DSRN  02. 

C 

IF(  (NB.GE.N1.0R.NB.GE.N3.0R.NB.GE.N5) .AMD. 
A(NB.LE.N2.0R.H3.LE.N4  .CJR.N3  .L5  .N6)  )  A4=10.0 
W  RI TE  <  2 , 6  00 )  A0,NP,A4,SS,NBRG2 
600  FORMAT ( 3A4, 2A2 I 
100  CONTI NJE 
STOP 


BRG10300 
BRG10310 
BRG10320 
BRG10330 
BRG10340 
BRG10350 
BRG10360 
BRG10370 
BRG10380 
BRG10390 
3RG10400 
BRG 1 041 0 
BRG10420 
BRG1 0430 
BRG 10440 
BRG10450 
BRG10460 
BRG10470 
BRG10480 
P  RG 1 0490 
BRG10500 
BRG10510 
BRG10520 
BRG10530 
BRG10540 
BRG 10550 
3 RG 105 60 
BRG10570 
BRG10580 
BRG10590 
BR  G 1 0600 
BRG10610 
BRG10620 
BRG 1 0630 
BRG10640 
BRG 1 0650 
BRG10660 
BRG 1 0670 
BRG10680 
BRG10690 
BRG10700 
BRG 1 071 0 
BRG10720 
BRG10730 
BR  G 1 0740 
BRG10750 
BRG10760 
BR  G 1 0770 


noon 


END 


THIS 

THAT 


SU8R3UT I 
DEFINES 


NE  DETERMINES  THE  MIN 
A  WIN33W  OF  WIDTH  NW I 


SUBROUTINE  WIN 3 3 W ( NBR G » NW I3TH , 

NI=NWI DTH/2 

MIN=NBRG-N1 

I F ( MI N. LT . 0  l  MIN  =  MIN  +  360 
MAX=NBRG+N1 

I  F(  MAX.  GT  .360)  MAX=MAX-350 

RETURN 

END 


BEARING  AND 
DTH  CENTERED 

MIN, MAX) 


THE  MAX  BEARING 
AT  NBRG . 


BRG10780 

BRG10793 

BRG10800 

BRGI0810 

BRGI0820 

BRG10830 

BRG1084D 

3RG10850 

BRG10860 

BRG10870 

BRG10380 

BRGI0890 

BRG10900 


BRGCNT 

C  THIS  PROGRAM  SUMS  THE  NJMBER  3F  BEARINGS  B7  INTEGER  BEARING  VALUE 
C  3F  A  DATA  =I_E  0=  10,000  RECORDS.  ONE  CONSTRAINT  CAN  BE  APPLIED 
C  TO  THE  DATA  SET.  THE  A  TERN  3-  S  WR I  DATA  SETS  CAN  BE  SPECIFIED  TO 
C  HAVE  A  MAXIMUM  TOLERABLE  VALJE. 

C 

C  BRGMAX  IS  A  360  ELEMENT  LINEAR  ARRAY  IN  WHICH  EACH  BLOCK 
C  CORRESPONDS  TO  A  NJMBER  OF  OGREES.  THE  CONTENTS  OF  AN  ELEMENT 
C  ARE  THE  NJMBER  OF  RECORDS  WITH  THAT  BEARING. 

r 

DIMENSION  BRGMXI360) 

INTEGER*?  NBR32 
DATA  3RGMX/360*0.0/ 

C  JSE  DSDSE  T  FOR  CP'CMS 

CALL  DSDSETd,  1600,2,  161 
CALL  DSDSET (2, 1 600,2, 16) 

CALL  DSDSET(3, 1600,2, 16) 

CALL  DSDSET (A, l 600,2, 16) 

WRITE! 6,600) 

200  WRITE! 6,601  ) 

RE AO ( 5 , 500) NF I „ E 
WRITE! 6,608) 

R EAD <  5 , 5 D 4 )  NR  EC 
WRITE! 6,602) 

READ! 5, 501 ) A4M A  X 
NDEL=0 
C 

C  READ  THRJ  THE  FILE.  IF  *4  NOT  TOO  LARGE,  AJGMENT  APPROPRIATE 
C  ELEMNET  IN  BRGMX 
C 

DO  101  11=1  ,N  RE  C 

READ! NFI„E,502,ERR=100,END=13D) A4,NBRG2 
NBRG=NBRG2 

I F ( A4.GT. A4MA< ) GO  TO  100 
IF(NBRG.EQ.  0  >  N  3  R  3  =NBR  G«-360 
BRGMX!  NBRG)  =BR3  MX  (NBRG)  M 
GO  TO  101 

100  NDEL  =  N  DEL  +-1 

101  CONTI N JE 
C 

C  WRITE  OUT ?  J  T  TO  FI.E  6 
C 

WRITE! 5, 605  )N: ILE 
WRITE! 6,6  06)44  MAX 
WRITE! 5,603)! I , 3 RGMX ( I ) , I = L ,360) 

WRITE! 5,504)NDEL 
NVAl=NREC-NQEL 


BRG00040 
B  RG00050 
BRG00060 
BRG00070 
BRG00080 
BRG00090 
3RG00100 
BRGOOliO 
BRG00120 
BRG001 3  0 
BRG00140 
BRG001 50 
BRG00160 
BRG00170 
BRG00180 
BRG00190 
BRG00200 
BRG00210 
BRG00220 
B  RG00230 
BRG00240 
BRG00250 
BRG00260 
BRG00270 
3RG00280 
BRG00290 
BRG00300 
BRG0031 0 
BRG00320 
BRG00330 
BRG00340 
B  RG  00350 
3  R  GOO 360 
BRG00370 
BRG00380 
BRG00390 
BRG00400 
BRG00410 
3RG00420 
BR  GOO 430 
BRG00440 
BRG0045  0 
BR  G  00460 
BRG00470 
BRG00480 
BRG00490 
BRG00500 
BRG00510 


c 

C  QUERY  FOR  ANOTHER  RUM.  REWIMD  IF  NECESSARY 
C 

WRITE! 6,607) 

REWIND  NFILE 
READ! 5, 503) MI 
IF(Ml.EQ.l)  CO  TO  200 
STOP 

500  FORMAT (II) 

501  F  ORM AT ! F5 . 3 ) 

502  FORMAT! 8X,A4, 2X, A2) 

503  FORMAT  < 1 1 ) 

504  FORMAT ( I  6 ) 

600  FORMAT! • ORROGRAM  SUMS  INDIVIDUAL  I  NT  5  3  =  R  BRGS  OF  A  DATA  SET') 

601  FORMAT! • OENTER  PI LE  NJM3ER  ill)') 

602  FORMAT! • OENTER  MAX  A4  VALUE  (F5.3)') 

603  FORMAT! 10! 1X»I3,1X,F5.1.2X)/) 

604  FORMAT!  •  0NUMI3ER  OF  DELETES  =  ’,15) 

605  FORMAT! ' OFILE  NUM8ER  *,I1) 

606  FORMAT ! ' 0 A4  MAX  =  *,F5.3) 

607  FORMAT!  '  OTYPE  A  1  FOR  ANOTHER  RUN,  TYRE  A  DIFFERENT  INT  TO  END') 

o08  FORMAT!  •  OENTER  NUMBER  OF  RECORDS  !I6)M 

END 


BRG0052D 
BRG00530 
3RG00540 
BRG00550 
BRG00560 
BRG0057  0 
BRG00580 
RRG00590 
3RG00600 
BRG00610 
BRG00620 
BRG00630 
BRG00640 
BRG0065  0 
BRG00660 
BRG00670 
PRG00680 
BRG00690 
BRG00700 
BRG00710 
BRG00720 
BRG00730 
BRG00740 


BRGH2 


C 

C 

c 

c 

c 

c 


c 

c 

c 


THIS  PROGRAM  DISP.AYS  A  HISTOGRAM  OF  FILTERED  BEARINGS  IN  FILE  DSRN  1 
FILTERING  IS  DONE  ON  THE  A4  TERM.  MAXIMUM  LEVEL  IS  SET  IN  A4MAX. 

A 4  ADMISSIBLE  BRGS  ARE  USED  TO  DETERMINE  MEAN  AND  STD.  A  WINDOW 
IS  CREATED  ABOJT  THIS  m=AN.  WIDTH  OF  WINDOW  IS  STD  TIMES  AN  INTEGER 
ENTERED  6V  THE  USER 

LOGICAL’!'!  SOURCE!  30)  fib*'  */ 

INTEGERS  NBRG2 

DIMENSION  BR3MATI10DDD),SnT(5),SCALE(2»,I0PT(5),  DATMAT  (  1 0000  ) 
DATA  IOPT/1  ,0,3,0,!/ 

JSE  DSDSET  FOR  CP  V  E  MS 


CALL  DSDSET  (1,  1600,2,  16) 

WRITE! 10, 600) 

WP.I  TE ( 1 0, 603) 

READ!  5,  5D  2  )  SOJRCE 
201  WRITE(10,607) 

R  EAD!  5,  500)  A 4 M  A X 
WRI  TE!  10,  610) 

READ ! 5 , 500)  X2 
N  VAL  =  0 
C 

C  READ  10000  RECORDS  AND  FILTER  ON  A4.  FILL  3RGMAT  WITH 
C  A 4  ADMISSIBLE  BEARING  VALUES. 

C 

DO  100  11=1,10000 

READ! 1, 501 , ENO=200)  A4,NBR32 

NBRG=NBRG2 

IF! A4.3T.A4MAO  30  TO  100 
NVAL=NVAL*1 
BRGMAT ( NVAL ) =NBRG 
CONTI NJE 
CUNTINJE 


100 
200 
C 


C 

c 


c 

c 

c 


FIND  MEAN  AND  STD  OF  BRGMAT  AND  SET  PARAMETERS  OF  WINDOW. 

CALL  3EIUGR!BRGMAT,NVAL, 103T,STAT , IER) 

X1  =  FL3AT( I  FIX! STAT! 1 )  ) ) 

X3=  FL 0 AT!  I F I X(  X2 * SQR T ! STA T ( 5)  )  )  ) 

X4=XlfX3 
X  5=X1 - X3 
M  1  =  0 

FILTER  BRGMAT  T HR J  WINDOW.  PJT  WINDOW  ADMISSIBLE  VALJES  IN  D ATAM AT . 


BRG00770 
BRG00780 
.3RG00790 
BRG00800 
BRG00810 
RRG00820 
BRG00830 
BRG00840 
BR  00085  0 
BRG00860 
BRG00670 
BRG00880 
BRG00890 
BRG00900 
BRG00910 
BRG00920 
BRG00930 
BRG00940 
3RG00950 
BRG00960 
3RG00970 
BRG00980 
RRG00990 
BRG01000 
BRG01010 
BRG01020 
8RG01030 
BRG0I040 

3  R  GO  10  50 
BRG01060 
BRG01070 
BRG01080 
BRG0L090 
BR  GO  1 1 0  0 
BRG01 11 0 
BRG01120 
BRG01130 
BRG01140 
BRG01150 
BRG01160 
BRG011/0 
BRG01180 
BRG01190 
BRG01200 
BR GO  1210 
BRG01220 
BRG01230 
BRG01240 


DO  101  12*1,  M/M 
X6=BR3MAT(  12)  +  3  60  . 

IF!  (BR3MATI  I2).LT.X5.0R.BRSMAT(  I2KGT.X4)  .AND.X6.GT.X4l  GO  TO  101 
M1  =  M1H 

I F ( BR3M AT ( I  2 ) .  L  T  .  X5  )  BRGM AT ( 12) =  BRGMA  T( 12) +36  0 
D ATMAT ( Ml ) =8RGMAT ( 12) 

101  CONTINUE 

WRITE<6,605)  COJRCE 
W  RI  TE  (  5  »  6  08  )  A'tMAX 
C 

C  CREATE  AND  OUTPUT  HISTOGRAM.  SET  JP  FOR  ANOTHER  RUN. 

C 

CALL  HI STF! DATMAT  ,M1 , 0) 

REWIND  1 
WRITE! 10,609) 

R  EADI  5,  505  )  M2 

I  FI  M2 .EQ. 1 )  SO  TO  201 

STOP 

500  FORMAT  I F5 . 3 ) 

501  FORMAT (8X,A4,2X»A2) 

502  FORMAT ( 30A1  ) 

504  F0RMATIF5.1) 

505  FORMAT!  ID 

600  FORMAT!  'OTHIS  PROGRAM  DISPLAYS  A  HISTOGRAM  OF  BEARINGS  IN  FILE*, 

A  *  DSRN  01  * ) 

603  FORMAT! • OENTER  RADIO WAVE  SOURCE  (3041)') 

605  FORMAT! '1RADI0WAVE  SOJRCE:  *,30A1) 

607  FORMAT! • OENTER  MAXIMUM  A4  VALUE  (F5.3)') 

608  FORMAT! //,' OA'*M  AX  :  *,F5.3J 

609  FORMAT! • OENTER  A  1  FOR  ANOTHER  RUN, IF  NOT,  DIFFERENT  INTEGER*) 

610  FORMAT!  • OENTER  PERMISSIBLE  NJMBER  OF  STD  (F5.3)*) 

END 


BRG01250 
BRG012S0 
BRG01270 
BRG01280 
**RG01290 
BRG01300 
BRG0131 0 
BRG0132  0 
BRG01330 
3RG01340 
BRG01350 
BRG01360 
BRG01370 
BRG01380 
3RG0I390 
BRG01400 
BRG01410 
3RG01420 
BR GO 1430 
BRG01440 
BRG01450 
BRG01460 
3RG01470 
BRG01480 
BRG01490 
BRG01500 
BRG01510 
BRG01520 
BRG01530 
BRG01540 
9RG01550 
BRG01560 


BRGHIS 


C  THIS  PR0GR4M  DISP.AYS  A  HISTOGRAM  OF  THE  BEARINGS  IN  FILE  DSRN  1. 

C  FILTERING  IS  DONE  3  M  EACH  R  EG  OR  D  ON  THE  A \  tfrm. 

C  SCALING  OF  THE  HISTOGRAM  CAN  BE  SET  BY  THE  JSER. 

C 

LOGICAL*!  SOURCE! 30) /30**  •/ 

D I  MENS  I  ON  BRGMAT!  10000) , SCALE! 2  ) 

C  USE  DSDSET  FOR  CP/CMS 

CALL  DSD$ET(1, 1600,2, 16) 

WRI  TE(  L  0,  6001 
WRITE! 10, 603) 

READ ( 5 »  592 )  SDJRCE 
WRITE!  LO,  607) 

READ! 5«  500)  A  V  M  A  X 
WRITE (10, 601) 

READ! 5,  504 )  SGALE(l) 

WRITE! 10,602) 

READ! 5 »  50  4 )  SG ALE ( 2) 

NVALO 

C 

C  READ  RECORDS,  FILTER  OUT  HIGH  A4  VALUES.  5'JT  A4  ADMISSIBLE 
C  VALUES  IN  LINEAR  ARRAY  BRGMAT. 

C 

DO  100  11=1,10000 
READ! 1» 501, EN)=200)  A4, N3RG 
IF! A4.GT.A4MAX)  GO  TO  100 
NVAL=NV AL+1 
BRGMAT! MVAL )=NBRG 
100  CONTI N JE 
200  CONTINUE 

WRI TE ( 5 , 605 )  SOURCE 
C 

C  DETERMINE  HISTOGRAM  SCALE  AND  THEN  CREATE  IT  WITH  HISTF. 

C 

CALL  FIX!  SCALE  ) 

CALL  HI STF! BRGMAT ,NVAL, 0) 

STOP 

500  FORMAT ! F5 .3  ) 

501  FORMAT!  8X,2A4) 

502  FORMAT { 30A1  ) 

504  FORMAT! F5.1  ) 

600  FORMAT!  • OTHIS  R R 3 GRAM  DISPLAYS  A  HISTOGRAM  OF  BEARINGS  IN  FILE*, 
A'  DSRN  01  *  ) 

601  FORMAT (' OENTER  MINIMUM  VALJE  OF  SCALE  (F5.1)‘) 

602  FORMAT! * OENTER  MAXIMUM  VALJE  OF  SCALE  (F5.1)*) 

603  FORMAT!  ' OENTER  RADIOWAVE  SOURCE  (30A1)*) 

605  FORMAT! /////,* DRADIOWAVE  SOURCE:  *,30Al) 


BRG01590 
BRG01600 
BRG01610 
3RG01620 
BR  GO  1 63  0 
BR  GO  1640 
BRG01650 
BRG01660 
3RG01670 
RR  GO 168  0 
BRG01690 
BRG01700 
BRG01710 
BRG01720 
BRG01730 
BRG01740 
BRG01750 
BRG01760 
BRG01770 
8PG017B0 
BRG01790 
BRG01800 
BRG0181 0 
BRG01820 
3RG01830 
BRG01840 
BRG01850 
BRG01860 
BRG01870 
BRG01880 
BRG01890 
BDG01900 
BRG01910 
RRG01920 
BRG01930 
BRG01940 
BRG01950 
BR  GO  1 960 
3RG01970 
BRG01980 
BRG01990 
BRG02000 
BRG02010 
PRG02020 
BRG02030 
BRG02040 
BRG02050 
BRG02060 


607 


FORMA  T(  'DENTED  M4XIMJM  44  7  4  L  'J  5  (F5.3)') 
END 


BRGD2070 

BRG02080 


noon  non  ro  oooo  oooo 


T I  ME  1 

THIS  PROGRAM  PLOTS  THE  43  ,  A4 , 3=  ARI  NG  AND  MASE  ME AS'JREMENTS  FROM 
THE  SWRI  SPACED  L33P  ANTENNA  VERSJS  TIME. 


LOGIC 
AC  B » CP 
INTE3 
DATA 
ACB/' B 
BC  SA/ ' 
CALL 

FOLLOWING 

THE  BEARI 

A  OMIN 
AOMAX 
A4MIN 
A4MAX 
NPMIN 
NPMAX 
MB  MI  M 
N6MAX 
WRITE 
WRITE 
READ! 

6  WRITE 
READ! 

SET  ARR AT 


4L 

*1  sour: 

“ 

(20) 

, 3R3ARY ( 45 ) 

, PHSARYI 45 

) 

, AOAR  Y( 

20) , A4ARY (20 

t  O 

A3  ,CA4,  C 

S 

, CINVAL 

ER 

*2  NBRG2 

BR 

GARY/45* 

« 

’/, 

PHSARY/45*' 

•  /, AOARY/ 

20*'  '/, 

A4ARY/2 

0*'  * 

*  / 

, CP/ *P» / 

CAO/ 

'  3*  /, 

C  A  4/ ’ 4 

*/,CS/'  •/ 

f 

Cl R/'LR 

' / , CSR / 

•SR'  / 

SA 

•  / , C IN/A 

L 

/  t  ?, 

/ 

DS 

DSETU,  1 

6 

30,2 

,16) 

0 

ECLARATI 

ONS  R 

EPRES 

ENT  TH 

E  RANGE  OF 

VALUES 

EXPECTE 

D  FOP 

NG 

, PHASE,  A 

0 

AND 

A  4  T 

S  R  MS  . 

=  0.0 
=  2.0 
=  0.0 
=  1  .0 
=  -180 
=  180 
=  0 

=  3  60 

(1 0,600) 

(  10,  603) 

5,531)  S3JRCE 
( 10,604) 

5,  532)  NO  J  T 

3  ARAMETE  R  S  Ml, M2  AND  M3  FOR  EITHER  PRT  OR  CON  OUTPUT 


M  1=45 
M2  =  23 
M3  =  20 

I  F  <  NO  LIT • E3 . 1 )  SO  TO  23B 
M  1  =  24 
M  2=1  5 
M3  =  10 


0ET5RMINF  NUMBER  OF  RECORDS  TO  BE  PLOTTED  IN  TERMS  OF  TIME. 
ONE  SECOND  EQUALS  =  IFTY  RECORDS. 

208  WRITEI10.601) 

R  5A0(  5,  530 )  NSTART 
WRITE(10,602) 

R  EAD(  5 ,  500)  NSTDP 


BRG02110 
RRG02120 
BRG02130 
BPG02140 
BRG021 50 
BR302160 
BRG02170 
BRG02180 
BRG02190 
RRG02200 
BRG02210 
BRG0222  0 
BRG02230 
BRG02240 
BRG02250 
BRG02260 
BRG02270 
BRG02280 
BRG02290 
BRG02300 
BRG02310 
BRG0^320 
RRG02330 
BRG02340 
BRG02350 
BRG023  60 
BRG02370 
BRG02380 
RRG02390 
BRG02400 
BRG0241 0 
BRG02420 
BRG02430 
BRG02440 
BRG02450 
BRG02460 
BRG02470 
BRG02480 
BPG02490 
BRG02500 
BRG02510 
BRG02520 
BRG02530 
BRG02540 
BRG02550 
RRG02560 
BRG02570 
BRG02580 


noon  ooooo  ooo  non  onnivM  ooo 


NADV=0 

! F ( NST  A  RT . EQ. ) ) 

N ADV=  N  S  TART  *53 

READ  FILE  TO  START  POINT . 

DO  10  0  Il=lj»NA3v/ 

00  READ<  1,  503, ENO  =  205,ERR  =  234» 
00  NHALT  =NST3P*5D 


DETE 

RMINE 

NUMBER  0= 

RECO 

NRUN  = 

NHAL 

T-NAD  V 

I  F  ( N  R 

JN.L 

E.O)  S 

0  TO 

PRIN 

IT  THE 

OUT 

PUT  LA 

BEL. 

CALL 

LABEL! NO J  T 

,  SOUR 

READ 

i  THRU 

THE 

REOJE 

STED 

DO  101  12 

=  1  ,NRJ 

N 

READ! 

1 , 50 

3 , ENO= 

101, E 

NBRG  = 

NBRS 

2 

USE 

MAPF 

AND 

mapi  r 

0  MAP 

AOARY,  PH 

S  ARY 

,  A4AR 

Y  AND 

MAPI 

IS  F 

OR  I 

NT  EGER 

MAPP 

CALL 

MAPF 

(AO, AO 

MIN,  A 

AOARY 

(  MAI 

)  =CAD 

CALL 

MAPI 

( NPHS, 

N°  MI  N 

PHSARYI MA 

2)  =  C5 

CALL 

MAPF 

(  A4,  A 4 

MIN,  A 

A4ARY 

(  MA3 

i  =  C  A  V 

CALL 

M&PI 

( N  RRS , 

NBMIN 

BRGAR 

Y I  MA 

4>  =C3 

IFISS 

.EQ. 

ClR.OR 

.  SS.E 

RECORDS  TO  BE  READ, 


SO  TO  210 

IF  LR,  SR  OR  SA  FLAGS  ARE  READ  IN  A  RECORD,  rfRITE  THE 
?  TO  THE  FIRST  ELEMENT  OF  OJTPJT  ARRAYS. 

209  AOARY  I L )=  CIN/AL 
PHSAR Y ( 1)=  CINVAL 
A4AR Y ( 1  )  =  CINYAl 
BRGARY(1)=  Cl  NY  AL 

210  I F( NO JT .EO . 1 )  SO  TO  202 


ARRAYS 
POINT  AND 


SYMBOL 


BRS02590 
BRG02600 
PRG02610 
BRG02620 
BRG02630 
B  RG02640 
BRG02650 
BRG02660 
BRG02670 
BRG02680 
BRG02690 
BRG02700 
BRG0271 0 
RRG0272  0 
BRG02730 
BRG02740 
BRG02750 
BRG02760 
BRG02770 
BRG02780 
BRG02790 
BRG02800 
BRG02810 
BRG02820 
BRG0283  0 
BRG02840 
BRG02850 
BRG02860 
BRG02870 
BRG02880 
BRG02890 
BRG0290O 
BRG02910 
BRG0292  0 
BRG02930 
PRG02940 
BRG02950 
BRG02960 
BRG02970 
BRG02980 
BRG02990 
BRG03000 
BRG03010 
BRG03020 
BRG0303  0 
RRG03040 
BRG03050 
BRG03060 


c 

C  THIS  PATH,  WRITE  ro  CRT. 

C 

WRITE!  3,605  )(  (  3  RGARY  ( I  )  ,  I  =L  ,  2  4  )  ,  I  PHS  A  R  Y<  I  )  ,  1=  1 , 24)  , 

A!  AOARY!  II  *  I  =1,15)  ,(  A4A*Y(  I  )  ,1  *1,  101) 

GO  TO  203 
C 

C  THIS  PATH,  WRITE  TO  PRINTER. 

C 

202  WRITE (8, 606)  3 R 3 AR Y , P HS A R Y , A  0 A R Y , A 4A R Y 
20  3  AOARY ( H  A1 )=  CS 

PHSARY! MA2 ) =C$ 

A4AR  Y  (  MA3  )  =  CS 
BRGARY  (  MA4)  =CS 

101  CONTINUE 
C 

C  QUERY  FOR  ANOTHER  RUN. 

C 

207  WRI  TEC  1  0,  607) 

READ!  5,  502)  NR E RUN 
IF(N.RERUN.NE.l)  GO  TO  204 
NCL0Sr  =  10000-\HALT 
C 

C  IF  ANOTHER  RJN,  READ  TO  END  OF  FILE  AND  REWIND. 

C 

DO  102  I3  =  l,Ni:_3SE 

102  READ! I, 503, ENO=205) 

205  REWIND  1 

GO  TO  206 
204  STOP 

600  FORMAT!  ' OTH  I S  3  3  OGRA M  GRAPHS  SWRI  DATA  VERSUS  TINE') 

601  FORMAT! 'OFILE  CONSISTS  0=  10000  RECORDS  REPRESENTING  200  SECS.', 
4/,'OENTER  START  TIME  IN  SECONDS  <I3):*) 

602  FORMAT (' OENTER  STO»  TIME  IN  SECONDS  (13)  (200  MAX):') 

603  FORMAT ( 'OENTER  SOURCE, FREQ  AND  TIME:') 

604  FORMAT!  •  OOUTPJT  TO  PRINTER  ENTER  1,  TO  CONSOLE,  ENTER  2:') 

60  5  FORMAT! 1X.24A1,  • **' ,24Al, , 15AI, •**• , 10A1) 

606  FORMAT! IX, 4541, '** ,45 Al, •*'  ,20A1,20A1) 

607  FORMAT!  « OENTER  1  FOR  ANOTHER  RUN, IF  NOT,  ENTER  OTHER  INTEGER: ') 

500  FORMAT (13) 

501  FORMAT! 20A1) 

502  FORMAT!  ID 

503  FORMAT! 3A4, 2A? ) 

END 

C 

C  THIS  SUBROUTINE  WRITES  A  LABEL  TO  A  CRT  OR  A  PRINTER. 

C 

SUBROUTINE  LA3 EL ( NOJT, SOJRC E) 


BRG03070 
BRG03080 
BRG03090 
BRG03100 
BRG031I0 
BRG03120 
3RG03130 
BRG03140 
BRG03150 
BRG03160 
BRG03170 
BRG03180 
BRG03190 
8RG03200 
BRG032I 0 
BRG03220 
BRG03230 
3RG03240 
BRG0325  0 
BRG03260 
BRG03270 
BRG03230 
3  R  G0^290 
BRG03300 
BRG0331 0 
BRG03320 
BR  00333  0 
BRGD334D 
BRG03350 
BRG03360 
BRG03370 
BRG03380 
BRG03390 
BRG03400 
BRG0341 0 
BRG03420 
BRG03430 
BRG03440 
BRG03450 
BRG034b0 
BR003470 
BRG03480 
BRG03490 
BRG03500 
BRG035]  0 
BRG03520 
BRG0353  0 
BRG03540 


logical*i  soir:e<20) 

WRITE!  3 ,600) 

WRI  TE  (  8 , 601  )  SOURCE 
I F( NO JT .E3. 1)  SO  TO  200 
W  RI  TE (  8 , 6  02  ) 

WRI  T E(  8,603  ) 

RETUR  N 

200  WRITE! 3,604) 

WRITE! 3,605  ) 

RETURN] 

600  FORMAT! ' 1'  , 20X, ' SWRI  SPACED  LOOP  ANT  EMMA  PARAMETERS',/, 

A21X, • AO, PHASE, A4  AND  BEARING  VERSUS  TIME*,////) 

601  FORMAT! 21X, 'RADIOWAVE  SOURCE:  *,20Al,//) 

602  FORMAT! L 1 X , • BE AR I NG'  , 19X , • PH AS E • , 18X , • AO ' , 1 IX , • A4 • , / , 

A*  00*  ,20X,  *3  60*  ,  2X,  '-130'  ,  7X,  •  00*  ,  7X,*  U80*  ,  2X,'  O'  , 

B 1 3X  , ' 2 '  , 2  X ,  '  0 ' , 3  X  , '  1 '  ) 

60  3  FORMAT!  1X,7 9! •«•)♦/, IX, 79! •*•) ) 

604  FORMAT! 13X, 'BEARING'  ,41X,  • PHASE* , 33X, • AO' , 19X,' A4* ,/, 

A'  00'  »41X,  '360',  1  X  ,  *  - 1 80  '  ,  3  7  X,  '  +180'  ,  'O'  ,  1 8X  ,  '  2'  , 

BIX, 'O',  18 X ,  '1*  ) 

605  FORMAT!  IX,  132!  •  *’  )  ) 

END 

C 

C  THIS  SUBROUTINE  MA3S  A  VARIABLE  (VAR)  WHICH  LIES  ON  OR  BETWEEN 
C  A  MINIMUM  VA„!JE  (VARMIN)  AND  A  MAXIMUM  VA.JE  (VARMAX)  TO  AN 
C  ELEMENT  !M2)  OF  AN  ARRAY  OF  LENGTH  (Ml). 

C 

SUBROUTINE  MA3  =  ( V AR , V AR MI N , V A RMAX , Ml ,  M2 ) 

IF(VAR.GT. VARMAX )  VAR= VARMAX 

M2=I  NT!  EL  OAT  (Ml  )  #  (  VAR-VARMI  INI)  /  (  VARMAX- V  ARM  I  N)  )  +1 
I F ! M2 • ST • Ml )  M  2  =  M  1 
RETURN 
END 
C 

C  THIS  SUBROUTINE  IS  AN  INTEGER  VERSION  OF  the  ABOVE  ROUTINE  MAPF. 
C 

SUBROUTINE  MAPI ( N , NM I N, NMAX , Ml , M2 ) 

I F! N. ST .NMAX)  N  =  NMAX 

M2=INT( FLOAT! Ml J AFLOAT! N-NMIN) /FLOAT  (  NMAX -N MI  M)  )  *1 

IF(M2.ST.M1)  M  2  =  M 1 

RETURN 

END 


BRG0355  0 
BRG03560 
BRG03570 
3RG03580 
BRG03590 
BRG03600 
BRG03610 
BRG03620 
BRG03630 
BRG03640 
B  R  GO 3650 
BRG03660 
BRG03670 
BRG03680 
BRG03690 
BRG03700 
BRG03710 
BRG03720 
B  RG  03730 
BRG03740 
BRG03750 
BRG03760 
BRG03770 
BRG03780 
BRG03790 
BRG03800 
BRG0381 0 
BRG03820 
BRG03830 
BRG03840 
BRG03850 
BRG03860 
BRG03870 
BRG03880 
BRG03890 
BRG03900 
BRG0391 0 
BRG03920 
RRG03930 
BRG03940 
3RG03950 
B  RG  03960 
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